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ABSTRACT 


The  objective  of  this  research  is  to  investigate  the  dynamic  elastic  properties  of  materials  and 
their  temperature  and  frequency  dependence  using  an  acoustic  resonance  based  method.  In  this 
technique,  die  torsional,  flexural,  and/or  longitudinal  vibrational  modes  of  a  "free-free"  bar  are 
selectively  excited  and  tracked  as  a  function  of  temperature.  The  resonance  frequency  of  the  torsional 
mode  is  used  to  obtain  the  dynamic  shear  modulus.  The  dynamic  Young's  modulus  is  obtained  from 
either  the  resonance  frequency  of  the  flexural  or  longitudinal  mode  of  the  bar.  The  quality  factor, 
Q,  of  each  mode  is  measured  to  obtain  the  damping  properties  of  the  material.  A  two  channel  phase- 
locked  loop  (PLL)  is  used  to  track  the  resonance  frequency  as  a  function  of  changing  temperature  for 
the  particular  resonance  mode  selected.  Using  this  technique,  the  storage  modulus  and  loss  tangent 
may  be  obtained  in  a  continuous  fashion.  Materials  tested  in  this  thesis  include:  Polyurethane  PR- 
1592,  a  common  sonar  encapsulant.  Polymethyl  Methacrylate  (PMMA)  or  plexiglass,  and 
Polycarbonate.  This  research  encompasses  the  theory,  accuracy,  limitations,  and  applications  of  this 
measurement  technique. 
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Z.  INTRODUCTION 


A.  BACKGROUND  AMD  MOTIVATION 

Elastomeric  materials  are  important  in  many  engineering 
applications.  This  research,  in  part,  was  undertaken  to 
support  continuing  efforts  to  develop  fiber  optic  acoustic 
sensors.  The  design  and  fabrication  of  fiber  optic 
hydrophones  is  particularly  dependent  upon  the  mechanical 
properties  of  the  materials  used  in  the  hydrophone.  In 
order  to  enhance  the  sensitivity  of  fiber  optic  devices, 
various  compliant  transducers  have  been  designed  to  generate 
large  amounts  of  strain  or  phase  shift  within  the  arms  of  an 
optical  fiber  interferometer  [Ref.  1,2].  The  transducers 
have  been  produced  in  the  shapes  of  plates,  mandrels  and 
ellipsoids.  These  techniques  often  involve  bonding  the 
fibers  to  an  elastomeric  material  in  such  a  way  as  to  induce 
strain  within  the  fiber  itself.  More  common  interest  in  the 
acoustic  properties  of  materials  originates  from  the 
automotive  and  industrial  communities  in  which  the 
elastomers  play  an  important  role  in  the  damping  of  noise, 
shock  and  vibrations. 

It  is  the  static  elastic  moduli  that  determines  the 
deformation  of  a  body  to  a  constant  force.  For  example,  in 
a  hydrophone  the  pressure  at  a  particular  operating  depth 
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products  a  force  on  the  transducer  that  results  in  a 
corresponding  deformation.  However,  the  dynamic  moduli  are 
used  to  determine  the  acoustic  sensitivity  or  responsivity 
of  a  transducer  to  time-varying  loads. 

The  technique  used  here  for  measuring  the  dynamic 
Young's  and  shear  modulus  is  a  refinement  of  one  first 
developed  by  Barone  and  Giacomini  [Ref.  3]  with 
modifications  made  by  Rudnik  at  UCLA.  Additional  revisions 
to  the  technique  and  apparatus  have  been  made  by  Garrett  and 
Brown  at  the  Navy  Postgraduate  School  [Ref.  4,5].  This 
method  was  recently  used  by  Wetterskog,  et.  al,  [Ref.  6]  to 
measure  the  temperature  dependence  and  dynamic  moduli  of 
materials  considered  for  use  in  their  fiber-optic 
hydrophones.  The  technique  was  also  used  by  Tan  [Ref.  7] 
while  investigating  the  viscoelastic  transitions  of 
elastomeric  materials. 

B.  METHOD  07  MEASUREMENT 

The  dynamic  moduli  of  materials  are  often  measured  with 
forced  vibration  techniques  which  are  generally  classified 
as  either  resonant  or  non-resonant  procedures.  Here,  a 
resonance-based  method  employing  a  freely  supported  rod  is 
used.  The  resonant  frequency  is  proportional  to  the  square 
root  of  the  modulus  of  the  material.  In  general,  the 
displacement  of  a  sample  at  resonance  is  larger  by  a  factor 
equal  to  the  quality  factor  (Q)  than  it  is  under  a  static 
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load.  Accordingly,  when  the  Q>1,  resonant  techniques  have  a 
higher  signal-to-noise  ratio  than  nonreeonant  techniques. 

The  (quality  factor  is  the  ratio  of  the  energy  stored  to  the 
energy  lost  within  a  material  during  one  cycle  and  is 
inversely  proportional  to  the  loss  tangent.  The  loss 
tangent,  tan  6,  is  defined  as  the  tangent  of  the  phase  angle 
between  the  applied  stress  and  resulting  strain.  The  Q, 
therefore,  is  a  measure  of  the  absorptivity  of  the  material. 
The  complex  dynamic  moduli  (storage  and  loss  properties  of 
the  material)  can  be  determined  using  the  resonant  frequency 
of  each  mode,  the  physical  dimensions  and  the  mass  of  the 
sample  rod. 

The  method  discussed  here  employs  a  forced  vibration 
created  by  electrodynamic  transduction.  The  force  is 
generated  on  a  small  wire  transducer  coil  mounted  on  a 
sample  rod  placed  between  a  pair  of  permanent  magnets.  The 
detection  of  the  vibration  at  the  receiver  is  also  by  means 
of  electrodynamic  transduction.  A  receiver  coil  is  mounted 
on  the  opposite  end  of  the  sample,  and  located  between  an 
identical  pair  of  permanent  magnets .  An  important  feature 
of  this  technique  is  that  the  torsional,  flexural,  or 
longitudinal  mode  of  the  sample  rod  can  be  selectively 
excited  depending  upon  the  orientation  of  the  transducer 
coil  within  the  magnetic  field.  Thus,  both  the  shear  and 
Young's  modulus  can  be  determined  independently  using  this 
single  apparatus .  The  resonant  technique  utilized  in  this 
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investigation  is  discussed  and  illustrated  in  sore  detail  in 
References  4  and  5  as  well  as  Chapter  III  of  this  Thesis. 

C.  SCOTS 

The  scope  of  this  project  is  to  further  refine  this 
aeasureaent  technique  to  pernit  independent  neasurenent  of  Q 
as  a  function  of  teaperature,  and  develop  new  aethods  to 
test  very  coapliant  and  hiqh  loss  naterials.  Further,  the 
capabilities  and  liaitations  of  this  method  as  to 
restrictions  on  damping  and  stiffness  properties  of 
candidate  materials  were  to  be  investigated.  Materials 
tested  in  this  thesis  include;  Polyurethane  PR-1592,  PMMA, 
and  Polycarbonate. 
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a.  nmtopocTicM 


1.  Clauioal  Theory 

The  classical  thaory  of  alas tie! ty  describes  the 
behavior  of  Materials  categorically  as  either  ideal  elastic 
solids  or  viscous  fluids.  The  Mechanical  properties  of  the 
elastic  solid  are  described  by  Hooke's  law:  For  small 
def orMations ,  stress  is  directly  proportional  to  strain  and 
independent  of  the  rate  of  strain.  Algebraically  that  is. 


o  «  e,  o  =  G  e  ;  (1) 

where  c  is  the  applied  stress,  «  is  the  resulting  strain, 
and  G  represents  the  appropriate  modulus  for  the  system  (in 
this  case  G  is  the  shear  modulus) . 

This  implies  that  the  moduli  of  elastic  materials 
have  no  frequency  dependence.  Ideal  viscous  fluids  are 
governed  by  Newton's  viscosity  law:  The  applied  stress  is 
proportional  to  the  rate  of  strain  and  independent  of  the 
strain  itself.  Algebraically  that  is. 


dV. 


dV 


(2) 


where  o  is  the  applied  stress,  dV/dn  is  the  rate  of  strain 
and  n  is  the  coefficient  of  viscosity.  The  energy 
transmitted  to  a  viscous  liquid  is  dispersed  as  heat. 


5 


Contrastingly,  in  an  ideal  elastic  solid  the  energy  of 
deformation  is  entirely  recoverable.  [Ref.  8,9,10] 

Under  certain  stress  and  temperature  conditions  many 
materials  exhibit  a  combination  of  elastic  and  viscous 
behavior.  This  suggests  that  materials  may  exhibit  either 
ideal  Newtonian  viscosity  or  ideal  elasticity  properties  in 
certain  temperature  regions.  These  materials  are  classified 
as  being  viscoelastic.  The  mechanical  properties  of  these 
materials  are  specified  by  a  frequency  and  temperature 
dependence  of  the  complex  elastic  moduli  and  loss  modulus. 

2.  Response  of  a  Driven  Famonic  Oscillator 

If  a  harmonic  oscillator  is  driven  with  a  sinusoidal 
force  such  that, 

F(t)  =F  sin  (cot)  ,  (3) 

then  the  displacement  for  the  system  is  of  the  form, 

x(t)  »  X  sin («£:-♦)  .  (4) 

0  is  the  phase  difference  with  respect  to  the  driving  force. 
The  amplitude  of  displacement  (X)  in  terms  of  the  resonance 
frequency  and  the  quality  factor  can  be  expressed  as. 


where  f,*»f/f0  (f0  is  the  system  resonance  frequency.  The 
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dimensionless  normalized  displacement  response  of  the  system 
is  therefore  expressed  as, 


The  phase  difference  between  the  driving  force  and 
displacement  is, 

♦  y°,)  . 


(«) 


Equations  6  and  7  are  plotted  verses  the  frequency  ratio  for 
various  Q's  in  Figures  2.1  (a)  and  (b) . 
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Phase  Angle  (degrees)  Normalised  Displa 


(a) 


(b) 

Figure  2.1  (a)  Normalized  displacement  response  vs. 

f*,  for  Q's  of  10,  5,  2,  1  and  .5. 

(b)  Phase  (displacement)  vs.  f*,. 
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Note  that  the  maximum  displacement  of  che  system 
does  not  necessarily  correspond  with  the  resonance  frequency 
(f„,»l) .  This  is  especially  true  for  lov  Q  systems  (Q<  5) . 
However,  it  is  seen  from  figure  2.1(b)  that  the  maximum 
rate-of -change  of  the  phase  always  occurs  at  resonance, 
which  corresponds  to  0=90°  for  the  displacement  response. 

The  corresponding  velocity  of  the  system  is  obtained 
by  taking  the  derivative  of  the  displacement.  The 
amplitude  of  the  velocity  response  (v)  is, 

v  «  <■> 

The  normalized  amplitude  of  the  velocity  response 
can  be  expressed  as. 


The  phase  of  the  velocity  response  then  becomes, 

»  =  tan-1<  fz0l/O  )  -  90°  .  (10) 

1  -  fr.l 

Equations  9  and  10  are  plotted  verses  the  frequency 
ratio  for  various  Q's  in  Figures  2.2  (a)  and  (b) . 
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Frequency  Ratio 

(a) 


(*) 

Figure  2.2  (a)  Normalized  velo 

f for  Q's  of  1 


o  o 


It  can  be  seen  from  these  two  plots  that  the  maximum 
velocity  amplitude  does  occur  at  resonance  (f-f0  or  f^«l) , 
no  matter  what  the  Q  for  the  system.  Further,  the  maximum 
rate -of -change  of  the  phase  always  occurs  at  resonance  which 
corresponds  to  <£»0°  for  the  velocity  response. 

B.  DTH1MTC  MODULI 

1.  Complex  Moduli 

The  technique  used  in  this  thesis  determines  both 
the  shear  and  Young's  complex  dynamic  moduli  by  measuring 
the  resonance  frequency  and  Q  of  the  fundamental  modes  of 
vibration  of  a  rod.  If  the  dimensions  of  the  rod  are  known, 
the  shear  modulus  can  be  calculated  from  the  speed  of 
torsional  waves  in  the  rod  and  the  Young's  modulus  can  be 
calculated  using  the  speed  of  either  the  longitudinal  or 
flexural  waves  in  the  rod. 

If  a  periodic  shearing  stress  is  applied  to  an 
elastomeric  material,  then  a  periodic  strain  having  the  same 
frequency  will  develop  in  the  material.  Figure  2.3  is  a 
sketch  of  an  oscillatory  stress  (<r)  and  the  resulting  strain 
(c),  between  which  there  exists  a  phase  angle  (0) .  The 
phase  angle  is  between  0  and  90°  and  is  dependent  upon  the 
amount  of  damping  present  in  the  material. 

The  complex  shear  modulus  (6*) ,  which  relates  the 
complex  stress  and  strain,  can  be  represented  by  two  vector 
components  [Ref.  8,9]  and  is  illustrated  in  Figure  2.2. 
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G*  a  G'  +  Gu 


(11) 


The  storage  or  elastic  shear  modulus ,  G',  is  defined 
as  the  ratio  between  the  stress  and  the  in-phase  component 
of  the  strain  [Ref.  8].  The  storage  modulus  is  proportional 
to  the  energy  stored  during  one  cycle  of  the  applied  stress. 


Figure  2.3  Time  profile  of  applied  periodic  shear 
stress  and  related  strain. 

The  viscous  or  loss  shear  modulus  (G”)  is  defined  as 
the  ratio  between  the  applied  stress  component  and  the 
quadrature  component  of  the  strain  [Ref.  9,10].  The  loss 
shear  modulus  is  a  measure  of  the  energy  dissipated  or  lost 
as  heat  during  one  cycle.  The  storage  and  loss  modulus  are 
related  by  the  tangent  of  the  phase  angle  (6)  as  shown  in 
Figure  2.4. 
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G' 


figure  2.4  Vector  diagram  of  the  complex  shear 
modulus . 

If  a  sinusoidal  tensile  stress  is  applied  to  the 
elastomeric  sample  instead  of  a  shear  stress,  then  the 
resulting  longitudinal  strain  will  be  out  of  phase  with  the 
applied  tensile  stress.  The  resulting  complex  Young's 
modulus  (E*)  is  related  to  its  storage  and  loss  components 
by, 

E'  *  E'  +  Eu  .  (12) 

The  complex  modulus  can  be  explained  using  a  forced 
damped  harmonic  oscillator  model  as  shown  in  Figure  2.5. 
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Figure  2.5  Damped  simple  harmonic  oscillator  with 
applied  force. 


The  equation  of  motion  for  this  system  is 

F(  t)  =  M9* +Ra  +K  X(  t)  ,  (13) 

where  X  is  the  displacement  from  equilibrium,  M  is  the  mass 
of  the  oscillator,  F.  is  the  mechanical  resistance,  K  is  the 
spring  constant,  and  F(t)  is  the  applied  force  [Ref.  11]. 

If  the  force  applied  and  resulting  motion  is  sinusoidal, 
then  the  expression,  after  appropriate  differentiation,  can 
be  written  in  complex  variable  form,  [Ref.  5] 

F  *  -w2MX*  +  juRX*  +  KX*  ,  (14) 

or 


F  =  -v>2MX*  +  K(  1  + 

A 


(15) 
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It  can  ba  saan  from  Equation  (15),  that  tha 
coefficient  of  tha  sacond  tarn  ia  tha  affective  stiffness 
for  tha  system.  This  affective  stiffness  is  analogous  to 
tha  complex  modulus  and  can  ba  expressed  as  follows. 


JT  -  *<i  ♦  -  *(i  ♦  jtanb)  .  (16) 

A 

Whan  the  driving  frequency  coincides  with  the 
mechanical  resonance  of  the  system  (the  frequency  at  which 
the  mechanical  reactance  vanishes)  then,  The 

storage  and  loss  components  of  the  effective  stiffness  are, 

K'  -  ;  K"  -  j-lL  .  (17) 

2.  Quality  Factor  and  Loss  Tangent 

The  Q,  as  seen  in  Figures  2.1  (a)  and  2.2  (a),  is  an 
indicator  of  the  sharpness  of  the  system  response  at 
resonance.  The  Loss  Tangent  is  the  ratio  of  the  energy 
dissipated  (EN  or  6")  to  the  average  rate  of  energy  stored 
(E#  or  G')  per  cycle  of  vibration,  or  the  reciprocal  of  Q, 


a» 

tanfi  «  — 
G' 


(18) 


The  quality  factor  for  a  system  can  be  defined  in  a 
variety  of  ways,  as  follows: 

1.  The  Q  is  the  ratio  of  the  reactance  of  the 
system  to  the  resistance. 
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(19) 


0  - 


2  %fjo 


2.  As  sssn  in  Figurs  2.1  (a) ,  Q  is  the  ratio  of 
the  displacement  amplitude  at  resonance  to  that  at  zero 
frequency. 


0  - 


x  </«f0) 

X(f~  0) 


(20) 


Equation  (20)  is  true  for  simple  harmonic  oscillators 
independent  of  the  value  for  Q,  including  Q's  less  than  1, 
assuming  the  stiffness  of  the  system  is  frequency 
independent.  For  the  case  of  a  viscoelastic  material,  a 
multiplicative  factor,  E(u**0)/E(ci>>‘<i>e) ,  is  required  to  make 
Equation  18  valid  for  dynamic  mechanical  systems.  Here,  E 
is  the  modulus  of  elasticity  corresponding  with  the  mode  of 
excitation. 


3.  The  Q  is  also  proportional  to  the  rate  of 
change  of  the  phase  at  resonance,  and  is  determined  by 
taking  the  derivative  of  0  (Eg.  11)  with  respect  to 
frequency: 


0-T#1".  ■  (21) 

The  derivation  [Ref.  11]  of  this  definition  is  found  in 
Appendix  A. 

4.  In  terms  of  the  relaxation  time  (r)  of  the 
system  the  Q  is. 
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(22) 


0  ■  */0T  , 

where  r-2»/R*  or  the  time  that  it  takes  for  the  free  decay 
amplitude  to  decrease  to  l/e  of  its  initial  value,  where 
In  e»l  [Ref.  10]. 

5.  Q  can  also  be  approximated  by, 

.  (23) 

where  f.  is  the  resonant  frequency  and  At  is  the  full 
bandwidth  over  which  the  displacement  has  dropped  to  V2  of 
its  value  at  resonance  (valid  for  Q>5) . 

6.  The  Q  can  also  be  found  from  a  pole  and  zero 
plot  of  the  mechanical  admittance  of  the  oscillator.  If  the 
characteristic  roots  (poles)  of  the  complex  admittance  of 
the  oscillator  are  expressed  as  y  «  a+jb  and  7*  -  a-jb, 
from. 


-2a 


(24) 


The  details  of  this  definition  of  Q  are  in  Appendix  B. 

7.  The  Q  shown  in  terms  of  the  critical  damping 


(25) 


where  f  is  the  fraction  of  critical  damping  for  the  system. 
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3.  ZaflUMM  of  Temperature  and  Frequency  on  tha 
Properties  of  Materials 

Tha  mechanical  propartias  of  viscoelastic  materials 
are  greatly  influenced  by  temperature  and  frequency.  In 
general ,  as  the  temperature  of  a  material  increases  the 
storage  modulus  decreases. 

The  transition  of  a  material  between  ideal  elastic 
solid  behavior  to  that  of  a  viscous  liquid  is  typically 
designated  by  the  glass  transition  temperature  (Tt) .  At 
temperatures  well  above  T,,  the  molecular  chains  in  a 
viscoelastic  material  are  able  to  slip  past  one  another  and 
begin  to  flow  due  to  their  increase  in  energy.  The 
increased  temperature  causes  the  material  to  behave  more 
like  a  viscous  liquid  provided  the  molecular  weight  of  the 
material  is  low.  For  materials  which  have  a  high  molecular 
weight  or  molecules  which  are  slightly  crosslinked,  the 
material  will  behave  in  a  rubber-like  fashion.  In  general, 
the  higher  the  molecular  weight  for  a  material  the  higher 
the  degree  of  intermolecular  bonding.  Therefore,  more 
energy  will  be  needed  to  break  those  bonds  resulting  in  a 
higher  melting  temperature  and/or  less  ability  for  the 
molecular  chains  to  flow  past  one  another.  [Ref.  8,9,12,13] 
The  rubber-like  behavior  occurs  because  the 
molecular  chains  within  the  material  are  extremely  twisted 
in  their  natural  (unstressed)  state,  and  become  untwisted 
whenever  a  stress  is  applied  to  the  material.  The  chains 
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revert  back  to  thair  natural  (twisted)  state  when  the  stress 
is  resowed.  This  process  will  continue  whenever  the  stress 
is  applied  or  reaoved.  Therefore,  the  material  is  capable 
of  hi9hly  elastic  deformation  when  in  this  rubber-like 
state.  [Ref.  8,9,12,13]  Below  T,  the  elastomer  is  said  to 
be  in  a  glassy  state.  This  glassy  state  of  the  material 
originates  in  the  atomic  structure  of  the  material.  \t  low 
temperatures,  the  thermal  and  free  volume  energy  of  the 
material  is  low  and  the  molecular  chains  can  no  longer  slip 
past  one  another.  In  other  words,  the  molecular  chain 
configurations  within  the  material  are  frozen  into  place 
causing  the  material  to  behave  more  like  an  ideal  solid. 
[Ref.  8,9,12,13] 

The  temperature  and  frequency  dependence  of  the 
material's  properties  may  be  significant  in  the  rubber  and 
glass  regions.  However,  this  dependence  is  most  pronounced 
when  the  viscoelastic  material  is  in  transition  between  the 
glass-like  (elastic  solid)  and  rubber-like  (viscous  liquid) 
states.  [Ref.  9,12,13] 

A  viscoelastic  material  typically  goes  through  four 
different  regions  or  transitions  of  dynamic  mechanical 
behavior  as  a  function  of  temperature  and  frequency.  These 
regions  are  defined  as  [Ref.  9,12]: 

1. )  A  glassy  region 

2. )  A  glass-rubber  transition  region 

3. )  A  rubbery  region 

4. )  A  flow  region 


19 


These  four  regions  are  illustrated  in  Figure  2.6  below. 


Figure  2.6  Temperature  and  frequency  dependence  of  the 
storage  modulus  and  loss  tangent  for 
viscoelastic  materials. 

Figure  2.6  shows  that  in  the  glassy  region,  the 
storage  modulus  has  a  maximum  value  and  is  relatively 
constant.  In  the  glass-rubber  transition  region,  the 
storage  modulus  decreases  rapidly  and  the  loss  tangent 
reaches  a  maximum  value.  In  the  rubbery  region  the  loss 
tangent  decreases  rapidly  and  then  levels  off.  The  storage 
modulus  stabilizes  and  is  relatively  independent  of 
temperature  and  frequency.  The  storage  modulus  decreases 
rapidly  in  the  flow  region;  however,  only  those  materials 
which  have  low  molecular  weight  and  virtually  no  cross 
linking  will  actually  begin  to  flow  [Ref.  9,12]. 


e.  moots  or  a  "trii-free"  bu 

A  long,  thin  rod  Bade  of  an  isotropic,  homogeneous  solid 
with  length  (L)  and  constant  cross-sectional  area  will 
propagate  vibrational  waves  in  three  independent  modes, 
provided  the  wavelength  of  vibration  (X)  is  much  greater 
than  the  rod  diameter.  The  boundary  conditions  existing  on 
the  rod,  coupled  with  the  physical  dimensions,  mass,  and  the 
elastic  moduli  will  determine  the  particular  resonance 
[Ref.  4]. 

In  this  application,  each  end  of  the  rod  is  "free", 
which  corresponds  to  that  of  2ero  moment  and  zero  stress  at 
the  boundaries  [Ref.  10].  The  physical  dimensions  of  the 
rod  are  normally  easy  to  measure.  By  measuring  the 
resonance  frequency,  the  elastic  properties  of  the  bar  can 
easily  be  determined  as  will  be  explained  next. 

1.  Longitudinal  Mode 

The  phase  speed  for  longitudinal  waves  is  [Ref.  8] 


where  p  is  the  mass  density  of  the  material  and  E  is  the 
Young's  modulus  of  elasticity. 

The  longitudinal  resonances  of  a  rod  having  free- 
free  boundary  conditions  are  easily  found  by  noting  that  the 
speed  (c) ,  wavelength  (X)  and  frequency  (f )  are  related  by, 
c~Xf  and  that  the  fundamental  wavelength  is  equal  to  twice 
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the  length  of  the  rod.  The  higher  order  nodes  for  a 
nondispersive  material  are  harmonically  related  and  can  be 
expressed  as  integer  multiples  of  the  fundamental.  Thus, 
the  longitudinal  resonances  of  a  "free-free"  rod  are: 


f*m~2L:  n“1'2*3 .  (27) 

where  n  is  the  mode  number  which  corresponds  to  the  number 
of  displacement  nodes  in  the  bar  standing  wave.  The  dynamic 
Young's  modulus  (EL)  can  be  obtained  from  the  previous  two 
equations: 

El  -  4pL*(-^jp)  (28) 

2 .  Torsional  Mode 

The  torsional  resonances  of  a  rod  having  free- free 
boundary  conditions  are  similar  to  that  of  the  longitudinal 
mode,  in  that  both  are  nondispersive  and  higher  order  modes 
will  also  be  integral  numbers  of  the  fundamental.  The 
torsional  resonance  frequencies  of  a  "free-free"  rod  are 


_  j*  OCf  ' 

a  “  2L  ' 
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The  phase  speed  for  the  torsional  waves  is 


cT  - 


where  G  is  the  shear  modulus  for  the  material. 


(29) 


(30) 


22 


The  dynamic  shear  modulus  can  be  obtained  from  the 
previous  two  equations: 

S-4pLa(  — )  (31) 

n 

From  equations  26  and  29  it  can  be  seen  that,  for  the 
longitudinal  and  the  torsional  modes,  both  the  dynamic 
Young's  and  shear  moduli  are  independent  of  the  cross- 
sectional  shape  of  the  rod. 

3.  Flexural  Mode 

The  phase  speed  for  the  flexural  wave  (cF)  is 
dispersive  and  varies  with  the  square  root  of  the  frequency. 
The  phase  speed  for  the  flexural  waves  is  [Ref.  10] 

cr  *  ^2*fk ct  ,  (32) 

where  the  radius  of  gyration  (k)  is  given  by  [Ref.  10] 

k*  =  <-| )fz3ds  ,  (33) 

where  S  is  the  cross-sectional  area  of  the  rod  and  z  is  the 
distance  of  an  element  above  the  neutral  axis  in  the 
direction  of  flexure.  For  a  circular  rod  of  diameter  d,  x 
is  d/4.  For  a  rectangular  bar  of  thickness  t,  x  is  t />/U 
[Ref.  10].  For  a  hollow  rod  the  radius  of  gyration  is 
V(oda-id*) /2,  where  od  is  the  outer  diameter  and  id  is  the 
inner  diameter  of  the  rod. 
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The  application  of  the  free-free  boundary  condition 
for  flexural  vibrations  leads  to  a  series  of  resonances 
given  by  [Ref.  8] 

f /  «  ;  a  -  3.0112,4.9994,7,9,11 .  (34) 

8L3 

where  L  is  the  length,  cF  is  the  flexural  wave  speed,  n  is 
the  node  number  and  x  is  the  radius  of  gyration.  This 
result  is  accurate  at  low  frequencies  where  the  effects  of 
rotary  inertia  and  shear  deformations  associated  with  the 
flexure  can  be  neglected.  [Ref.  4,14,15] 

The  dynamic  Young's  modulus  can  be  obtained  from 
Equations  32  and  34: 


Both  the  longitudinal  and  flexural  modes  are 
proportional  to  the  Young's  modulus.  Typically,  the 
resonant  frequencies  corresponding  to  the  longitudinal 
(compressional)  modes  are  an  order  of  magnitude  higher  than 
the  resonant  frequencies  of  the  flexural  mode.  Because  the 
redundancy  in  the  measurement  technique  is  easily  exploited, 
using  both  modes  provides  both  an  internal  check  of  the  data 
and/or  an  extended  frequency  range  over  which  the  Young's 
modulus  can  be  determined. 
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4.  Correction  for  Added  lad  Messes 

The  coils  and  the  epoxy  add  additional  mass  and 
stiffness  to  the  sample  which  is  under  investigation.  The 
additional  stiffness  is  usually  small  and  can  be  neglected 
for  those  samples  which  are  stiff  enough  to  support  their 
own  weight  and  shape.  The  added  mass  of  the  transducer 
coils  is  generally  only  a  few  percent  of  the  total  mass  of 
the  rod,  but  causes  slight  changes  in  the  resonance 
frequency.  This  change  can,  to  first  order,  be  accounted 
for  by  the  substitution  of  an  effective  length  (L^)  term 
into  the  equations  for  the  dynamic  moduli  [Ref.  4]. 

The  effective  length  term  can  be  found  from 
Rayleigh's  method  of  equating  the  maximum  kinetic  and 
potential  energies  of  a  vibrating  system  to  obtain  the 
resonance  frequencies  of  the  rod  [Ref.  16].  The  rod  is 
considered  to  be  a  simple  harmonic  oscillator  with  potential 
energy  of  l/2Kx0a  and  with  kinetic  energy  of  l/2mv0a .  Noting 
v^wXo  [Ref.  10]  and  equating  the  potential  and  kinetic 
energies  yields, 

=  -|i Dc|  -  <0=^  ,  (36) 

where  K  is  the  spring  constant,  M  is  the  mass  of  the  system 
and  <■>  is  the  angular  frequency.  The  change  in  resonance 
frequency  in  terms  of  the  change  in  the  kinetic  energy  is 
expressed  as  [Ref.  4] 
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(37) 


bf  _  6u>  _  1  SKB 
f  «  2  KB  ' 

Garrett  has  shown  [Ref.  4]  that  the  relative 
shift  in  the  kinetic  energy  is 

bKE  2m 
KB  AT  * 

Thus,  the  effective  length  for  the  longitudinal  node 
of  vibration  having  mass  added  to  each  end  (1/2  n)  is 

Lige  -  LU*-g)  •  (39) 

and  the  Young's  modulus  is  now  expressed  as, 

E  «  4p  (!,&,>*(  —  >  (40) 

J3 

The  same  procedure  is  used  to  find  the  effective 
length  in  the  torsional  case.  The  difference  is  that  the 
vibration  in  the  torsional  mode  is  generated  by  the 
restoring  force  within  the  material  and  the  moment  of 
inertia  due  to  the  mass  of  the  rod.  The  kinetic  energy 
related  to  the  rotation  of  the  rod  is 

KB  -  ±I(»eo)2  .  (41) 

It  has  also  been  shown  [Ref.  16]  that  the  relative  change  in 
kinetic  energy  for  the  torsional  mode  is 
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As  was  the  case  with  the  longitudinal  node,  the  length  of 
the  rod  is  inversely  proportional  to  the  torsional  resonant 
frequency.  The  relative  shift  in  the  frequency  and 
consequently  the  change  in  the  relative  length  of  the  rod  is 


6 L  m  -6f  ,  1  6KB  m  2m 
L  f  2  KB  M 


(43) 


and  the  effective  length  for  the  torsional  mode  becomes. 


L*gt  -  L(  1-4?)  . 


(44) 


The  corrected  shear  modulus  thus  becomes. 


G  •  4p(L,//)a(f//n)*  .  ( 

The  relative  shift  in  the  kinetic  energy  for  the 
transverse  (flexural)  mode  of  vibration  is  [Ref.  4] 


4.1m 

M 


(46) 


The  fundamental  resonant  frequency  for  the  flexural  mode  is 
proportional  to  the  square  root  of  the  energy  ratio  and 
inversely  proportional  to  the  square  of  the  length  of  the 
rod.  The  relative  change  in  length  is 


6L  m  -1  bf  m  -1  .  -1  6KBx  m  to 
L  2  f  2  2  KE  '  M 


(47) 


This  leads  to  an  effective  length  correction  term  for  the 
flexural  mode  which  is  identical  to  the  longitudinal  mode. 


2' 


(48) 


L«t  -  Ml  *  -g)  . 

and  tha  Young's  modulus  for  the  flexural  node  becomes, 

Bm  1024  .  (49) 

«*  d2  n 2 

The  effective  length  correction  accounts  only  for 
the  added  nass  presented  by  the  addition  of  the  transducer 
coils.  For  compliant  samples  the  additional  stiffness  added 
by  the  coils  would  need  to  be  accounted  for  as  well. 

However,  to  a  first  approximation,  the  argument  could  be 
made  that  the  added  stiffness  would  offset  the  correction 
for  the  added  mass,  making  neither  correction  necessary. 
Details  on  how  the  transducer  wire  coils  are  attached  to  the 
rod  in  order  to  excite  and  detect  the  vibrations  in  the 
sample  will  be  presented  in  Chapter  III. 

D.  MASTER  CURVES:  PRESENTATION  OF  DATA  AS  A  FUNCTION  OF 

SEDUCED  FREQUENCY 

1.  Introduction 

When  a  sinusoidal  stress  is  applied  to  a  material, 
the  frequency  and  temperature  dependence  of  both  the  storage 
moduli  and  loss  tangent  are  related.  This  relationship  is 
known  as  the  time-temperature  superposition  principle,  and 
is  illustrated  by  the  fact  that  the  same  effect  upon  the 
dynamic  material  properties  can  be  generated  either  by 
varying  the  temperature  and  holding  the  frequency  constant 
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or  by  varying  the  frequency  in  an  opposing  Banner  while 
aaintaining  constant  teaperature. 

It  is  difficult  to  distinguish  between  the  frequency 
and  the  teaperature  dependence  of  the  viscoelastic 
properties  of  the  aaterial.  As  was  seen  earlier  in  Figure 
2.6,  the  dependence  of  the  aaterial  properties  upon 
teaperature  and  frequency  changes  dramatically  in  the  glass 
transition  zone.  This  shows  that  the  relationship  between 
frequency  and  teaperature  is  even  nore  dramatic  in  this 
region.  A  method  of  reducing  the  variables  (frequency  and 
temperature)  was  developed  to  simplify  the  separation  of  the 
frequency  and  temperature  dependence  of  the  material 
properties  to  a  dependence  upon  a  single  variable.  This 
dependency  is  normally  expressed  in  terms  of  a  shift  factor, 
aT,  which  combines  the  temperature  and  frequency  into  one 
variable.  The  shift  factor  is  a  monotonically  decreasing 
function  of  temperature.  [Ref.  8] 

2 .  Williams-Landel-Ferry  Equation 

For  many  viscoelastic  materials,  the  general  form 
for  the  shift  factor  (ar)  is  expressed  by  the  Williams- 
Landel-Ferry  (WLF)  equation  [Ref.  8]  as, 
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logar 


cj{T-T9) 
C2  +  T  -  T0 


(50) 


where  c,°  and  cf  ara  empirically  datarainad  constants,  T  is 
tha  aaasurad  absoluta  temperature  in  kelvin  (K) ,  T0  is  a 
chosan  rafaranca  taaparatura  also  givan  in  K. 

Tha  constants  c,°  and  c,*  hava  baan  datarainad 
praviously  for  Polyaathyl  aathacrylata  (PHMA)  by  Parry 
kRaf .  8]  and  for  PR-1592  by  Capps  [Raf.  9].  Tha  aquation 
for  tha  corrasponding  shift  factor,  as  givan  by  Ferry,  for 
PUNA  is 


Log  ae 


-21.5CT-211) 
43.1  *  T  -  211 


and  as  givan  by  Capps  for  PR-1592  is 


(51) 


Log  a-  -  --1^(^83^151 
~  T  107  +T-  283.15 


(52) 


Tha  approxiaation  raprasantad  by  tha  WLF  aquation 
works  wall  for  aatarials  whose  glass  transition  taaparatura 
happens  to  fall  within  the  frequency  range  of  interest  such 
as,  PR1592 .  For  other  aatarials  such  as,  PHMA  and 
polycarbonate,  this  is  not  tha  case.  Without  T(,  the 
reference  taaparatura  becoaes  soaewhat  arbitrary  and  leads 
to  confusion  when  analyzing  data.  The  other  problea  with 
the  WLF  equation  can  be  seen  in  Equations  51  and  52  above. 
That  is  the  reference  temperatures  used  in  the  two  equations 
are  different.  Even  if  the  materials  had  been  the  same,  the 
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data  presented  by  each  could  act  bo  roadily  compared.  Tho 
data  arc  presented  bora  for  all  tliraa  aatorials  tested  in 
both  rev  fora  aad  in  terms  of  the  shift  factor  for  the 
broadest  possible  comparison. 

There  are  alternative  methods  for  approximating  a 
material's  dependence  upon  temperature  and  frequency.  One 
such  method  is  that  of  a  fractional  calculus  model  as 
described  by  Bagley  in  [Ref  17].  However ,  until  an 
approximation  is  demonstrated  to  be  suitable  for  the 
majority  of  viscoelastic  materials,  it  may  be  better  to 
simply  present  the  raw  data. 
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a.  zantooocTZov 

In  the  following  sections  tho  technique  for  measuring 
the  dynamic  properties  of  materials  will  be  outlined.  It 
should  be  noted  that  the  sample  rod  is  assumed  to  be  a 
lumped  acoustical  system,  and  all  the  equations  for  the 
resonance  and  Q  described  in  Chapter  II  pertain  to  the 
following  discussion. 

This  measurement  technique  is  both  precise  and 
relatively  inexpensive  to  perform,  and  since  it  is  resonant** 
based,  improves  the  signal-to-noise  ratio.  The  resonant 
frequency,  combined  with  knowledge  of  the  dimensions  and 
mass  of  a  homogeneous  and  isotropic  sample,  allow 
computation  of  the  material's  dynamic  moduli. 

The  resonant  frequencies  of  the  longitudinal,  torsional 
and  flexural  modes  of  vibration  are  measured  independently 
using  a  single  sample.  Thus,  both  the  dynamic  shear  and 
Young's  modulus  can  be  obtained.  This  is  quite  convenient, 
since  only  two  independent  moduli  are  needed  to  completely 
define  the  entire  set  of  the  mechanical  properties  of  an 
isotropic,  homogeneous  elastomeric  material  [Ref.  11]. 

The  basic  experimental  instrumentation  consists  of  the 
following  components: 
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1 •  A  test  supla  ud«  froa  uttrial  of  interest. 

2.  A  transducer  coil  nounted  oa  each  ends  of  the 
tost  Maple. 

3.  A  suspension  apparatus  used  to  support  the  test 
saaple  in  Mgnetio  fields  produced  by  a 

pair  of  peraanent  aagnets. 

4.  Electronic  instruaentation  used  to  excite  and 
detect  acoustic  vibration  of  the  saaple. 


1.  Saaple  Preparation 

Any  nonferrous  notarial  which  can  be  east,  drawn,  or 
aaohined  into  the  shape  of  a  long,  thin  rod  with  constant 
eross  section  can  be  investigated.  This  is  true  regardless 
of  what  that  cross  section  is. 

2.  Transducer  Coils 

The  transducer  coils  aounted  on  either  end  of  the 
saaple  are  identical.  Reciprocity  dictates  that  either  end 
of  the  rod  can  be  used  as  the  driver  or  the  receiver.  The 
coils  are  conprised  of  Mo.  32  gauge  insulated  copper  wire 
which  is  tightly  wound  into  a  4  turn  coil.  Each  4  turn  coil 
weighs  approxinately  0.3  g.  While  the  actual  coil 
paraneters  are  not  of  particular  inportance,  the  wire  used 
should  be  snail  to  reduce  mss  loading  effects  on  the  end  of 
the  saaple  but  have  sufficient  current  capacity  (typically 
50—200  Ma„J  (Ref.  4]. 

The  coils  are  shaped  to  fit  the  saaple  and  aligned 
in  a  siailar  orientation  on  each  end  of  the  saaple.  The 
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coils  art  then  attached  using  a  general  purpose  five -minute 
epoxy.  The  ends  of  the  insulated  coil  are  stripped  of  their 
insulation  and  solder  is  applied  to  ensure  good 
conductivity,  and  a  simple  continuity  check  is  performed. 

The  4  turn  coils  typically  have  a  resistance  of  1.1  n. 

It  should  be  noted  that  there  are  compromises 
associated  with  the  size  and  number  of  turns  used  in  the 
transducer  coils.  For  example,  if  the  turns  ratio  is 
increased  by  a  factor  of  N,  the  received  signal  will  be 
increased  by  a  factor  of  N2.  This  results  from  a  factor  of 
N  more  force  delivered  to  the  sample  and  a  factor  of  N 
higher  receive  sensitivity.  This  will  be  accomplished, 
however,  at  the  expense  of  added  mass  to  the  ends  of  the 
sample.  The  larger  the  mass  of  the  coil,  the  greater  the 
end  mass  correction.  [Ref.  5] 

3.  Physical  Properties  of  Test  Samples 

The  physical  properties  for  each  of  the  finished 
test  samples  are  listed  in  Table  3.1. 
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TABLR  3.1  PHYSICAL  PROPERTIES  OF  TEST  SAMPLE 


MATERIAL 

LENGTH 

DIAMETER 

RODMASS 

DENSITY 

TOT ALMAS S 

L  (m) 

d  (m) 

m  (g) 

p (kg/m3 ) 

*MT  (g) 

POLYCARBONATE 

0.3604 

0.0127 

54.72 

1192.9 

55.33 

PMMA  1 

0.3305 

0.0124 

49.41 

1237.9 

50.05 

PMMA  2 

0.3000 

0.0125 

43.4 

1186.4 

44.05 

PR-1592 

0.3602 

0.0123 

49.5 

1154.9 

50.60 

±. 0005m 

±. 0005m 

±.05g 

±.5 

±.05g 

*  MT  *  M  (rod)  +  m  (coils) 

C.  SUSPENSION  APPARATUS 


The  test  apparatus  used  in  this  research  is  a 
modification  of  one  used  previously  [Ref.  6,7],  and  is  shown 
in  Figure  3.1.  It  provides  for  sample  excitation,  position 
adjustment,  and  suspension.  The  excitation  system  is 
comprised  of  one  set  of  permanent  U-shaped  magnets  located 
at  each  end  of  the  apparatus,  between  which  the  driving  and 
receiving  coils  are  placed.  The  gap  between  the  magnet 
pairs  is  adjustable  from  0  to  approximately  3  cm.  The 
magnetic  field  strength  within  the  gap  is  approximately  2.4 
±  0.1  KOe  (0.24  ±  0.002  Tesla)  over  a  temperature  range  of  - 
17 ®C  to  85WC  [Ref.  5].  The  magnetic  remanence  for  these 
magnets  increases  at  lower  temperatures  [Ref.  18]. 
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Pigura  3,1  Resonance  based  dynamic  moduli  measurement 
apparatus . 


D«  ENVIRONMENTAL  CHAMBER 

An  environmental  chamber  is  used  to  house  the  samples 
thereby  permitting  temperature  control.  The  BHD-408  Bencl 
Top  Temperature  and  Humidity  Test  Chamber  (Associated 
Environmental  Systems)  used  [Appendix  C]  has  been  modifier 
to  accept  an  external  C02  canister  used  to  lower  the  minim  m 
temperature  to  approximately  -65°C.  The  upper  limit  of  this 
system  is  145°C.  The  temperature  control  stability  of  this 
system  is  specified  at  ±0,5°F  («0.3°C) .  The  internal 
dimensions  of  the  chamber  are  24  "x  2  4"x  24".  The 
environmental  chamber  is  capable  of  either  manual  or 
automated  operation  using  an  IEEE  488  interface. 
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■.  nucnvi  excitatiov  am d  detbctiom  07  desired  rksoxavt 

H00B8 

On*  of  tlio  most  important  features  of  the  free-free  bar 
method  of  measuring  the  dynamic  modulus  is  the  ability  to 
excite  the  longitudinal,  torsional  and  flexural  modes  of 
resonance  independently  using  the  same  inexpensive 
transducer.  The  sample  rod  of  the  material  is  positioned  in 
the  measurement  apparatus  such  that  the  transducers  mounted 
on  the  ends  of  the  rod  lie  within  the  gap  created  between 
the  permanent  magnet  poles.  This  places  the  transducers  in 
the  vicinity  of  the  maximum  magnetic  field  strength  during 
measurements.  The  orientation  of  the  transducer  coils 
within  the  magnetic  field  determines  which  resonance  mode  is 
excited.  The  induced  voltage  which  results  from  the 
movement  of  the  coils  within  the  magnetic  field  can  be 
expressed  mathematically  as  follows  [Ref.  4] 


where  Da  is  in  an  incremental  area  subtended  by  the 
transducer  wire  and  B,  the  magnetic  field.  The  induced  emf 
for  a  small  segment  of  wire,  length  1,  moving  with  velocity 
u,  in  B  is 


emf  =  B-IxU  . 


(54) 
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1.  Torsional  Mods 

The  torsional  node  of  vibration  is  excited  when  the 
transducer  coil  lies  within  the  nagnet  pair  as  shown  in 
Figure  3.2,  and  the  coil  driven  sinusoidally.  This  creates 
a  force  which  produces  a  torque  on  the  rod,  thus  generating 
torsional  waves  which  propagate  within  the  rod. 

A  similarly  placed  transducer  at  the  opposite  end  of 
the  sample  responds  to  the  torsional  vibration  and  generates 
an  electro-motive  force  (emf)  within  the  magnetic  field 
created  by  the  permanent  magnet.  The  flux  within  the  coil 
will  vary  as  a  result  of  the  change  in  the  angle  of  the  coil 
relative  to  the  magnetic  field. 


Figure  3.2  Orientation  of  the  transducer  coil  within 
gap  of  the  magnet  for  excitation  and 
detection  of  the  torsional  mode.  The  arrows 
indicate  the  direction  of  motion  created  by 
the  electromagnetic  force  on  the  coil. 
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2  .  Flexural  aeda 

The  flexural  node  is  excited  by  rotating  the  sample 
by  90°  and  translating  the  rod  up  or  down  approximately  a 
distance  equal  to  one  half  of  the  bar  diameter  such  that  one 
section  of  the  long  transducer  coil  sections  lies  in  a 
stronger  magnetic  field  as  shown  in  Figure  3.3.  The  forces 
generated  by  the  two  transducer  coil  sections  will  be  in 
opposite  directions;  however,  the  positions  of  the  coil 
sections  within  the  magnetic  field  gradient  are  such  that 
the  forces  will  not  be  equal.  The  two  opposing  forces  will 
cause  flexural  vibrations  to  propagate  along  the  axis  of  the 
rod.  The  flexural  motion  of  the  vibration  will  cause  the 
similarly-oriented  receiving  transducer  coil  at  the  opposite 
end  of  the  sample  rod  to  move  up  and  down  within  the  field. 
The  change  in  flux  through  the  transducer  coil,  induced  by 
this  movement,  generates  an  emf. 
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Figure  3.3  Orientation  of  the  transducer  coil  within 
gap  of  the  magnet  for  excitation  and 
detection  of  the  flexural  mode.  The  arrows 
indicate  the  direction  of  motion  created  by 
the  electromagnetic  force  on  the  coil. 

3 .  Longitudinal  Mode 

The  longitudinal  mode  of  vibration  is  excited  by 
increasing  the  distance  between  the  magnets  so  that  the 
magnetic  field  is  concentrated  on  the  small,  vertical, 
section  of  coil  that  spans  the  diameter  of  the  rod  at  the 
end  of  the  sample  as  shown  in  Figure  3.4.  The  currents 
within  this  small  section  of  the  transducer  coil  generate  a 
longitudinal  force  on  the  end  which  will  propagate  as 
longitudinal  vibration  along  the  length  of  the  rod. 
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Figure  3.4  Orientation  of  the  transducer  coil  within 
gap  of  the  magnets  for  excitation  and 
detection  of  the  longitudinal  mode.  The 
arrows  indicate  the  direction  of  motion 
created  by  the  electromagnetic  force  on  the 
coil. 


The  longitudinal  motion  of  the  vibration  will  cause 
the  receiving  transducer  coil  at  the  opposite  end  of  the 
sample  rod  to  move  in  and  out  of  the  magnetic  field.  The 
change  in  flux  through  the  transducer  coil,  induced  by  this 
movement,  generates  an  emf. 


1.  Introduction 


The  electron ic  instrumentation  used  to  detect  the 
induced  emf  can  range  in  sophistication  from  an  ordinary 
oscillator  and  voltmeter  to  a  computer-controlled  logging 
and  analysis  system.  There  are  advantages  and  disadvantages 
associated  with  the  chosen  instrumentation.  The  basic 
components  needed  are: 

1.  A  suspension  apparatus 

2 .  A  driver 

3 .  A  receiver 

4 .  Magnets 

2 .  Electronic  Instrumentation 

a.  Fundamental  Instrumentation 

The  fundamental  instrumentation  scheme  needed  for 
taking  measurements  consists  of  an  audio  oscillator  or 
function  generator  to  drive  the  sample  rod  and  a  voltmeter, 
frequency  counter,  and/or  an  oscilloscope  to  measure  the 
received  signal.  The  advantages  of  using  these  instruments 
are  the  relatively  low  cost,  availability,  and  simplicity  of 
use.  The  disadvantages  for  using  only  these  instruments  are 
the  limited  capabilities,  and  the  time  involved  in  taking 
point-by-point  measurements .  The  schematic  representations 
for  the  various  configurations  of  these  components  are  shown 
in  Figure  3.5. 
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Figure  3.5  Block  diagram  of  a  simple  instrumentation 
setup  for  room  temperature  measurements . 

Here,  this  simple  instrumentation  scheme  is  further  enhanced 

by  adding  additional  electronic  equipment.  An  Ithaco  1201 

pre-amplifier  (with  filter)  is  used  to  amplify  the  signal 

from  the  receiver  coil.  A  high  pass-filter  with  a  10  Hz 

rolloff  connected  to  the  preamplifier  effectively  removes 

the  fluctuations  on  the  oscilloscope  caused  by  the  soft 

suspension  system  or  by  external  environmental  vibration. 

On  the  driving  side,  an  HP492A  power  amplifier  improves  the 

coupling  of  the  driver's  output  impedance  to  the  low- 

impedance  transducer. 

t>_  Tn^t-n;  ,  rotation  Used 

More  sophisticated  equipment  is  available  which 

will  enhance  the  capabilities  and  decrease  the  time  needed 

to  make  these  measurements .  An  HP3562A  spectrum  analyzer  is 

> 
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used  for  this  purpose.  The  spectrum  analyzer  is  used  as 
both  the  signal  source  and  the  receiver,  so  that  the 
frequency  response  of  the  saterial  sample  is  readily 
displayed  on  the  CRT.  The  HP3562A  also  has  many  other 
useful  features,  such  as  pole-zero  curve  fit  capabilities 
and  the  ability  to  measure  the  frequency  response  (magnitude 
and  the  phase  change)  with  respect  to  a  reference  input. 

With  the  HP3562A,  one  can  manipulate  the  data  with  simple 
math  calculations  such  as  multiplying  by  i/jw  to  produce  the 
displacement  response  of  the  system.  The  spectrum  analyzer 
also  makes  it  very  simple  to  identify  the  various  resonant 
modes.  A  schematic  representation  of  the  instrumentation 
actually  used  for  the  room  temperature  evaluation  is  shown 
in  Figure  3.6. 


HMMA 


Figure  3.6  Block  diagram  of  the  FFT  based 

instrumentation  used  for  room  temperature 
measurements . 
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o.  Lock-in  Amplifier 

The  lock-in  amplifier  is  a  phass  sansitiva 
datactor  and  low  pass  filter  which  is  capabla  of  datacting 
small  signals  in  tha  prasanca  of  noise.  Tha  lock-in 
asplifiar  is  an  assantial  part  of  tha  Phasa  Locked  Loop 
(PLL)  and  is  described  in  detail  in  section  G  of  this 
chapter.  The  lock-in  asplifiar  is  also  an  excellent 
instrument  for  taking  room  temperature  measurements. 

3.  Determination  of  Resonance  Frequency 

As  shown  in  Figure  2.5(a)  the  maximum  velocity 
response  can  always  be  used  to  find  the  resonance  frequency. 
Another  useful  method  for  determining  the  resonance 
frequency  is  to  determine  where  the  phase  crosses  0°  (for 
velocity  measurements)  using  the  HP3562A  spectrum  analyzer 
as  described  in  Figures  2.4(b)  and  2.5(b).  Probably  the 
most  accurate  method  for  measuring  the  resonance  frequency 
is  to  use  the  pole-zero  curve  fit  from  the  HP3562A  as 
described  in  Appendix  B. 

4.  Measurement  of  the  Quality  Factor  (Q) 

The  main  method  used  for  determining  the  quality 
factor  (Q)  here  is  that  of  using  the  pole  zero  capabilities 
of  the  HP3562A  dynamic  signal  analyzer.  The  HP3562A  dynamic 
signal  analyzer  has  an  internal  curve  fitting  algorithm 
which  provides  a  synthesized  curve  fit  of  poles  and  zeroes. 
The  Q  and  resonance  frequency  for  the  fundamental  resonance 
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and  its  ovsrtoass  srs  sssily  obtained  from  thsss  poles  and 
aeroes.  This  technique  has  bean  described  by  Brown  [Bef. 

19]  and  is  presented  in  further  detail  in  Appendix  B.  As 
further  verification  of  the  Q  for  the  system  was  also  found 
by  Measuring  the  resonance  frequency  and  the  frequencies 
where  the  amplitude  is  vT  (3dB)  down  from  the  resonance 
value.  This  is  a  valid  approximation  whan  the  Q  is  greater 
than  S. 

a.  rbsobawcb  lucms 
1.  Introduction 

When  the  temperature  of  the  sample  material  is 
varied,  the  resonance  frequency  and  the  Q  of  the  material 
also  vary.  The  resonant  free-free  bar  technique  described 
earlier  in  Section  D,  can  be  used  to  measure  the  storage 
modulus  and  the  loss  modulus  components  as  function  of 
temperature  and  frequency  [Ref.  4,7].  Since  this  is  a 
resonant  based  method,  the  sample  must  be  kept  on  resonance 
while  the  temperature  is  changing.  This  can  be  accomplished 
manually  by  readjusting  the  drive  frequency  to  coincide  with 
resonance  every  time  the  temperature  changes.  Here,  a  two 
channel  phase-locked  loop  (PLL)  is  used  to  continuously 
track  the  resonant  frequency  as  a  function  of  temperature. 

As  will  be  shown  later  in  this  chapter,  the  quadrature 
signal  from  the  PLL  is  used  to  control  the  frequency  of  a 
Voltage  Controlled  oscillator  (VCO)  and  to  maintain  the 
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sample  at  rasonanca.  Tha  in-phase  amplitude  of  the  PLL  is 
shown  in  Chapter  III,  Section  G,  Subsection  2  to  be 
proportional  to  Q  and  is  used  to  determine  the  change  in  the 
damping  properties  of  the  material  as  the  temperature  and 
frequency  are  varied. 

2.  Phase  Looked  Loop 

The  phase  locked  loop  (PLL)  is  used  to  ensure  that 
the  sample  will  remain  at  resonance  while  the  temperature  is 
varied. 

a.  Introduction  to  PLL  circuits 

A  PLL  consists  of  three  main  components,  a  phase 
sensitive  detector  (a  mixer  or  multiplier) ,  a  low  pass 
filter,  and  a  voltage  controlled  oscillator  as  shown  in 
Figure  3.7. 


Figure  3.7  Block  diagram  of  a  basic  PLL  circuit. 

A  mixer  is  a  circuit  element  whose  output  is  the 
product  of  its  two  inputs.  Typically,  one  input  is  the 
signal  of  interest  and  the  other  is  a  reference  signal. 
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Whan  the  two  inputs  ara  sinusoidal,  the  output  frost  the 
sixer  is  two  sinusoids  having  frequencies  equal  to  the  sun 
and  difference  frequencies  of  the  two  inputs.  [Ref.  20] 

o (t)  ■  -|cos[  t+ -42) 3— |cos [(«!+«*)  fc+^+^a)]  #  (55) 

where  e(t)-e,+e2,  e,»sin(«,t+^,)  and  e^-sin^t+^J .  When  the 
two  input  signals  are  identical  in  frequency,  the  difference 
frequency  will  be  zero;  however,  there  could  still  be  a 
phase  difference  between  the  two  signals. 

The  output  signal  fron  the  mixer  is  delivered  to 
a  low  pass  filter  (LPF) .  The  low  pass  filter  removes  the 
sum  component  from  the  mixer  output  signal  leaving  only  the 
difference  component.  The  output  from  the  LPF  is  a  dc 
voltage  whose  amplitude  is  proportional  to  the  difference  in 
phase  between  the  two  signals.  A  secondary  effect  of  the 
low  pass  filter  is  to  allow  low  frequency  components  of  the 
mixer  output  signal  to  pass  through  while  the  higher 
frequency  components  are  attenuated.  This  has  the  effect  of 
reducing  the  noise  associated  with  the  mixer  output  signal. 
The  amount  of  attenuation  depends  upon  the  RC  time  constant 
and  the  equivalent  noise  bandwidth  (ENBW)  of  the  low  pass 
filter. 

Often,  the  reference  signal  used  with  a  lock-in 
amplifier  is  a  square  wave.  In  this  case,  the  output  from 
the  mixer  will  be  composed  of  a  large  number  of  sum  and 
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difference  frequencies.  In  this  case,  the  mixer  will  again 
produce  a  phase  sensitive  dc  output  when  the  signal 
frequency  is  equal  to  the  fundamental  or  odd-harmonic  of  the 
reference. 

The  voltage  controlled  oscillator  (VCO)  has  a 
free-running  frequency  which  may  be  modified  by  the 
application  of  an  external  control  voltage  [Ref.  20].  The 
change  in  frequency  is  proportional  to  this  control  voltage. 
In  other  words,  the  VCO  adjusts  its  output  frequency 
according  to  this  feedback  voltage. 

The  output  signal  from  the  VCO  is  used  as  the 
input  reference  signal  to  the  mixer.  The  circuit  is 
"locked"  if  the  frequency  of  the  reference  signal  (VCO) 
equals  the  signal  of  interest  and  a  fixed  phase  difference 
(0  degrees  in  our  case)  between  the  two  signals  is 
maintained. 

b.  PLL  utilizing  a  two  channel  lock-in  analyzer 
Figure  3.8  shows  a  schematic  diagram  of  a  basic 
PLL  circuit  utilizing  a  two  channel  lock-in  analyzer  for  the 
phase-sensitive  detection.  A  two  channel  lock-in  amplifier 
contains  two  individual  phase-sensitive  detector  (PSD) 
circuits.  Both  PSD  circuits  use  the  signal  driving  the 
excitation  coil  as  the  reference  and  the  received  signal  as 
the  input.  One  of  the  channels  uses  a  reference  circuit 
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phase  shifter  to  change  the  phase  of  the  signal  by  90°  in 
order  to  produce  a  quadrature  output. 


Figure  3.S  Block  diagram  of  the  basic  instrumentation 
used  to  track  the  resonance  of  the  rod  as  a 
function  of  temperature. 

Accordingly,  two  dc  output  signals  are  provided,  one  of 
which  is  in  phase  with  the  signal  of  interest  and  the  other 
in  quadrature . 

By  utilizing  a  feedback  loop  between  the 
quadrature  output  and  the  reference  circuit  phase  shifter, 
an  operational  integrator  can  be  formed.  The  output  from 
this  integrator  can  then  be  used  to  correct  any  phase  shift 
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b«tv««a  the  rcftrtnot  oirouit  ia  the  VCO  mad  the  iaput 
aigaal.  Whenever  a  phaac  shift  is  dstsetsd  between  tbs 
quadrature  iaput  sigaal  aad  ths  rsfsrsaos  sigaal,  tha 
integrator  feedback  loop  will  produoa  a  dc  arror  sigaal 
proport ioaal  to  tha  phasa  diffaraaoa  batwaaa  tha  two 
sigaals.  Tha  affaot  of  tha  faadbaok  loop  is  to  continually 
foroa  tha  output  from  tha  quadratura  mixer  to  ba  zaro. 

o.  Application  of  PLL  to  traok  rasoaaaoa  of  sample 
ha  HF3314A  Function  Ganarator  is  used  to  drive 
tha  sample  rod  at  a  predetermined  resonance  frequency  which 
corresponds  to  a  certain  tempera ture.  Tha  HP3314A  function 
ganarator  also  has  tha  capability  of  providing  tha  reference 
circuit  discussed  earlier  via  an  internal  vco.  Tha  VCO  has 
a  control  range  of  4-10%  to  -100%  of  tha  free-running 
frequency.  However,  tha  linear  operating  range  for  the  VCO 
function  appears  to  ba  only  4-10%  to  -80%  of  tha  preset 
output  frequency  (used  to  driva  tha  sample  rod) .  since  the 
resonance  frequency  of  the  same  material  increases  as  the 
temperature  decreases,  the  system  is  usually  locked  at  a  low 
temperature. 

The  velocity  of  the  motion  of  the  receiving 
transducer  coil  within  the  magnetic  field  at  the  opposite 
end  of  the  sample  rod  induces  an  emf  voltage  within  the 
coil.  Due  to  the  type  of  rocking  motion  generated,  the 
output  signal  will  have  a  frequency  and  a  phase  related  to 
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temperature  and  the  stiffness  of  the  sample  rod.  This 
voltage  signal  is  filtered  and  used  as  the  input  to  the  two 
channel  lock-in  amplifier.  A  square  wave  trigger  output 
from  the  frequency  generator  is  used  as  the  reference  signal 
in  the  lock-in  amplifier.  The  quadrature  output  signal  of 
the  lock-in  is  use*'  as  the  feedback  loop  to  the  vco  via  an 
external  low  pass  filter.  When  the  system  is  locked,  the 
temperature  of  the  sample  rod  can  be  varied  and  the 
resulting  change  in  its  resonance  frequency  will  produce  a 
phase  error  at  the  quadrature  mixer.  The  feedback  loop  will 
force  the  reference  frequency  to  change  to  that  of  the  input 
signal.  The  quadrature  signal  will  be  forced  back  to  zero 
and  the  in-phase  signal  will  be  a  maximum. 

The  new  reference  frequency,  maximized  in-phase 
voltage,  and  the  temperature  of  the  sample  rod  can  be 
recorded.  Knowing  these  three  parameters  and  the  physical 
dimensions  of  the  sample  rod,  the  dynamic  modulus  and  the 
loss  properties  of  the  sample  rod  material  can  be  determined 
as  a  function  of  both  temperature  and  frequency. 

3 .  Electronic  Instrumentation 

The  components  which  were  used  for  the  resonant 
tracking  measurement  include: 

1.  Suspension  apparatus  (w/  sample  and  attached 
transducers) 

2.  Phase-locked  loop  consisting  of 

a.  voltage  controlled  oscillator 
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b.  integrating  circuit 

c.  two  channel  look-in  analyser 

Controllable  tenperature  chamber 

4 .  Tenperature  gage 

5 .  Frequency  counter 

6 .  Voltmeter 

Figure  3.9  is  a  block  diagram  of  the  system  of  the 
instrumentation  used  for  the  resonance  tracking  measurement. 
The  temperature  control  chamber  stabilises  the  test  sample 
at  the  selected  temperatures  so  the  measurements  can  be 
taken  while  the  sample  rod  is  in  thermal  equilibrium.  A 
thermistor  and  a  voltmeter  are  used  to  record  the  actual 
tenperature  on  the  sample  rod.  The  frequency  counter  is 
used  to  record  the  resonance  frequency  of  the  sample  rod  at 
various  temperatures.  The  voltmeter  is  used  to  record  the 
in-phase  output  voltage  which  is  proportional  to  the  Q.  A 
power  amplifier  and  a  pre-amp  were  added  to  improve 
performance.  To  facilitate  the  data  collection  process,  the 
entire  measurement  system  was  automated  using  an  HP9836 
computer.  Each  voltmeter  and  the  frequency  counter  are 
equipped  with  an  IEEE488  bus  interface. 
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Figure  3.9  Block  diagram  of  the  instrumentation  used  to 
track  the  dynamic  moduli  as  a  function  of 
temperature  and  frequency. 

4.  Measurement  of  the  Dynamic  Moduli  using  a  PLL 

A  sample  is  first  driven  at  a  selected  resonant  mode 
by  adjusting  the  VCO  manually  with  the  integrator  shorted 
and  the  error  signal  feedback  path  open  ( quadrature  output 
from  the  lock-in)  so  that  no  error  signal  is  present  at  the 
input  to  the  VCO.  The  resonance  frequency  of  the  sample  rod 
is  determined  by  either  the  maximum  velocity  amplitude  or 
the  zero  crossing  of  the  phase  between  the  driving  frequency 
and  received  signal.  When  the  drive  frequency  from  the  VCO 
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and  tha  rasonant  fraquency  of  the  sample  rod  are 
synchronized,  the  quadrature  error  siqnal  will  be  a  phase 
sensitive  dc  voltage  as  described  earlier.  The  output  from 
the  quadrature  signal  channel  of  the  lock-in  amplifier  is 
set  to  zero  by  adjusting  an  internal  phase  shifter.  The 
integrator  is  then  opened  and  the  feedback  loop  reconnected. 

Once  the  control  loop  becomes  locked  and  stable,  the 
temperature  of  the  sample  rod  can  be  altered.  As  the 
temperature  changes,  the  natural  resonance  frequency  of  the 
rod  shifts.  The  lock-in  detects  the  frequency  shift  and 
alters  the  quadrature  output  voltage.  This  voltage  is  then 
used  as  a  feedback  signal  to  "control"  the  VCO  output 
frequency  so  that  it  continually  matches  the  resonant ' 
frequency  of  the  sample  rod.  The  resonance  frequency  can 
either  be  recorded  directly  from  the  VCO  or  from  an  external 
frequency  counter.  (A  thermistor  connected  to  a  digital 
multimeter  is  used  to  record  the  sample  rod  temperature.) 

An  HP9000  series  computer  is  used  to  record  the  resonant 
frequency,  the  sample  rod  temperature  and  the  in-phase 
voltage.  The  control  program  was  originally  written  by 
Beaton  [Ref.  €],  modified  by  Tan  [Ref.  7],  and  modified 
again  by  this  author  [Appendix  0] . 

5.  Measurement  of  the  Q  (and  Loss  Tangent)  using  a  PLL 

The  in-phase  voltage  from  the  lock-in  amplifier  is 
proportional  to  the  Q  of  the  system  when  the  system  is 
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looked  onto  rtioBtnot.  The  q  ou,  therefore,  be  obtained  at 
a  function  of  temperature  by  raoording  tba  in-phase  voltaga 
from  tba  look-in  amplifier.  Tba  ralationsbip  between  tba 
in-pbasa  voltaga  and  tba  loss  tangent  can  be  derived  from 
tba  modal  of  a  damped  simple  barmonie  oscillator  driven  at 
resonance.  (Shown  previously  by  Brown  and  Garrett  in 
Reference  5.) 

Resonance  of  a  mechanical  system  is  defined  as  the 
frequency  at  which  the  input  mechanical  reactance  goes  to 
sero  (Ref.  10].  Therefore,  the  impedance  of  the  mechanical 
system  at  resonance  is  the  resistance  of  the  system.  Voting 
the  following  definition  of  the  quality  factor,  [Ref.  12] 


0  «  ,  (56) 

where  u0  is  the  angular  resonance  frequency,  a  is  the  mass, 
R  is  the  resistance  of  the  system  (at  resonance  is  equal  to 
the  applied  force  over  velocity,  R=P/u) ,  the  Q  then  equals. 


q  = 


u„  a  u 


(57) 


Rearranging  Equation  55  yields  the  following  expression. 


Q  _  m 
u«tf  F 


(58) 


If  the  reasonable  assumption  is  made  that  the 
driving  force  can  be  held  constant  (independent  of  both 
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temperature  and  fraquancy)  and  that  the  mass  of  the  sample 
doas  not  change ,  than  Equation  58  is  also  constant. 

The  velocity  of  the  system  is  proportional  to  the 
emf  generated  in  the  coils  of  wire.  The  velocity  in  the 
above  equation  for  Q  can,  therefore,  be  replaced  by  the  in- 
phase  voltage  (Vk)  from  the  lock-in  amplifier  for  a 
particular  temperature  and  frequency.  This  results  in  a 
proportionality  constant  (A)  which  is  derived  as  follows, 


in*"o 


(59) 


Therefore,  the  loss  tangent  as  a  function  of  temperature  and 
frequency  is. 


tanfi  »  4  * 


(60) 


Q  A  V^«e 

A  disadvantage  of  this  system,  as  can  be  seen  in 
Equation  59,  is  that  an  independent  measurement  of  the  Q  for 
the  system  must  be  taken  in  ord«r  to  determine  the 
proportionality  constant,  A.  The  accuracy  of  this 
calibration  is  directly  dependent  upon  maintaining  the  exact 
temperature  and  frequency  between  the  two  measurements. 

Even  so,  this  measurement  is  easily  accomplished  by  using 
the  room  temperature  method  discussed  in  Section  F  of  this 
chapter. 
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XP.  MRTBRXRL  IMVUTXOATXOVS 

A.  XMTRODOCTXOM 

Data  vaa  collected  for  Polyurethane  PR-1592,  Polymethyl 
methacrylate  (PMMA) ,  and  Polycarbonate.  Of  these  materials, 
PR-1592  is  the  only  one  which  has  a  glass  transition 
temperature  within  the  normal  operating  range  of  this 
system.  Both  the  PMMA  and  the  polycarbonate  are  relatively 
stable  materials  over  the  temperature  range  tested.  For 
each  material  tested,  the  following  sections  provide 
discussion  and  selected  data  for  both  samples  at  room 
temperature  and  resonant  tracking. 

B.  POLYURBTHAMB  PR-1592 

The  PR-1592  was  included  in  this  report  in  order  to 
demonstrate  the  effect  of  the  glass  transition  temperature 
upon  the  complex  elast  'oduli  and  loss  tangent  of  a 
material.  PR-1592  is  a  commonly  used  sonar  transducer 
encapsulant.  An  accurate  characterization  of  its  properties 
is  a  matter  of  importance.  This  material  was  discussed  in 
further  detail  by  Tan  [Ref.  7]. 

1.  Material  Properties  at  Room  Temperature 

Table  4.1  presents  the  results  of  the  torsional  mode 
(shear  modulus)  for  the  PR-1592  test  sample  at  room 
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temperature.  Figurs  4.1  clearly  shows  ths  first  3 
resonances  sad  thsir  corresponding  pole-sero  data  for  saeh. 

TABLE  4.1  FR-1592  SHEAR  MODULUS  RESULTS  (Ts2S.4*C). 


f 

f/n 

Q 

Q 

Q(Leff) 

(Hi) 

<KS) 

<3dB) 

(CURVE) 

(GPa) 

256.7 

256.7 

2.04 

2.95 

0.043 

532.2 

266.1 

2.04 

4.06 

0.046 

811.6 

270.5 

3.06 

2.07 

0.048 

2.  Essoaaat  Trade  lag  Evaluation  of  PR-1592 

Ths  in-phase  voltags  from  ths  rssoaanes  tracking  of 
ths  PMMA  saapls  was  calibratsd  by  dirsctly  obtaining  ths  Q 
at  ssvsral  diffsrsnt  tsmpsraturss  using  ths  HP3562A  Spsotrun 
Analyssr.  Ths  proportionality  constant  (A)  obtained  fro* 
this  calibration  was  than  used  to  derive  ths  Loss  tangent  as 
a  function  of  frequency  and  teaperature.  The  results  of 
these  measurements  for  the  first  two  torsional  resonances 
are  shown  in  Tables  4.2  and  4.3.  Separate  plots  of  this 
data  as  a  master  curve  for  both  of  the  first  two  torsional 
modes  are  presented  in  Figures  4.6  and  4.7  directly 
following  the  master  curves  plotted  from  the  actual 
resonance  tracking  data  as  further  verification  of  this 
measurement  technique. 
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Figure  4.1  PR-1592  frequency  response  curves  and  pole 

zero  information  for  the  torsional  modes  at 
T-25.4°C. 
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nuns  4.2  IN-PHASE  VOLTAGE/QUALITY  FACTOR  CALIBRATION 
FOR  THE  FIRST  TORSIONAL  NODE  OF  PR-1592. 


1m  alia 

FREQUENCY 

QUALITY 

IN-PHASE 

"A" 

T 

f 

FACTOR 

VOLTAGE 

Q/(fV) 

CC) 

(Hz) 

Q 

V  (DCV) 

(msec/V) 

-16.75 

787.5 

4.76 

N/A 

N/A 

-10.16 

741.5 

3.05 

-0.2006 

-19.31 

-5.67 

579.3 

3.40 

-0.2034 

-28.85 

4.50 

427.0 

2.87 

-0.2288 

-29.38 

14.23 

329.9 

2.40 

-0.2859 

-25.44 

25.40 

256.7 

2.95 

-0.3879 

-29.63 

37.02 

222.1 

3.54 

-0.5149 

-30.95 

46.42 

199.0 

4.17 

-0.6673 

-31.40 

55.85 

187.0 

4.82 

-0.8403 

-30.67 

67.01 

177.2 

6.08 

-1.0734 

-31.97 

74.60 

175.8 

6.23 

-1.2534 

-28.27 

84.60 

171.2 

7.40 

N/A 

N/A 

AVERAGE  (A) : 

-28.58 

Figures  4.2  and  4.3  are  two  of  the  actual  data 
curves  taken  at  the  opposite  extremes  of  the  temperature 
ranges.  These  two  figures  are  shown  to  demonstrate  the 
effect  of  the  glass  transition  temperature  upon  the  dynamic 
material  properties  of  viscoelastic  materials.  All  of  the 
curves  with  the  data  used  to  determine  the  proportionality 
constant  for  PR-1592  can  be  found  in  Appendix  E. 
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TAILS  4.3  IN-PHASE  VOLTAGE/QUALITY  FACTOR  CALIBRATION 
FOR  THE  SECOND  TORSIONAL  MODE  OF  PR-1592. 


FREQUENCY 

f 

(H2) 

QUALITY 

FACTOR 

Q 

IN-PHASE 

VOLTAGE 

V  (DCV) 

"A" 

Q/(FV) 

(BStC/V) 

-16.75 

1636.9 

5.25 

N/A 

N/A  ! 

-10.16 

1322.9 

2.93 

-0 

.1935 

-11.44 

-5.67 

1216.4 

3.49 

-0 

.1874 

-15.31 

4.50 

898.7 

2.99 

-0 

.1962 

-16.96 

14.23 

681.4 

3.40 

-0 

.2475 

-20.16 

25.40 

532.2 

4.06 

-0 

.3521 

-21.66 

37.02 

495.5 

4.61 

-0 

.4846 

-19.20 

46.42 

423.6 

4.84 

-0 

.6205 

-18.41 

55.85 

394.9 

5.88 

-0 

.8041 

-18.52 

67.01 

374.8 

6.47 

-1 

.0487 

-16.46 

74.60 

358.9 

10.15 

-1 

.2205 

-23.17 

84.60 

350.8 

9.79 

N/A 

N/A 

AVERAGE 

(A): 

-18.13 
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Figure  4.2  PR-1592  direct  measurement  of  Q,  T=-16.75°C. 
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rigor*  4.3  PR-1592  direct  measurement  of  Q,  T=84.6°c. 
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Figure 


4.4  Heater  curve  for  the  first  torsional  mode  of 
PR- 1592. 
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UJ 
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4.5  Hester  curve  for  the  second  torsional  mode 
of  PR-1592. 
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Figure  4.6  Master  curve  for  the  first  torsional  mode  of 
PR- 1592  using  data  obtained  from  direct 
measurement  of  Q  at  various  tempera ture. 


♦  MODULUS 
O  LOSS  TANGENT 
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REDUCED  FREQUENCY  (Hz) 


Figure  4.7  Master  curve  for  the  second  torsional  mode 
of  PR- 1592  data  obtained  from  direct 
measurement  of  Q  at  various  temperature. 


The  glass  transition  of  viscoelastic  materials 
explained  earlier  is  clearly  seen  in  Figures  4.4  and  4.5  for 
the  first  and  second  torsional  modes  of  the  PR-1592  sample. 
Figures  4.6  and  4.7,  which  are  plotted  from  the  data  taken 
in  the  independent  measurements  for  Q,  provide  further 
verification  of  this  phenomena.  During  the  transition 
between  glass-like  to  rubber-like  states  the  Q  reaches  its 
minimum  value  (thus,  the  loss  tangent  has  its  maximum) .  The 
plots  also  show  that  the  resonance  frequency  has  higher 
values  at  the  lower  temperatures  and  a  lower  value  at  the 
higher  temperatures.  In  these  four  figures  the  shear 
modulus  and  loss  tangent  are  plotted  verses  reduced 
frequency.  The  reference  temperature  used  to  derive  the 
reduced  frequency  is  283.15  K. 

While  Figures  4.6  and  4.7  essentially  verify  the 
measurements  taken  using  the  PLL,  they  are  also  an  excellent 
example  of  the  of  the  unique  differences  in  the  measurement 
techniques  used  to  obtain  the  data.  The  data  taken  using 
the  resonance  tracking  technique  involving  the  PLL  is 
continuous;  providing  a  complete  analysis  of  the  material 
for  the  single  resonance  over  a  wide  range  of  temperatures. 
The  spectrum  analyzer  data  synthesis  provides  complete 
information  about  all  of  the  harmonic  resonances  for  a 
particular  vibrational  mode  at  a  single  temperature. 
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hovtvar,  it  tbn  toooua  vary  difficult  to  distinguish 
between  tbs  first  and  second  nodes. 


Tables  4.4  and  4.5  are  a  sunnary  of  data  collected 
for  the  first  and  second  torsional  nodes  of  the  PR-1592 
sanple.  The  actual  and  complete  data  for  both  nodes  can  be 
found  in  Appendix  E. 


TABLE  4.4  PR-1592  RESONANT  TRACKING  DATA  FOR  THE 

FIRST  TORSIONAL  NODE. 


FREQUENCY 
f  <HS) 

SHEAR  NODULUS 

G  (MPa) 

IN-PHASE  VOLTAGE 
V  (DCF) 

-14.57 

733.56 

39.53 

-0.2046 

-  5.18 

580.63 

24.77 

-0.2092 

5.36 

437.20 

14.04 

-0.2261 

15.35 

338.70 

8.43 

-0.2951 

25.26 

278.03 

5.68 

-0.3870 

35.19 

237.21 

4.13 

-0.4877 

45.24 

211.77 

3.29 

-0.6631 

55.28 

203.17 

2.86 

-0.8283 

65.12 

187.85 

2.59 

-1.0252 

73.68 

182.11 

2.44 

-1.2440 

83.72 

178.52 

2.30 

-1.4435 

£0 


TAB LI  4.5  PR-1592  RESONANT  TRACKING  DATA  FOR  THE 

SECOND  TORSIONAL  NODE. 


FREQUENCY 

SHEAR  MODULUS 

IN-PHASE  VOLTAGE 

T  (#C) 

f  (HZ) 

G  (Mpa) 

V  (DCV) 

-13.77 

1479.76 

1447.1 

-0.2046 

-  5.24 

1184.13 

927.0 

-0.1872 

5.34 

877.51 

509.1 

-0.2001 

15.26 

678.94 

304.8 

-0.2552 

25.17 

553.51 

202.6 

-0.3495 

35.14 

477.74 

150.9 

-0.4544 

45.04 

429.38 

121.8 

-0.5984 

54.97 

396.75 

104.1 

-0.7859 

65.42 

376.62 

93.8 

-1.0044 

75.19 

363.87 

87.2 

-1.2205 

C.  POLYMETHYL  METHACRYLATE  (PMMA) 

Polymethyl  Methacrylate  (PMMA)  is  a  viscoelastic 
material  of  moderate  to  low  molecular  weight.  The  PMMA 
proved  to  be  relatively  stable  within  the  normal  operating 
temperatures  (-15°C  to  85°C)  of  the  unmodified  test  chamber. 
The  original  PMMA  rod,  however,  melted  and  the  transducer 
coils  drooped  out  of  the  magnetic  field  during  a  high 
temperature  test  at  about  90°C.  This  is  clearly  an 
indication  that  the  PMMA  has  a  relatively  rapid  glass  to 
viscous  liquid  transition,  probably  associated  with  the 
material's  low  molecular  weight,  as  explained  earlier  in 
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Chapter  II ,  Section  A.  A  second  sample  was  then  Bade  which 
was  used  for  the  remainder  of  the  test. 

1.  Material  Properties  at  Room  Temperature 

Tables  4.6  and  4.7  show  the  results  for  the  room 
temperature  evaluation  of  the  torsional  and  flexural  modes 
of  the  PMMA  sample.  Figures  4.9  and  4.10  are  the 
corresponding  plots  of  the  frequency  response  for  the 
torsional  and  flexural  modes,  respectively.  Figure  4.9 
shows  the  first  four  harmonically-related  resonances  and  the 
corresponding  pole  zero  data  for  each  resonance  of  the 
torsional  mode.  Figure  4.10  shows  the  first  three 
dispersive  resonances  of  the  flexural  mode  for  the  same 
sample. 

TABLE  4.6  PMMA  SHEAR  MODULUS  RESULTS  (T-23.3eC). 


f 

f/n 

Q 

Q 

E(Leff) 

(Hz) 

(HZ) 

(MEAS) 

(CURVE) 

(Gpa) 

2027.0 

2027.0 

17.29 

17.23 

1.862 

4093.3 

2046.6 

20.71 

20.09 

1.898 

6298.4 

2099.5 

21.66 

22.24 

1.997 

8378.9 

2094.7 

22.53 

23.24 

1.988 
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TAIL I  4.7  PUNA  YOUNG'S  MODULUS  RESULTS  (T-21.4°C) . 


f 

f/n* 

Q 

Q 

E(Leff 

(Hr) 

(Hr) 

(3db) 

(CURVE) 

(Gpa) 

241.4 

26.623 

11.28 

9.27 

4.824 

672.6 

26.910 

18.23 

13.80 

4.928 

1321.9 

26.978 

12.72 

11.01 

4.953 

The  density  for  the  PMMA  sample  vas  1186.43  kg/m*. 
The  effective  length  for  the  torsional  mode  vas  0.309  m  and 
for  the  flexural  mode  vas  0.305  m.  The  dynamic  Poisson's 
ratio  (p)  vas  approximated  for  the  PMMA  sample  to  be  .330 
(>*■£/ 2G-1)  uy  using  the  shear  modulus  for  the  first 
torsional  node  and  the  Young's  modulus  for  the  third 
flexural  mode.  This  is  only  an  approximation  because  the 
measurement  temperatures  differ  slightly. 


72 


FREQ  RESP 


Curvo  Fit. 
gglup _ £ng _ Z+ra*. 


POLES  12  ZEROS  14 


i  -ssi.  see 

-3.  13432K 

2  350.  082 

1.  SS33K*J  1.  S2728K 

3  e.  40785W 

-2.  31817W±J  3.  38455* 

4  -30.  8146  ±J  2.  OZBISk 

2.  12288K±J  4.  83S45K 

3  -101.848  *J  4.  OSlSBk 

-1.  82S78K*J  7.  41020* 

8  -141.  387  ±J  8.2888* 

1.  I8003K±J  7.  8B1S8K 

7  -180.242  ±J  8.  378S8K 

T 1  me  da  1  ay*  O.  O  S  Co  1  n«- 

-133.0/4  Scale-  1.0 

Figure  4.9  PMMA  frequency  response  curves  and  pole  zero 
data  for  the  torsional  inodes  T*23.3°C. 
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FREQ  RESP 


Figure  4.X0  PMMA  frequency  response  curves  and  pole  zero 
data  for  the  flexural  nodes  at  T«21.4#C. 
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2.  Etsoaanot  Track lag  Evaluation  of  PMMA 

The  PMMA  sample  showed  no  glass  transition  within 
the  temperature  range  tested  as  can  be  seen  in  Figures  4.11 
and  4.12  for  the  first  torsional  and  the  first  flexural 
modes  of  the  PMMA  sample.  The  reference  temperature  used  to 
derive  the  reduced  frequency  is  211  K  as  shown  previously  in 
Equation  51.  The  in-phase  voltage  for  the  PMMA  sample  was 
calibrated  by  using  the  value  for  the  Q  obtained  at  room 
temperature  for  the  first  torsional  and  flexural  modes.  The 
proportionality  constant  (A)  obtained  from  these 
calibrations  were  A=-0.0743  and  A*-0.1759,  respectively. 


IE-13  IE-12  IE-11  IE-10  IE-9  IE-8 

REDUCED  FREQUENCY  (Hz) 

Figure  4.11  Master  curve  for  the  first  torsional  mode 
of  PMMA. 
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REDUCED  FREQUENCY  (Hz) 


Figure  4.12  Master  curve  for  the  first  flexural  mode 
of  PMMA. 

The  frequency  and  temperature  dependence  of  the 
shear  and  Young's  moduli  are  plotted  three  dimensionally  in 
Figures  4.13  and  4.14,  respectively. 
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Figure  4.14  Young's  modulus  as  a  function  of  frequency 
and  temperature  for  PMMA. 


Tables  4.8  and  4.9  are  a  summary  of  the  data 
collected  for  the  first  torsional  and  flexural  modes  of  the 
PMMA  sample.  The  actual  and  complete  data  collected  for  the 


PMMA  sample  is  found  in  Appendix  F. 


TABLK  4.8  PMMA  RESONANT  TRACKING  DATA  FOR  THE  FIRST 
TORSIONAL  MODE. 


TEMPERATURE 
T  (°C) 

FREQUENCY 
f  (HZ) 

SHEAR  MODULUS 

G  (GPa) 

IN-PHASE  VOLTAGE 
V  (DCV) 

-13.58 

2170.23 

2.13 

-0.1862 

-  5.36 

2141.82 

2.08 

-0.1777 

5.39 

2101.45 

2.00 

-0.1562 

15.30 

2031.55 

1.92 

-0.1328 

25.10 

2019.47 

1.84 

-0.1120 

35.67 

1972.06 

1.76 

-0.0878 

44.95 

1911.01 

1.65 

-0.0757 

54.98 

1848.78 

1.55 

-0.0572 

65.31 

1804.99 

1.47 

-0.0318 
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1 1BLI  4.9  PMMA  RESONANT  TRACKING  DATA  FOR  THE  FIR8T 
FLEXURAL  MODE. 


Ev.-kWy 

FREQUENCY 

YOUNG ' S  MODULUS 

IN-PHA8E  VOLTAGE 

T  (*C) 

f  (HZ) 

E  (GPa) 

V  (DCV) 

-14.52 

271.50 

6.09 

-0.4797 

-  5.00 

266.72 

5.88 

-0.4424 

5.42 

260.86 

5.62 

-0.3713 

15.40 

254.96 

5.37 

-0.3062 

25.36 

248.79 

5.12 

-0.2475 

35.35 

241.26 

4.81 

-0.2104 

45.35 

233.92 

4.52 

-0.1741 

54.75 

228.34 

4.31 

-0.1054 

65.25 

226.86 

4.25 

-0.0612 

75.10 

243.99 

4.92 

-0.0304 

D.  POLYCARBONATE 

Poly carbonate  is  a  Material  which  had  not  previously 
been  tested  with  free-free  bar  technique  used  in  this 
thesis.  Polycarbonate  is  a  candidate  material  for  use  in 
hydrophone  applications,  thus  the  complex  moduli  and  loss 
properties  of  this  material  are  of  significance. 

1.  Material  Properties  at  Room  Temperature 

Tables  4.10  through  4.12  show  the  room  temperature 
material  properties  for  the  three  vibrational  modes  of  the 
polycarbonate  sample.  Figures  4.15  through  4.17  are  the 
corresponding  plots  of  the  frequency  response  and  pole-zero 


information.  Tho  rosults  of  a  hollow  polycarbonate  sample 

rod  tested  at  room  temperature  are  in  Appendix  6. 

TABLE  4.10  POLYCARBONATE  SHEAR  MODULUS  RESULTS  FOR 
TORSIONAL  MODES  AT  T=23.0°C. 


f 

(HI) 

f/n 

(HI) 

Q 

(3dB) 

Q 

(CURVE) 

G(Leff) 

(GPa) 

1180 

1180 

53.6 

73.31 

0.908 

2377 

1188.5 

69.91 

76.85 

0.918 

3560 

1186.7 

86.8 

88.36 

0.913 

4792 

1198 

94.96 

77.88 

0.929 

5983 

1196.6 

67.22 

80.60 

0.927 

7215 

1202.5 

72.88 

73.91 

0.936 

8416 

1202.3 

84.16 

72.81 

0.936 

TABLE  4.11  POLYCARBONATE  YOUNG'S  MODULUS  RESULTS  FOR 
THE  FLEXURAL  MODES  AT  T=23.0°C. 


f 

(HI) 

f/n* 

(HI) 

Q 

(3dB) 

Q 

(CURVE) 

E(Leff) 

(GPa) 

123.5 

13.62 

47.23 

40.84 

2.498 

344 

13.763 

63.59 

114.16 

2.551 

669.5 

13.663 

114.84 

119.48 

2.514 

TABLE  4.12  POLYCARBONATE  YOUNG'S  MODULUS  RESULTS  FOR 
THE  LONGITUDINAL  MODES  AT  T=23.5°C. 


f 

f/n 

Q 

Q 

B(Leff) 

(HI) 

(HD 

(3dB) 

(CURVE) 

(GPa) 

2000 

2000 

50 

38.15 

2.534 

4010 

2005 

12.94 

53.47 

2.547 

6000 

2000 

32.96 

81.96 

2.534 
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Figure  4. 1C  Polycarbonate  frequency  response  curves  and 
pole  zero  information  for  the  flexural 
modes  at  T»23.0#C. 


82 


;q 


-sm.  oi 


-eo.  a  I — 

-ea.  a  I 

Fxd  Y  IK 


Cur' vc*  Fit. 


POL.ES 


ZEROS 


-3.  9B872K 

-9.  S6826K 

-26.0098  *J  1.98417K 

1 5.  0463 

-37.  3778  *J  3.  99707K 

-54.  4383  *J  2.  D48S4K 

-38.  67 *7  ±J  6.  Ol 129K 

-18.  1307  ±J  3-  9S839K 

—SO.  3892  *J  8.  03S6K 

-SO.  1063  ±J  8.  03087K 

-13.  2737K*J  IS.  0744K 

-42.3628  *J  8.  01908K 

Tima  daloy"  O. O  S  Cotn» 

493.1  Seala-  1.0 

Figure  4.17  Polycarbonate  frequency  response  curves  and 
pole  zero  information  for  the  flexural 
modes  at  T=23.5°C. 


The  dynamic  Poisson's  ratio  ( v)  was  approximated  for 
the  polycarbonate  sample  to  be  .366  (p-E/2G-l)  by  using  the 
shear  modulus  for  the  first  torsional  mode  and  the  Young's 
modulus  for  the  third  flexural  mode.  This  is  only  an 
approximation  due  to  the  fact  that  the  resonance  frequencies 
involved  differ  slightly. 

2.  Resonance  Tracking  Evaluation  of  Polycarbonate 

Figures  4.18  and  4.19  are  plots  of  the  shear  moduli 
of  the  first  torsional  and  Young's  first  flexural  modes  of 
the  polycarbonate  sample  verses  reduced  frequency.  At  the 
time  the  data  was  taken ,  a  WLF  equation  for  the 
polycarbonate  sample  could  not  be  found.  The  closest  match 
available  was  that  which  was  used  for  the  PMMA  sample. 
Accordingly,  the  same  WLf  equation  was  used  for  the 
polycarbonate  sample.  The  reference  temperature  used  to 
derive  the  reduced  frequency  is  211  X. 

The  in-phase  voltage  for  the  polycarbonate  sample  was 
calibrated  by  using  the  value  for  the  Q  obtained  at  room 
temperature  for  the  first  torsional  and  flexural  modes.  The 
proportionality  constant  (A)  obtained  from  the  room 
temperature  evaluations  of  the  polycarbonate  sample  were 
A»- 0.2811  and  A*-0.7532,  respectively.  The  in-phase  voltage 
data  for  both  figures  4.18  and  4.19  has  only  been  verified 
at  one  temperature .  While  this  data  is  thought  to  be 
correct  there  is  at  this  time  no  verification  of  the  data. 
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Figure  4 .IS  Master  curve  for  the  first  torsional  node 
of  polycarbonate. 
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Figure  4.19  Master  curve  for  the  first  flexural  node 
of  polycarbonate . 

The  frequency  and  tenperature  dependence  of  the 
shear  and  Young's  noduli  are  plotted  three  dinensionally  in 
Figures  4.20  and  4.21,  respectively. 
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Tables  4.13  through  4.14  are  a  summary  of  the  data 
collected  for  the  first  torsional  and  flexural  nodes  of  the 
Polycarbonate  sample.  The  actual  data  recorded  for  the 
polycarbonate  sample  can  be  found  in  Appendix  G. 

TABLE  4.13  POLYCARBONATE  RESONANT  TRACKING  DATA  FOR  THE 
FIRST  TORSIONAL  MODE. 


TEMPERATURE 
T  (aC) 

FREQUENCY 
f  (HZ) 

SHEAR  MODULUS 

G  (GPa) 

IN-PHASE  VOLTAGE 
V  (DCV) 

-14.32 

1225.38 

0.972 

-0.1338 

-  5.21 

1208.87 

0.964 

-0.1237 

5.30 

1195.48 

0.926 

-0.1670 

15.30 

1186.29 

0.911 

-0.1793 

25.23 

1177.01 

0.897 

-0.2325 

35.29 

1167.37 

0.883 

-0.3043 

45.32 

1157.42 

0.868 

-0.3185 

55.33 

1150.67 

0.858 

-0.3588 

65.04 

1138.74 

0.840 

-0.3769 

74.84 

1126.85 

0.822 

-0.3784 

83.39 

1120.66 

0.813 

-0.3867 
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TABLE  4.14  POLYCARBONATE  RESONANT  TRACKING  DATA  FOR  THE 
FIRST  FLEXURAL  MODE. 


TEMPERATURE 
T  (*C) 

FREQUENCY 
t  (HZ) 

YOUNGS  MODULUS 
E  (GPa) 

IN-PHASE  VOLTAGE 
V  (DCV) 

-13.43 

127.38 

2.64 

-0.2143 

-  5.17 

125.77 

2.59 

-0.2572 

5.03 

124  ,44 

2.54 

-0.3070 

14.86 

123.46 

2.50 

-0.3807 

25.47 

122.60 

2.46 

-0.4606 

35.45 

121.78 

2.44 

-0.4874 

45.41 

120.79 

2.39 

-0.4830 

55.37 

120.52 

2.38 

-0.4647 

65.09 

119.91 

2.35 

-0.4574 

74.96 

118.93 

2.32 

-0.4740 

B.  TRANSFER  FUNCTION  OF  A  HIGHLY  DAMPED  RLC  CIRCUIT 

The  measurement  of  very  compliant  and  lossy  materials 
(very  low  Q)  poses  special  problems  for  this  measurement 
method.  The  frequency  response  for  this  material  at  room 
temperature  proved  to  be  virtually  flat.  This  is  an 
interesting  problem  from  a  resonance-based  measurement 
technique  perspective.  What  is  and  how  can  the  Q  actually 
be  measured  if  the  frequency  response  at  room  temperature  is 
essentially  flat?  How  can  the  resonance  be  tracked  if  it 
can't  first  be  detected?  These  are  just  two  of  the 
questions  which  had  to  be  answered. 


The  investigation  of  these  types  of  materials  began  by 
returning  to  the  basics.  From  the  previous  discussions  is 


Chapter  IX,  there  are  several  ways  to  approach  these 
problems.  The  first  problem  to  be  resolved  was  whether  or 
not  the  resonance  and  the  quality  factor  for  a  material  with 
low  Q  could  actually  be  measured  accurately  with  the 
equipment  available.  This  was  accomplished  by  simulating  the 
mechanical  system  with  a  simple  rlc  circuit.  The  inductance 
(L)  and  the  capacitance  (C)  of  the  circuit  are  analogous  to 
the  compliance  and  inertiance  of  the  system.  The  resistance 
is  analogous  to  the  resistance  or  real  part  of  the 
impedance.  The  resonant  frequency  can  be  obtained  from  the 
impedance 


Z2  -  R2  +  (joL—_ yi— ) 2  . 

JoC 


(61) 


The  resonant  frequency  occurs  when  the  impedance  reaches  its 


lowest  value 


The  Q  can  be  expressed  as. 


R 


The  results  of  this  experiment  are  summarized  in  Table 
4.15  below.  The  curves  for  both  the  phase  and  pole  zero 
methods  are  included  in  Appendix  H. 
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The  results  of  this  experiment  ere  summarized  in  Table 
4.15  below.  The  curves  for  both  the  phase  and  pole  zero 
methods  are  included  in  Appendix  H. 


TABLS  4. IS  RLC  EXPERIMENT  RESULTS. 


RES* 

(n) 

Fo (THEORY) 
(HZ) 

Fo( CURVE) 
(HZ) 

FO (PHASE) 
(Hz) 

Q (THEORY) 

Q (CURVE) 

21 

4654.9 

4097 

4072 

13.97 

4.34 

70 

N 

4194 

4147 

4.22 

2.48 

120 

N 

4299 

4233 

2.46 

1.71 

170 

n 

4381 

4307 

1.74 

1.82 

220 

N 

4439 

4357 

1.48 

1.104 

260 

N 

4487 

4394 

1.13 

0.97 

370 

If 

4590 

4456 

0.80 

0.70 
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V.  CONCLUSIONS  AMD  RECOMMENDATIONS 


A.  CONCLUSIONS 

The  main  focus  of  this  effort  was  to  define  in  detail 
the  "free- free"  bar  technique  as  it  pertains  to  the 
measurement  of  the  dynamic  moduli  and  loss  properties.  This 
research  encompasses  the  theory,  accuracy,  limitations,  and 
applications  of  the  "free- free"  bar  measurement  technique. 
This  research  addressed  the  inherent  advantages  and 
limitations  associated  with  using  the  PLL  resonance  tracking 
method  and  the  transfer  function  method.  In  the  transfer 
function  method,  the  parameters  of  a  computer  generated 
curve  fit  was  used  to  determine  the  resonance  frequency  and 
Q  of  the  sample.  It  was  pointed  out  in  this  research  that 
because  the  receive  transducer's  output  is  proportional  to 
velocity  as  opposed  to  displacement,  the  maximum  amplitude 
could  be  used  to  determine  the  resonance  frequency  of  the 
the  system.  This  is  especially  critical  in  the  application 
of  the  PLL  resonance  tracking  in  that  the  maximized  in-phase 
voltage  is  related  to  the  maximum  velocity  of  the  sample  rod 
and  is  used  to  determine  the  loss  properties  of  the  material 
as  a  function  of  frequency  and  temperature.  An 
investigation  into  the  ability  of  this  technique  to  measure 
materials  with  very  low  Q  was  also  initiated. 
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An  advantage  of  the  transfer  function  technique  is  that 
the  data  synthesis  performed  by  the  curve  fit  actually  tabes 
into  account  the  effects  of  overlapping  modes  for  any  given 
resonance.  Therefore,  the  pole  zero  curve  fit  appears  to  be 
the  more  accurate  way  of  determining  the  actual  Q  and  the 
resonant  frequency  of  a  given  material.  In  addition,  the 
spectrum  analyzer  is  capable  of  providing  information  on  as 
many  resonances  as  can  be  excited  for  each  mode  at  the  same 
time  in  contrast  to  the  frequency  tracking  in  the  PLL 
method.  The  disadvantage  of  using  the  spectrum  analyzer  is 
that  the  data  acquisition  was  time  consuming  since  an 
automated  data  acquisition  system  has  yet  to  be  developed. 

The  two  channel  PLL  provides  continuous  tracking  of  the 
resonance  frequency  of  any  given  mode  continuously  in  time. 
The  in-phase  output  is  proportional  to  the  maximum  velocity 
response  of  the  system.  This  ensures  that  the  damping 
properties  can  also  be  measured.  Therefore,  since  the 
resonance  of  the  system  occurs  at  the  maximum  amplitude  of 
the  velocity  response,  the  PLL  provides  an  accurate  method 
for  determining  the  actual  temperature  and  frequency 
dependence  of  both  the  dynamic  moduli  and  the  loss  tangent. 
Another  advantage  of  using  the  resonant  tracking  technique 
is  the  data  acquired  using  this  system  is  easily  manipulated 
by  outputting  the  data  into  a  tab  deliminated  format  and 
transferring  to  a  Macintosh  computer. 


92 


On*  disadvantage  with  this  technique  is  that  at  present 
the  systen  is  only  capable  of  tracking  one  resonance  at  a 
tine.  The  in-phase  voltage  still  has  to  be  calibrated  by 
direct  measurement  of  the  Q.  The  set  up  and  operation  time 
required  to  track  the  resonance  over  a  temperature  range  of 
-15°C  to  85 °C  for  a  particular  mode  in  this  research  was 
approximately  six  hours. 

The  PR1592  sample  shoved  a  definite  glass  to  rubber 
transition  within  the  -15 *C  to  85 °C  temperature  range  of 
this  test.  The  independent  measurement  of  the  q  at  various 
temperature  utilising  the  data  synthesis  technique  involving 
the  spectrum  analyser  provided  verification  of  the  PLL  data 
and  important  insights  into  the  differences  between  the  two 
measurement  techniques.  neither  the  PMMA  nor  the 
polycarbonate  samples  showed  any  visible  transitions  within 
temperature  range  of  the  apparatus.  Compliant  materials 
that  show  little  temperature  and  frequency  dependence  are 
candidates  for  hydrophone  applications. 

The  author  was  responsible  for  modifications  to  the 
environmental  chamber  encorporating  the  use  of  CO2  as  a 
coolant  to  decrease  the  low  temperature  capabilities  down  to 
-85 *C.  The  increased  temperature  range  will  allow 
determination  of  the  glass  transition  temperatures  of  highly 
damped  materials  and  in  tracking  the  resonance  of  materials 
through  the  glass  transition  region. 
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The  results  of  the  RLC  circuit  simulation  of  highly 
dampsd  (low  Q)  systsms  prove  that  it  is  possible  to 
electronically  measure  both  the  Q  and  the  resonance  of  such 
materials  with  the  existing  equipment.  The  apparent 
differences  between  the  theoretical  values  and  the  measured 
values  is  a  direct  result  of  the  influence  of  the  impedance 
associated  with  the  inductor.  As  can  be  readily  seen  from 
Table  4.15,  the  theoretical  values  and  the  measured  values 
match  more  closely  as  the  resistance  is  increased  and  the 
importance  of  the  inductor's  impedance  decreases.  In  spite 
of  this  flaw,  the  LCR  circuit  did  prove  the  basic  point  that 
the  spectrum  analyzer  is  capable  of  measuring  very  low 
values  for  Q  at  least  for  a  single  degree  of  freedom  system. 

These  facts  coupled  with  the  new  low  temperature 
capabilities  provide  the  basics  needed  to  enable  the  highly 
damped  materials  to  be  cooled  and  characterized  below  their 
glass  transition  temperatures  in  the  region  where  the 
dynamic  moduli  is  at  a  maximum  for  that  material.  This  is 
significant  in  that  in  this  region  the  resonance  frequency 
for  that  material  is  also  a  maximum.  Therefore,  the 
resonance  frequency  should  be  readily  apparent  and  capable 
of  being  both  detected  and  either  locked  onto  using  the 
phase  locked  loop  technique  or  monitored  using  the  spectrum 
analyzer. 
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B.  BBCCMMBMDATIOM8  TOR  TURTHRR  STUDY 

Questions  surrounding  tbs  sotual  ability  to  par fora 
resonance- based  aaasuraasnts  upon  low  q  materials  still 
remain  unanswered.  Investigations  into  automating  the 
acquisition  of  data  using  the  spectrum  analyser  would  be  a 
very  worthwhile  endeavor.  Using  several  PLL  circuits  in 
parallel  to  measure  more  than  one  mode  at  a  time  may  also 
worth  investigating  but  may  prove  to  be  cost  prohibitive. 
Another  possibility  is  using  the  spectrum  analyser  in 
conjunction  with  the  PLL  in  order  to  quantify  the 
relationship  between  the  Q  and  the  in-phase  voltage. 

The  mechanical  side  of  the  resonant  based  measurements 
for  compliant  materials  is  still  very  much  a  challenge. 
Certainly,  the  question  of  how  to  support  these  soft 
materials  is  a  problem  to  be  investigated  as  was  evidenced 
by  the  first  PMMA  rod  which  melted  causing  the  coils  to 
droop  out  of  the  magnetic  field.  Balancing  this  problem 
while  maintaining  the  freely  suspended  condition  which  now 
exists  could  prove  somewhat  tricky.  There  are  several 
possible  solutions  to  this  problem  including  additional 
mechanical  supports.  Freely  and  vertically  suspending  the 
rod  vertically  using  a  fine  wire  in  order  to  retain  the 
"free-free"  boundary  conditions  is  another  method  to  be 
considered.  Rigidly  mounting  the  sample  on  one  end  and 
mounting  both  the  transmitting  and  receiving  coils  to  the 
other  end  of  the  bar  is  another  viable  alternative.  With  > 
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this  alternative,  however,  cross  talk  created  by  the  induced 
current  in  the  transmitter  and  the  close  proximity  of  the 
receiver  coil  will  have  to  be  addressed.  This  effect  could, 
however,  probably  be  eliminated  by  using  a  third  subtraction 
coil  outside  of  the  magnetic  field.  Another  alternative 
worth  investigating  is  to  use  a  host  bar  with  known 
properties  to  support  the  softer  sample. 


DBLZVATIOV  VOS  DSYSBMIM&TXOS  OF  Q  FROM  SLOPS 


The  phut  relationship  between  applied  force  and  the 
resulting  displacement  of  a  damped  harmonic  oscillator  is 


where 


+  -  tan*1  [ - ^ — 3 

Voting  that  the  derivative  of  aretan  is 


A1 


3tan-1u  _  l  du 
Sr  l+ua  3br 

The  derivative  of  equation  A1  is 


l^r  ia  0<l-»2«i> 

Qd-oli) 


0*<l-«&i)a 


A2 


A3 


0*  (l-<£»l)  a  c  g(l-O^l)  +<*r*l  <20“~1>  3 

o*  d-»ki>  +«Lt  o*  d-«»j)  * 


Simplification  of  this  equation  yields. 


a+  =  0(1 -»Li) 

Q*  (1  2  +*>1*1 


A4 


AS 


When  «=«.,  (i.e.  resonance)  equation  A5  reduces  to 


*— ■* ‘ 20  ' 


Thus, 


afci  =  20  .  0  =  a  aj»£ 
If1***  x  '  0  altf  °  * 


A6 


A7 


where  4  is  in  radians.  This  is  an  exact  result.  [Ref.  11] 
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The  rttonanca  frequency  and  Q  of  the  modes  of  a 
free-free  bar  can  be  found  from  a  pole-zero  plot  of  the 
mechanical  admittance  of  the  oscillator  as  has  been 
described  by  Brown  [Ref.  19].  The  theory  behind  this 
concept  is  derived  based  on  approximating  the  sample  rod  as 
a  damped  driven  harmonic  oscillator.  The  complex  admittance 
for  the  oscillator  is  [Ref.  12] 


_ 2“ _ 

[W(^ca)2  ♦  ~(j«)  + 

M  M 


B4 


where  N  is  the  mass  of  the  oscillator,  R  is  the  mechanical 
resistance,  and  K  is  the  stiffness.  By  substituting  s=ju> 
into  the  equation,  the  admittance  can  be  manipulated  into 
pole - zero  format  as  follows  [R?f.  19] 


Y*ls)  *  3  -  ^  f  3  ) 

M(s2  ♦  £s  +  £)  M  <S~Y)  '  1 

M  M 

where  y  and  y*  are  complex  conjugates.  The  characteristic 
roots  for  the  above  equation  are  [Ref.  10,  19] 


Yi.a 


is 1 


B3 


Therefore,  the  transfer  admittance  function  has  a  pole 
at  y  and  y *  and  a  single  zero  at  «=0 .  The  pole -zero  plot 
is  illustrated  in  Figure  B.l.  The  poles  are  located  on  a 
circle  of  radius  <•>.. 
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Figure  B.l  Pole  zero  plot  of  complex  admittance 
function  [Ref  19]. 

The  quality  factor  can  be  found,  if  the  poles  are  expressed 
as  7  -  a+jb  and  7*  *  a-jb,  as  follows 

Q  «  .  B4 

-2a 

The  resonance  frequency  for  the  oscillator  is  found  by 

«0  -  jb*  *  aJ  -  2>( l+-|p)  -  *>  *  B5 

where  the  second  approximation  is  only  valid  for  low  or 
moderate  damping.  [Ref  19] 
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APPENDIX  D.  RESONANT  TRACKING  COMPUTER  PROGRAM 
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10! 

20! 

30! 

40! 

50! 

60! 

70! 

80! 

90! 

100! 

110! 

120! 

130! 

140! 

150! 

160! 

170! 

180! 

190! 


RT20DLB  :  Resonance  Tracking  System 
Written  by  :  Brian  Beaton, H*ckT«n,  iWusi  A.  8r»on, 
Date  :  18  Mar  90 
Last  modified  :  21  FEB  92 

Features  added  :  1.  End  mass  correction  using  effective 
length.  In  order  to  calculate 
corrected  moduli.  D.  Bartlett, 

21  Feb  92 

2.  Input  and  display  total  mass. 

D.  Bartlett,  21  Feb  92 

3.  Change  ruuning  time  to  program 

to  5  hours  with  measurements  taken 
every  2  minutes.  D.  Bartlett, 

21  Feb  92 


210!  Program  initialization 


230! 

240  MASS  STORAGE  IS  "INTERNAL, 4,0" 

250  RE-STORE  "RT20B:INTERNAL,4,r 

260! 

270  OPTION  BASE  0 

280  DIM  Thermistor(0:240),  Resfreq(0:240),  Time(0:240), 

Inphmag(0:240) 

290  DIM  Mod(0:240),  X(0:240),  Y(0:240),  Alphat(0:240), 

Redfreq(0:240) 

300  DIM  Temp(0:240) 

310  DIM  Gtor(0:240),Elong(0:240),Eflex(0:240) 

320  DIM  Label$(0: 25)  [50] 

330  DIM  Main_title$[50],Sub_title$[50] 

340  DIM  X_axis_name$[50],Y_axis_name$[50] 

350  DIM  Blockl$[50)31ock2S[50],Block3$[50],Block4$[50] 

360  DIM  Cl$[25],C2$[25],C3S[25],C4a$[25],C4b$[25],C5$[25] 

370  DIM  C6$[25],C7$[25],C8$(25],C9${25],C10$[25] 
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380 
390  ! 
400 
410 
420  ! 
430 
440 
450 
460 
470  ! 
480 
490 
500 
510 
520  ! 
530 
540 
550 
560 
570 
580 
590 
600 
610 
620 
630 
640 
650 
660 
670 
680 
690! 


DIM  L1$[50]£2$[50],L3$[50]X4${50] 

acarS=CHRJ(255)&CHR${75) 

Home$=CHRS(255)&CHR$(75) 


!  CLEAR  SCRN  key 
!  HOME  key 


Flag=0 

Tmode=l 

Lmode=l 

Fmode=3.0112 


!Flag=l  to  indicate  editing  data 
!  FIRST  torsional  mode 
!  FIRST  longitudinal  mode 
!  FIRST  flexural  mode 


Sample_rate=2 
Arg  1  =60/Sample_rate 
Pi=PI 

Default_grid$="Partial" 


!  1  sample  every  2  minutes 
!  samples  per  hour 
!  3.14159265 
!  Full,  Partial  or  No 


Label$(  1  )="Resonance  Frequency  vs  Time"!  Gra. 4  Main  Titles 
Label$(2)="Temperature  was  varied"!  Sub  Tilk;  Axes  names 
Label$(3)="Time,  hours" 

Label$(4)="Resonance  frequency,  Hz" 

Label$(5)="Temperature,  C" 

Label$(6)="Resonance  Frequency  vs  Temperature" 
Label$(7)="Temperature  vs  Time" 

LabeI$(8)="Frequency  variation  was  observed" 
Label$(9)="Young's  Modulus  vs  Temperature" 
Label$(10)="Shear  Modulus  vs  Temperature" 
LabeI$(lI)="Shear  modulus,  Pa,  xl0A" 

Label$(  12)=" Young’s  modulus,  Pa,  xl0A" 

Label$(13)="In  Phase  Magnitude,  DCV" 

LabelS(14)="In  Phase  Mag  vs  Time" 

Label$(15)="In  Phase  Mag  vs  Temp" 

Label$(16)="In  Phase  Mag  vs  Freq" 


700 

Pen  1=1 

!  White 

Default  colors 

710 

Pen2=2 

!  Red 

720 

Pen3=3 

!  Yellow 

730 

Pen4=4 

!  Green 

740 

Pen5=5 

!  Cyan 

(greenish  blue) 

750 

Pen6=5 

!  Blue 
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780  Options:  !  Allow  user  to  process  existing  data  or  collect  oew  data 


800! 

810  OUTPUT  KBD;Qear$;  !  Gear  the  CRT 

820  OUTPUT  KBD;Home$;  !  Home  display 

830  STATUS  l,9;Screen  !  Get  screen  width 

840  Center=(Screen-28)/2  !  Leading  spaces  for  centering 

850  GCLEAR 

860  PRINTER  IS  CRT  !  Use  CRT  for  displaying  menu 

870  PRINT  TABXY(1,1)  !  Start  at  top  with  blank  line 

880  PRINT  TAB(Center);"Key  Purpose" 

890  PRINT  TAB(Center);" - " 

900  PRINT  TAB(Center);"  0  Process  existing  data" 

910  PRINT  TAB(Center);"  4  Collect  new  data" 

920  PRINT  TAB(Center);"  5  Editing  existing  data" 

930  FOR  Keynumber=0  TO  9  !  Off  all  keys 

940  ON  KEY  Keynumbcr  LABEL  ""  GOSUB  Commentl 

950  NEXT  Keynumber 

960  ON  KEY  0  LABEL  "Existing  Data"  GOTO  Existingjiata  !  Key  0 

970  ON  KEY  4  LABEL  "New  Data”  GOTO  New.data  !  Turn  on  Key  4 

980  ON  KEY  5  LABEL  "Edit  Data"  GOTO  Editing_data  '.Turn  on  Key  5 

990  BlinkO:  WAIT  1 

1000  DISP  "Select  an  option" 

1010  WAIT  1 

1020  DISP 

1030  GOTO  BlinkO 

1040  RETURN 

1050!******************************************************************* 
1060  New.data:  !  Main  program  driver  for  collecting  new  data 
1070!******************************************************************* 


1080 

GOSUB  Clear_keys 

1090 

GOSUB  Init.param 

1100 

GOSUB  Init_periph 

1110 

GOSUB  Input_nm_data 

1120 

GOSUB  Create_new_file 

1130! 

GOSUB  DrawJreqjime 

1140 

GOSUB  Draw_mag_time 

1150 

GOSUB  Take_plot_data 

1160 

GOSUB  Close.file 
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1170  GOSUB  Free_periph 

1180  GOSUB  Postjjroccss 

1190  STOP 


iilil  •  T  *  •  T  r,  *  TT”T  •  »  »  .»(•  »  TTTTT'TIITTITT”I':"rTT 


1210  Existing_data:  !  Main  program  driver  for  process  existing  data 

1230  GOSUB  Clearjceys 

1240  GOSUB  Open_old_file 

1250  GOSUB  Post_process 

1260  STOP 

\210\~*”*™*”**“**™*****™™**^******************************* 

1280  Editing_data:  !  Main  program  driver  to  compute  Redfreq  for  old  data 

1300  Flag=l 

1310  GOSUB  Gearjceys 

1320  GOSUB  Open_old_file 

1330  GOSUB  Calcjredfreq 

1340  GOSUB  Create_new_file 

1350  GOSUB  Qosejile 

1360  STOP 


1380  Clear_keys:  !  "Turns  off  all  softkeys  and  clears  the  screen 


1400  FOR  Keynumber=0  TO  9 

1410  ON  KEY  Keynumber  LABEL  ""  GOSUB  Commentl 

1420  NEXT  Keynumber 

1430  OUTPUT  KBD;Clear$; 

1440  RETURN 

1460  Imt_param:  !  User  sets  the  minimum  and  maximum  values  for  initial  plot 


1480  Center=(Screen-46)/2  !  Leading  spaces  for  centering 

1490  PRINT  TABXY(1,1)  !  Start  at  top  with  blank  line 
1500  PRINT  TAB(Center);"Key  Purpose" 

1510  PRINT  TAB(Center);" - " 

1520  PRINT  TAB(Center);"  0  Set  minimum  resonance  frequency" 

1530  PRINT  TAB(Center);"  1  Set  maximum  resonance  frequency" 

1540  PRINT  TAB(Center);"  2  Set  minimum  in  phase  magnitude” 

1550  PRINT  TAB(Center);"  3  Set  maximum  in  phase  magnitude" 
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1560  PRINT  TAB  (Center);’'  5  Set  minimum  temperature  in  ,  C 

1570  PRINT  TAB(Center);"  6  Set  maximum  temperature  in  ,  C 

1580  PRINT  TAB(Center);"  7  Set  the  default  values” 

1590  PRINT  TAB  (Center);"  8  Review  the  assigned  values" 

1600  PRINT  TAB(Center);"  9  To  proceed" 

1610! 

1620  ON  KEY  0  LABEL  "Min  Freq,  Hz"  GOSUB  Minjrequency 

1630  ON  KEY  1  LABEL  "Max  Freq,  Hz"  GOSUB  Maxjirequency 

1640  ON  KEY  2  LABEL  "Min  Mag,  DCVH  GOSUB  Minjnag 

1650  ON  KEY  3  LABEL  "Max  Mag,  DCV"  GOSUB  Maxjnag 

1660  ON  KEY  4  LABEL""  GOSUB  Commentl 
1670  ON  KEY  5  LABEL  "Min  Temp,  C"  GOSUB  Minjemperature 

1680  ON  KEY  6  LABEL  "Max  Temp,  C  GOSUB  Maxjemperature 

1690  ON  KEY  7  LABEL  "Default  Values"  GOSUB  Default_values 

1700  ON  KEY  8  LABEL  "Review  values"  GOSUB  Review.values 

1710  ON  KEY  9  LABEL  "Proceed"  GOTO  Moveon 

1720! 

1730  Blinkl:  WATT  1  !  Wait  for  user  selection  and 

1740  DISP  "Select  an  option"  !  then  take  appropriate  action 

1750  WATT  1 

1760  DISP 

1770  GOTO  Blinkl 

1780  Moveon:  GOSUB  Gearjceys 

1790  DISP  "Proceeding 

1800  WAIT  1 

1810  DISP 

1820  RETURN 


1840  Commentl:  !  Alerts  user  when  an  unassigned  soft  key  is  selected 

1860  BEEP  300.1 

1870  DISP  "This  soft  key  is  unassigned" 

1880  WAIT  1 

1890  DISP 

1900  RETURN 

19101******************************************************************* 

1920  Default.values:!  Assigns  default  values  for  minimum  and  maximum  TJF 
1930!*************************  ****************************************** 

1940  User_freq_min= 1 65  0 

1950  Uscr_freq_max=1900 


in 


UserjempjninsO 
Uscr_tcmp_max=25 
DISP  "The  Default  Values  are  Set" 
WATT  1 

GOSUB  Review_values 
RETURN 


30  Minjrequency:  !  Accepts  user  input  for 


L  x  _x  x  x  xj 


hihtUiMM 


frequency 


INPUT  "The  min  freq  for  the  plot,  in  Hz  ?",  UserJreqjnin 
DISP  "The  min  frequency  is  set  at:  ";User  Jreqjnin;"  Hz" 
WATT  1 
DISP 
RETURN 


2110  Maxjrequency:  !  Accepts  user  input  for  maximum  frequency 
21201^w^^**^»***»***********^*^***********»*********»*****»>** 
2130  INPUT  "The  max  fireq  for  die  plot,  in  Hz  ?",User_ffeq_max 

2140  DISP  "The  max  freq  is  set  at*  ";UserJreqjnax;"Hz" 

2150  WATT  1 
2160  DISP 

2170  RETURN 


2190  Min.mag:  !  Accepts  user  input  for  minimum  magnitude 


rAMl 


T  •  *  • T  * TTTTtTTTtT 


INPUT  "The  min  mag  for  die  plot,  in  DCV  ?",Userjnagjnin 
DISP  "The  minimum  magnitude  is  set  at:  ";User_mag_min 
WAIT  1 
DISP 
RETURN 


2270  Maxjmag:  !  Accepts  user  input  for  maximum  magnitude 


ft  Lxxi. 


"T . tttttttttt 


2290  INPUT  "The  max  mag  for  the  plot,  in  DCV  ?",Userjnag_max 

2300  DISP  "The  maximum  magnirnde  is  set  at:  ”;Userjnag_max 

2310  WATT  1 
2320  DISP 

2330  RETURN 

234Qj»»Ma»»«»**  ********************************************************* 


2350  Review_values:  !  Presents  currently  assigned  values  for  review 

2370  DISP  "Minimum  frequency:  ";User_freq_min;"  Hz" 

2380  WAIT  1 

2390  DISP  "  Maximum  frequency:  ";User_freq_max;"  Hz" 

2400  WAIT  1 

2410  DISP  "Minimum  magnitude:  ";User_mag_rain;"  DCV" 

2420  WAIT  1 

2430  DISP  "  Maximum  magnitude:  ";User_magjmax;"  DCV" 

2440  WAIT  1 

2450  DISP  "Minimum  temperature:  ";Userjemp_min;"  C" 

2460  WAIT  1 

2470  DISP  "  Maximum  temperature:  ";User_temp_max;"  C" 

2480  WAIT  1 

2490  DISP 

2500  RETURN 


kirn 


2520  Min_temperature:  !  Accepts  user  input  of  minimum  temperature 


T_t  •  i  T 


2540  INPUT  "Minimum  temperature,  degrees  C  ?",User_temp_min 

2550  DISP  "The  min  temp  is  set  at:  ";Userjemp_min;"Deg  C" 

2560  WAIT  1 

2570  DISP 

2580  RETURN 

2600  Max.temperature:  !  Accepts  user  input  of  maximum  temperature 

2620  INPUT  "Maximum  temperature,  degrees  C  ?",User_temp_max 

2630  DISP  "The  max  temp  is  set  at:  "  ;User_temp_max;"Deg  C" 

2640  WAIT  1 

2650  DISP 

2660  RETURN 

2680  Init_periph:  !  Initializes  the  voltmeter  and  the  frequency  counter 


27001  720  -  HP5316A  Universal  Counter 

2710!  722  -  HP3456A  Digital  Voltmeter 

2720!  724  -  HP3478A  Digital  Multimeter  (added  on  23  Jan  91) 

2730! 

2740 


OUTPUT  720;  "IN” 


!  Default  state 
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2750  OUTPUT  722;"HF4R1M66STG100STT  !  2  wire  ohms,  THMS  deg  C 
2760  !  100  line  cycles  integration 

2770  OUTPUT  724;"F1RAN5T1Z1D1" 

2780  RETURN 


2800  Free_periph:!  Frees  the  voltmeter  and  the  counter  from  the  HPIB  bus 


2820  LOCAL  720 

2830  LOCAL  722 

2840  LOCAL  724 

2850  RETURN 


2870  Input_run_data:  !  Accepts  user  input  of  selected  run  data 

2$ffl[**************+************************************************+*+* 

2890  Center=(Screcn42)/2  !  Leading  spaces  for  centering 

2900  PRINT  TABXY(1,1)  !  Start  at  top  with  blank  line 

2910  PRINT  TAB(Center);"Key  Purpose" 

2920  PRINT  TAB(Center);" - B 

2930  PRINT  TAB(Center);"  0  Enter  rod  identification" 

2940  PRINT  TAB(Center);"  1  Enter  the  run  number  for  this  mode" 

2950  PRINT  TAB(Center);"  2  Eater  the  mode  for  this  run" 

2960  PRINT  TAB(Center);”  3  Enter  the  date  for  this  run" 

2961  PRINT  TAB(Center);"  4  Enter  the  total  mass  of  rod  &  coils" 

2970  PRINT  TAB(Center);"  5  Enter  the  mass  of  the  rod  only" 

2980  PRINT  TAB(Center);"  6  Enter  the  length  for  this  rod" 

2990  PRINT  TAB(Center);"  7  Enter  the  diameter  for  this  rod" 
3000  PRINT  TAB  (Center);"  8  Review  the  information  entered" 

3010  PRINT  TAB(Center);"  9  To  proceed" 

3020! 

3030  ON  KEY  0  LABEL  "Filename  (Rod  ID)"  GOSUB  Rod.id 
3040  ON  KEY  1  LABEL  "Run  No."  GOSUB  Run_number 
3050  ON  KEY  2  LABEL  "Mode"  GOSUB  Mode 
3060  ON  KEY  3  LABEL  "Date"  GOSUB  Date 
3070  ON  KEY  4  LABEL  "Total  mass"  GOSUB  Total  jnass 

3080  ON  KEY  5  LABEL  "Rod  Mass"  GOSUB  Rodjnass 

3090  ON  KEY  6  LABEL  "Rod  Length"  GOSUB  Rodjength 

3100  ON  KEY  7  LABEL  "Rod  Diameter"  GOSUB  Rod.diameter 
3110  ON  KEY  8  LABEL  "Review  entries"  GOSUB  Review_entrys 
3120  ON  KEY  9  LABEL  "Proceed"  GOTO  Onward 
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3130  Blink2:  WAIT  1 

3140  D1SP  "Select  an  option  or  proceed" 

3150  WATT  1 

3160  DISP 

3170  GOTO  Blink2 

3180  Onward:  Volume=(Length/100.0)*Pi*.25*(Diameter/100.0)A2 

3190  Densitys(Mass/1000)/Volume 

3191  M=Totalmass-Mass 

3200  GOSUB  Clear.keys 

3210  DISP  "Proceeding " 

3220  WAIT  1 

3230  DISP 

3240  RETURN 


3260  Rod_id:  !  Accepts  user  input  of  filename  &  rod  identification 


3280  INPUT  "Filename  (Rod  ID)  block  label  (i.e.  £0*4)  ?",Blockl$ 

3290  WHILE  LEN(Blockl$)>12 

3300  DISP  "Limit  Rod  Identification  to  12  characters” 

3310  BEEP  300,. 1 

3320  WATT  1 

3330  DISP 

3340  INPUT  "Rod  identification  block  label  (i.e.  ECP4) 

?"31ockl$ 

3350  END  WHILE 

3360  DISP  "The  Rod  ID  is  set  at  ";BIockl$ 

3370  WAIT  1 

3380  DISP 

3390  RETURN 


3410  Run_numben  !  Accepts  user  input  of  run  number 


3430  INPUT  "Run  number  for  this  mode  (i.e.  4)  ?"31ock2$ 

3440  DISP  "The  Run  number  is  set  at:  ";Block2$ 

3450  WATT  1 
3460  DISP 

3470  RETURN 


3490  Mode:  !  Accepts  user  input  of  the  mode 
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3510 

3520 

3530 

3540 

3550 

3560 

3570 

3580 

3590 

3600 

3610 

3620 

3630 

3640 


INPUT  "Mode  (Tor»ooal,Longjtudinal,Flexural)  ?"JBlock3$ 
SELECT  UPC$(Block35) 

CASE  "FLEXURAL", "TORSIONAL", "LONGTTUDINAL" 

DISP  "The  mode  is  set  as:  ";Block3$ 

WAIT  1 
DISP 

CASE  ELSE 
BEEP  300, .1 

DISP  "Choices  are:  Torsional  .Longitudinal,  Flexural" 
WAIT  1 
DISP 

GOTO  Mode 
END  SELECT 
RETURN 


3660  Date:  !  Accepts  user  input  of  run  date 

36701******************************  v**^******************************* 
3680  INPUT  "The  date  for  this  sample  (i.e.  10  SEP  89)  ?"JBlock4$ 

3690  WHILE  LEN(Block4S>9 

3700  DISP  "Limit  date  entry  to  9  characters" 

3710  BEEP  300.1 

3720  WATT  1 

3730  DISP 

3740  INPUT  "Enter  the  date  (i.e.10  SEP  89)  ?"3Iock4$ 

3750  END  WHILE 

3760  DISP  "The  Date  is  set  au  ";Block4$ 

3770  WAIT  1 

3780  DISP 

3790  RETURN 

37911************  ***************************************  **************** 
3792  Total.mass:  !  Accepts  user  input  of  total  sample  mass 
3793!******************************************************************* 

3794  INPUT  "The  total  mass  for  sample  rod  (grams)  ?",Totalmass 

3795  DISP  "The  total  mass  is  set  at:  ";Totalmass;"  grams" 

3796  WAIT  1 

3797  DISP 

3798  RETURN 

3800!******************************************************************* 

3810  Rodjnass:  !  Accepts  user  input  of  the  rod  mass 
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3830  INPUT  "The  mass  for  this  rod  (units:  grams)  ?",Mass 
3840  DISP  The  mass  is  set  at:  ";Mass;"  grams" 

3850  WAIT  1 
3860  DISP 

3870  RETURN 

3890  Rodjength:  !  Accepts  user  input  of  the  rod  length 


y/wi 

3910  INPUT  "The  length  for  this  rod  (units:  cm)  ?",Length 

3920  DISP  "The  length  is  set  at:  ";Length;"centimeters" 
3930  WAIT  1 

3940  DISP 

3950  RETURN 


3970  Rod_diameter:  !  Accepts  user  input  of  the  rod  diameter 


3990  INPUT  "The  diameter  for  this  rod  (units:  cm)  ?"JDiameter 

4000  DISP  "The  diameter  is  set  at:  ";Diameter,"  centimeters" 

4010  WAIT  1 
4020  DISP 

4030  RETURN 

4QQfi\***+*+*+**++++********++*+**+***+*+i,****ir***+*************+******+* 

4050  Review_entrys:  !  Presents  currently  assigned  values  for  review 
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4190 

4200 

4210 

4220 


4240 

4250!* 

4260 

4270 

4280 

4290 

4300 

4310 

4320 

4330 

4340 

4350 


DISP  "Diameter  "-.Diameter,"  centimeters" 
WAIT  1 
DISP 
RETURN 


*** 

Draw.mag.time:  !  Produces  Magnitude  vs  Time  graph  w/o  curve 


Main_title$=Label$(14) 

Sub_title$=Label$(2) 

X_axis_name$=Label$(3) 

Y.axis_name$=Label$(13) 

Xmin=0 

Xmax=5 

YminsActual.mag.min 
YmaxsActual.mag.max 
GOSIJB  Generic.plot 
RETURN 


4370  Draw.freq.time:  !  Produces  Frequency  vs  Time  graph  w/o  curve 


4390 

4400 

4410 

4420 

4430 

4440 

4450 

4460 

4470 

4480 


Main_tirleS=LabelS(l ) 
Sub.title$=Label$(2) 
X_axis_name$=Label$(3) 
Y_axis_name$=Label${4) 
XminsO 
Xmax=5 

Ymin=Actual_freq_min 
Ymax=Actual_fireq_max 
GOSUB  Generic_plot 
RETURN 


!  Actual  minimum  frequency 
!  Actual  maximum  frequency 


4500  Young_mod:  !  Compute  Young's  modulus  in  flexure  modes 

4520  D=Diameter/100 

4530  LsLeagth/100 

4531  LeffsL*(l+(M/Mass)) 

4540  Arg6=Density*((32*LeffA2V(Pi*D*FmodeA2))A2 

4550  Eflex(l)=Arg6*Resfreq(I)A2 

4560  RETURN 

4570J******************************************************************* 
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4580  Shear.mod:  !  Compute  Shear  Modulus 


4590!******************************************************************* 


4600  D=Diameter/100 

4610  L=Length/100 

4611  Lefs=L*(l+<2*M/Mass)) 

4620  Arg5=(Density*4*LefsA2)/(TxrtodeA2) 

4630  Gtor(D=(Arg5*Resfreq(J)A2) 

4640  RETURN 


4660  Lyoungjnod:  !  Compute  Young's  modulus  in  longitudinal  mode 

^yQ{**«««*«*****»4i***««* *********************************************** 

4680  D=Diameter/100 

4690  L=Length/100 

4691  Lefl=L*(l+(M/Mass)) 

4700  Arg7=(Density*4*LeflA2)/(LmodeA2) 

4710  Elong(K)=Arg7*Resfreq(K)A2 

4720  RETURN 

47301******************************************************************* 


4740  Compute_redfreq:  !  Compute  redfreq  of  new  data 
4750!******************************************************************* 
4760  Temp(I)=Thermistor(l)+273. 1 5 

4770!  X(I)=-12.9*(Temp(I)-283.15)/(107+Temp(I)-283.15)  !  PR  1592 

4780  X(I)=-21.5*(Temp(I)-21  l)/(43.1+Temp(I)-21 1)  !  Plexi-glass  o 

4790  Y(l)=Xa)*LOG(10) 

4800  Alphat(I)=EXP(Y(D) 

4810  Redfreq(I)=Resfreq(I)*Alphat(I) 

4820  RETURN 

4830!******************************************************************* 


4840  Calc.redfreq:  !  Convert  RT20B  to  RT20C  format 

^^\^m*¥*****i^*^************************  *************************** 


4860  FOR  1=0  TO  5*Argl 

4870  Temp(I)=Thennistor(I)+273. 1 5 

4880  X(I)=-12.9*(Temp(I)-283.15)/(107+Temp(I)-283.15)!PR  1592 

4890  Y(I)=X(I)*LOG(10) 

4900  Alpfaat(I)=EXP(Y(I)) 

4910  Redfreq(I)=Resffeq(I)*Alphat(I) 

4920  NEXT  I 

4930  RETURN 

4940'******************************************************************* 
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4950  Computejnodulus:  !  Computes  appropriate  modulus  based  on  mode 


4970 

4980 

4990 

4991 

4992 

4993 
5000 
5010 
5020 
5030 
5040 
5050 
5060 
5070 
5080 
5090 
5100 
5110 
5120 
5130 
5140 
5150 
5160 
5170 
5180 
5190 
5200 
5210 
5220 
5*30 
5240 
5250 
5260 
5270 
5280 
5290 
5300 
5310 


DISP  "Computing  appropriate  modulus . " 

D=Diameter/100  !  convert  to  meters 

L=Length/100  !  convert  to  meters 

Leff=L*(l+(M/Mass)) 

Le£s=L*(l+(2*M/Mass)) 

Lefl=L*(l+(M/Mass)) 

SELECT  UPC$(Block3$) 

CASE  "FLEXURAL" 

Arg2=Density*((32*LeffA2)/(Pi*D*FmodeA2))A2 
Eflexjnax=Arg2*Actual_freq_maxA2 
Power=LOG(Eflex_max)  DIV  LOG(IO) 

Scale_factor=10APower 
Scaie_eflex_max=Eflex_max/Scale_factor 
Scale_eflex_min=(Arg2*AcrualJreq_minA2)/ScaIe_factor 
FOR  1=0  TO  5*Ar;i 

Eflex(I)=(Arg2*Resfreq(I)A2)/Scale_factor 
NEXT  I 

CASE  "TORSIONAL" 

Arg3=(Density*4*LefsA2)/(TmodeA2) 

Gtor_max= Arg3  *Acrual_freq_maxA2 
Power=LOG(Gtor_max)  DIV  LOG(IO) 

Scale Jactor=10APower 
Scale_gtor_max=Gtor_max/Scale_factor 
Scale_gtor_min=(Arg3*ActuaJ_freq_minA2)/Scale_factor 
FOR  J=0  TO  5*Argl 

Gtor(J)=(Arg3*Resfreq(J)A2)/Scale_factor 
NEXT  J 

CASE  "LONGITUDINAL" 

Arg4=<Density*4*LeflA2)/(LmodeA2) 
Elong_max=Arg4*Acmal_freq_maxA2 
Power=LOG(Elong_raax)  DIV  LOG(IO) 

Scale_factor=l  0APower 
Scale_elongjnax=EIong_max/Scale_factor 
Scale_eIongjnin=(Arg4*Actuai_freq_minA2)/Scale_factor 
FOR  K=0  TO  5*Argl 

Elong(K)=(Arg4*Resfreq(K)A2)/Scale_facior 
NEXT  K 
END  SELECT 
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DISP 

RETURN 


5350  Percent jnoduli:!  Determines  the  percent  change  in  modulus  per  deg  C 

5370  Delta_temp=Actual_temp_max-Actual_temp_min 

5380  Fsum2pt=Actual_freq_max+Actual_freq_min 

5390  Fdiff2pt=Actual_freq_max-Actual_freq_min 

5400  Fminls=Acnialjemp_max*Slope+Intercept 

5410  Fmaxls=Acnial_temp_min*Slope+Intercept 

5420  Fsumls=Fmaxls+Fminls 

5430  Fdiffls=Fmaxls-Fminls 

5440  Prcnt_per_c_2pt=100*4*(Fdiff2pt/Fsum2pt)/Delta_temp  !  2  pt 

5450  Prcnt_per_cJs=100*4*(Fdiffls/Fsumls)/E)eIta_temp  !L  squares 

5460  RETURN 

^flQl******************************************************************* 

5480  Magjemp:  !  Produces  In  phase  mag  vs  Temperature  graph  w/  curve 


****** 


5500 

5510 

5520 

5530 

5540 

5550 

5560 

5570 

5580 

5590 

5600 


Mam_title$=Label$(15) 

Sub_dtleS=Label$(2) 

X_axis_name$=Label$(5) 

Y_axis_name$=Label${  13) 

Xmin=ActuaI_temp_min 
Xmax=Actual_temp_max 
Ymin=Actual_mag_min 
Ymax=Actual_mag_max 
GOSUB  Generic _plot 

CALL  Generic_curve(Pen3,Thermistor(*)tInphmag(*)t5*Argl) 
RETURN 


5620  Magjreq:  IProduces  In  phase  mag  vs  Resonant  frequencies 
56301******************************************************************* 
5640  Main _title$=Label$(l  6) 

5650  Sub_titleS=" " 

5660  X_axis_name$=Label$(4) 

5670  Y_axis_name$=Label${13) 

5680  Xmin=Actual_freq_min 

5690  Xmax=Actuai_freq_max 

5700  Ymin=Actual_mag_min 
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5710  Y  max= Actual_mag_max 

5720  GOSUB  Generic_plot 

5730  CALL  Genericj:urve(Pen3,Resfreq(*),Inphmag(*),5*Argl) 

5740  RETURN 


5760  Shear_temp:  !  Produces  Shear  modulus  vs  Temperature  graph  w /  curve 


5780 

5790 

5800 

5810 

5820 

5830 

5840 

5850 

5860 

5870 

5880 


Main_title$=Labei$(  1 0) 

Sub jitle$=" " 

X_axis_name$=Label$(5) 
Y_axis_name$=Label$(l  1  )&VAU(Power) 


Xmin=Actual_temp_min 
Xmax=ActuaI_temp_max 
Ymin=Scale _gtor_min 
Ymax=Scale_gtor_max 
GOSUB  Generic_plot 

CALL  Generic_curve(Pen3,Thennistor(*),Gtor(*)y5*Argl) 
RETURN 


!  Actual  minimum  temperature 
!  Actual  maximum  temperature 
!  Scaled  minimum  shear  modulus 
!  Scaled  maximum  shear  modulus 


5900  Young_temp:!  Produces  Young's  modulus  vs  Temperature  graph  w /  curve 


5920  Main_titleS=Label$(9) 

5930  Sub_title$="  " 

5940  X_axis_name$=Label$(5) 

5950  Y_axis_name$=Label$(  1 2)&V  AL$(Po  wer) 

5960  Xmin=Actual_temp_min  !  Actual  minimum  temperature 

5970  Xmax=Actualjempjnax  !  Actual  maximum  temperature 

5980  SELECT  UPC$(Block3$) 

5990  CASE  "FLEXURAL" 

6000  Ymin=Scale_eflex_mm  -  !  Scaled  minimum  Y  (flexural) 

6010  Ymax=Scale_eflex_max  !  Scaled  maximum  Y  (flexural) 

6020  GOSUB  Generic_plot 

6030  CALL  Generic_curve(Pen3,Thermistor(*)£flex(*),5*Argl) 

6040  CASE  "LONGITUDINAL" 

6050  Ymin=Scale_elong_min  !  Scaled  min  Y  (longitudinal) 

6060  Ymax=Scale_elong_max  !  Scaled  max  Y  (longitudinal) 

6070  GOSUB  Generic_plot 

6080  CALL  Generic_curve(Pen3  .Thermi stor(  *)  ,Elong(  * )  ,5  *  Arg  1 ) 

6090  END  SELECT 

6100  RETURN 
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61  ioi*****************************^********************  **************** 
6120  Generic_plot:  !  Produces  design  and  layout  of  line  graph 


6140 

6150 

6160 

6170 

6180! 

6190 

6200 

6210 

6220 

6230 

6240 

6250 

6260 

6270! 

6280 

6290 

6300 

6310 

6320 

6330! 

6340 

6350 

6360 

6370 

6380 

6390! 

6400 

6410 

6420 

6430 

6440 

6450 

6460 

6470! 

6480 

6490 


GINIT  !  Initialize  various  graphics  parameters 

GCLEAR  !  Gear  the  graphics  display 

GRAPHICS  ON  !  Turn  the  graphics  display  on 

OUTPUT  KBD;Gear$;  !  Gear  the  CRT 

IF  Plot_device$="Plotter"  THEN  !  Route  to  HP7475A  plotter 
GRAPHICS  OFF 
BEEP  1000.1 

DISP  "Press  CONTINUE  when  the  plotter  is  ready" 

PAUSE 

OUTPUT  705;"IP  1300,1000,9000,6750;"!  Set  scaling  points 
PLOTTER  IS  705,"HPGL" 

END  IF 

Xgdumax=100*MAX(lJRATIO)  !  How  many  gdu's  wide  screen  is 
Ygdumax=100*MAX(l, 1/RATIO)  !  How  many  gdu's  high  screen  is 
CALL  Scale(.02,.98,.10,.98,TopJLeft^Ccntr,Ycntr^Cgdumax,Ygdumax) 

PEN  Pen2 
FRAME 

IF  LEN(Main_title$)>30  THEN  !  Size  main  title 

CALL  Label(4,.6,0,6,Pen2,Xcntr,Top-2,Mam_title$) 

ELSE 

CALL  Label(6,.6,0,6JPen2rXcntr,Top-2Main_title$) 

END  IF 

CALL  Label(4,.6,0,6J>en2rXcntr,.9*Top,Sub_titleS) 

CALL  BlockJnfo(Blockl$,Block2$,Block3$,Block4$,Xgduinax,Pen5,Pen2) 
CALL  Scale(.15,.90,.25,.8,TopJ-eftrXcntr,Ycntr^Cgdumax,Ygdumax) 

CLIP  OFF 

CALL  Labei(4,.6,0,4,Pen5,Xcntr,-l  8,X_axis_name$) 

CALL  Label(4,.6,90,4J>en5  ,-12,Ycntr,Y_axis_name$) 

CLIP  ON 


!  Auto  scale  X  axis 


SELECT  XscaleS 
CASE  "Auto  scale  X" 
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6500 

6510 

6520 

6530 

6540 

6550 

6560 

6570 

6580 

6590 

6600 

6610 

6620 

6630 

6640 

6650! 

6660 

6670 

6680 

6690 

6700 

6710 

6720! 

6730! 

6740 

6750 

6760 

6770 

6780 

6790 

6800 

6810 

6820 

6830 

6840 

6850 

6860 

6870 

6880 


CALL  Xscale(XmmJ(maxJ£minor,Xmajor) 

CASE  "User  scale  X"  !  User  scale  X  axis 

SELECT  PassS 

CASE  "Iniiial"  !  Initial  pass 

Xmin=0  !  Minimum  time.  Hours 

Xmax=5  !  Maximum  time.  Hours 

Xmajors.5 
Xminors.l 

CASE  "Follow  on"  !  Use  user  modified  values 

YminsYminmannal 

Xmax=Xmax_manual 
Xmajor=Xmajor_manual 
Xminor=Xmmor_manual 
END  SELECT 
END  SELECT 

SELECT  YscaleS 

CASE  "Auto  scale  Y"  !  Auto  scale  Y  axis 

CALL  Yscale(Ymin,Ymax,Yminor,Ymajor) 

CASE  "User  scale  Y"  !  User  scale  Y  axis 

SELECT  PassS 

CASE  "Initial"  !  Initial  pass 

Ytnin-User.fireq.min  !  Use  user  entered  min  fireq 
YmaxsUser.fireq.max  !  Use  user  entered  max  fireq 
YminsUser.mag.min  !  Use  user  entered  min  mag 
YmaxsUser.mag.max  !  Use  user  entered  max  mag 
Ymajor=(Y  max-Ymin)/5 
YminorsYmajor/5 
PassSs"Follow  on" 

Xscale$="Auto  scale  X" 

YscaleS="Auto  scale  Y" 

CASE  "Follow  on"  !  Use  user  modified  values 
YminsYmin.manuai 
YmaxsY  max.manual 
YmajorsY major.manual 
Y minors  Y minorjnanual 
END  SELECT 
END  SELECT 

WINDOW  Xmin,Xmax,Ymin,Ymax 
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6890! 

6900 

6910 

6920 

6930 

6940! 

6950 

6960 

6970 

6980 

6990 

7000 

7010! 

7020 

7030 

7040 


CALL  Lbl_axes(2,.  6  JPen4  XniinJCmax^Xmajor,  Ymin,  Y max, Y major) 
PEN  Pen2 

AXES  Xmirar.YminorXniin.YniinXmajor/Xiniiwr.Ymajor/Y minor ,3 
AXES  Xminor,YminorXni«,Ymax,Xmajor/Xininor,Ymajor/Yminor,3 

SELECT  Grid.typcS  .'Grid 

CASE  "Partial" 

GRID  Xmajor/YmajcrXmin/Ymin 
CASE  "Full" 

GRID  Xminor.YminOTXmin.YminJCmajcff/Ximnor.Ymajor/Yininor.l 
END  SELECT 

PENUP 

GRAPHICS  ON 
RETURN 


7060  Take_plot_data;  !  Collects  fireq  vs  temperature  data  during  the  run 

7080  ON  KEY  0  LABEL  "Take  Data"  GOTO  Take 

7090  FOR  Keynumber=l  TO  9 

7100  ON  KEY  Keynmnber  LABEL  ""  GOSUB  Commend 

7110  NEXT  Keynumber 

7120  Wait:  GOTO  Wait 

7130  Take:  ON  KEY  0  LABEL  ""  GOSUB  Commentl 

7140  PEN  Pen3 

7150  OUTPUT  KBD;"Taking  data 

7160  T0=TIMEDATE 

7170  SELECT  UPC$(Block3$) 

7180  CASE  "FLEXURAL" 

7190  FOR  1=0  TO  5*Argl 

7200  ENTER  720;Resfreq(I) 

7210  ENTER  722;Thennistor(I) 

7220  ENTER  724;Inphmag(I) 

7230  Time(I)=I/Argl 

7240  PLOT  Time(I)Jnphmag(I) 

7250  GOSUB  Computejedfreq 

7260  GOSUB  Youngjnod 

7270  OUTPUT  @Path_l  ;Enex(I),Thermistor(I), 

Resfreq(I),  Inphmag(I),Redfreq(I) 
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7280 

7290 

7300 

7310 

7320 

7330 

7340 

7350 

7360 

7370 

7380 

7390 

7400 

7410 

7420 

7430 

7440 

7450 

7460 

7470 

7480 

7490 

7500 

7510 

7520 

7530 

7540 

7550 

7560 

7570 

7580 

7590 

7600 

7610 

7620 

7630 


DISPI 

WAIT  118.315  (Adjust  for  1  sample  every  2  min 
NEXT  I 

CASE  TORSIONAL" 

FOR  J=0  TO  5*Argl 
ENTER  720-^esfreq(J) 

ENTER  722;Tbennistof(J) 

ENTER  724;Inphjnag(J) 

T*ime{J)=J/Argl 
PLOT  Time(J)Jnphmag(J) 

GOSUB  Compme_redfreq 
GOSUB- Shear jnod 

OUTPUT  @Path_l;Gtw^J),Thermistor(J)^esfreq(J), 
Inphmag(J),  Redfreq(J) 

DISPJ 

WATT  118.315 
NEXT  J 

CASE  "LONGITUDINAL" 

FOR  K=0  TO  5*Argl 
ENTER  720*JRes£rcq(K) 

ENTER  722;Thermistor(K) 

ENTER  724;Inpfamag(K) 

Time(K)=K/Argl 
PLOT  Time(K)4nphmag(K) 

GOSUB  Compute  jredfireq 
GOSUB  Lyoungjnod 

OUTPUT  @Path_l  ^ long(K),Thennistor(K) JResfreq(K), 
Inphmag(K),Redfreq(K) 

DISPK 
WAIT  118.315 
NEXTK 
END  SELECT 
T1=TIMEDATE 

DISP  "TT  TOOK  ";DROUND(T1-TO,4);"SECONDS" 

WAIT  10 

OUTPUT  KBD;Qear$;  !  Gears  the  CRT 

BEEP 

ON  KEY  0  LABEL  "Post  Process"  GOTO  7650 


7640  Waitl:  GOTO  Waitl 
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7650 


7670  Post_process:  !  Permits  user  to  extract  plots,  tables  &  otter  info 


7690  GCLEAR 

7700  CALL  Least_squares(5*Argl  ,Thennistor(*)tResfreq(*), 

A3,C MW?) 

7710  Slope=A  !  Keeps  call  to  Least_squares  to  one  program  line 

7720  Intercept=B 

7730  Correlation^ 

7740  Slope_error=D 

7750  Intrcpterr=E 

7760  Tmean=F 

7771  Fmean=G 

7780! 

7790  Actual_tmp_min=MIN(Therinistor(*))  !  Find  max  and  min 

7800  Actual jemp_max=NlAX(Thennistor(*))  !  temp  and  frequency 

7810  Actual  Jreqjnin=MIN(Resfreq(*)) 

7820  Actual  Jreq^max=MAX(Resfreq(*)) 

7830  Actual jnag_min=MIN(Inphmag(*)) 

7840  Actual_mag_max=MAX(Inphmag(*)) 

7850! 

7860  Volume=^5*(Length/100)*PinDiameter/100)A2 

7870  GOSUB  Percentjnoduli 

7880  GOSUB  Compute_modulus 

7890! 

7900  OUTPUT  KBD;Qear$;  !  Gear  the  CRT 

7910  OUTPUT  KBD;Hotne$;  !  Home  display 

7920  GCLEAR 

7930  PRINT  TABXY(l.l)  !  Start  at  top  with  blank  line 

7940  Center=(Screen-42)/2  !  Leading  spaces  for  centering 

7950  PRINT  TAB(Center);"Key  Purpose" 

7960  PRINT  TAB(Center);" - " 

7970  PRINT  TAB(Center);"  0  Plot  Frequency  vs  Time" 

7980  PRINT  TAB(Center);"  1  Plot  Temperature  vs  Time" 

7990  PRINT  TAB(Center);"  2  Plot  Frequency  vs  Temperature" 

8000  SELECT  UPC$(Block3$) 

8010  CASE  ="TORSIONAL" 

8020  PRINT  TAB(Center);"  3  Plot  Shear  modulus  vs  Temp" 

8030  CASE  ="FLEXURAL" 
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8040  PRINT  TAB(Center);"  3  Plot  Young's  modulus  vs  Temp" 

8050  CASE  ="LONGITUDIN  AL" 

8060  PRINT  TAB(Center);"  3  Plot  Young's  modulus  vs  Temp" 

8070  END  SELECT 

8080  PRINT  TAB(Center);"  4  Send  Information  to  the  Printer" 

8090  PRINT  TAB(Center);"  5  Select  the  graph  output  device" 

8100  PRINT  TAB(Center);"  6  Plot  In  phase  mag  vs  Temperature" 

8110  PRINT  TAB(Center);”  7  Plot  In  phase  mag  vs  Temperature" 

8120  PRINT  TAB(Center);"  7  Select  the  type  of  X  axis  scaling" 

8130  PRINT  TAB(Center);"  8  Select  the  type  of  Y  axis  scaling" 

8140  PRINT  TAB(Center);"  9  Exit  this  program" 

81501 

8160  ON  KEY  0  LABEL  Treq  vs  Tune"  GOSUB  Freqjime 
8170  ON  KEY  1  LABEL  "Temp  vs  Time"  GOSUB  Tempjime 

8180  ON  KEY  2  LABEL  Treq  vs  Temp"  GOSUB  Freqjemp 

8190  SELECT  UPC$(Block3S) 

8200  CASE  ="TORSIONAL" 

8210  ON  KEY  3  LABEL  "G  vs  Temp"  GOSUB  Shearjemp 

8220  CASE  =TLEXURAL" 

8230  ON  KEY  3  LABEL  T  vs  Temp"  GOSUB  Youngjemp 

8240  CASE  ="LONGITUDINAL" 

8250  ON  KEY  3  LABEL  T  vs  Temp"  GOSUB  Youngjemp 

8260  CASE  ELSE 

8270  ON  KEY  3  LABEL""  GOSUB  Commentl 

8280  END  SELECT 

8290  ON  KEY  4  LABEL  "Prim  Info"  GOSUB  Dumpjnfo 

8300  ON  KEY  5  LABEL  "Output  Device"  GOSUB  Output jievice 

8310!  ON  KEY  6  LABEL  "Grid  Option"  GOSUB  Grid.option 

8320  ON  KEY  6  LABEL  "Imag  vs  Temp"  GOSUB  Magjemp 
8330!  ON  KEY  7  LABEL  "X  scale  option"  GOSUB  Xscalejjption 

8340  ON  KEY  7  LABEL  "Imag  vs  Freq"  GOSUB  Mag_freq 

8350  ON  KEY  8  LABEL  "Y  scale  option"  GOSUB  Yscale.opdon 

8360  ON  KEY  9  LABEL  "Exit  Program"  GOSUB  Program_end 
8370! 

8380  Blink3:  WAIT  1 

8390  DISP  "Make  a  DECISION" 

8400  WAIT  1 

8410  DISP 

8420  GOTO  Blink3 
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8430  RETURN 


'..IlMiaaMiaiMiiiiiaMiitMiatMiaiaaifliaiMMiMMiMMiiMiiiiiMii 


8450  Grid_option:  !  Accepts  the  user's  choice  for  the  plot  grid 


8470  GRAPHICS  OFF  !  Turns  off  graphics  display 

8480  OUTPUT  KBD;Clear$;  !  Clean  the  CRT 

8490  DISP  Grid_type$;"  grid  is  currently  selected " 

8500  INPUT  "Enter  F  -  full;  P  -  partial;  N  -  No  grid 
?",ResponseS 

8510  SELECT  Responses 

8520  CASE  "F 
8530  Grid_type$=TulT 

8540  CASE  "P"  ' 

8550  Grid_type$="Partial" 

8560  CASE"N" 

8570  Grid_type$="No" 

8580  CASE  ELSE 

8590  DISP  "No  change" 

8600  WAIT  1 

8610  DISP 

8620  WAIT  1 

8630  END  SELECT 

8640  DISP  Grid.typeS;"  grid  is  selected" 

8650  WAIT  1 

8660  DISP 

8670  RETURN 


8690  Dumpjnfo:  !  Sends  selected  data  and  table  info  to  the  printer 


8710 

8720 

8730 

8740 

8750 

8760 

8770 

8780 

8790 

8800 

8810 


PRINTER  IS  701 

Perfskip$=CHRS(27)&CHR$(38)&CHRS(108)&CHR$(49)&CHR$(76) 

Fonnfeed$=CHR$(12) 

PRINT  PerfslripS  !  Skip  on  Perforation 

PRINT  Formfeeds 

PRINT  USING  "3/,#"  !  Three  line  feeds 

PRINT  "Filename  (Rod  ID):  "&Blockl$;TAB(70);"Page  1  of  2" 

PRINT  "Run:  "&Block2S 
PRINT  "Mode:  "&Block3$ 

PRINT  "Date:  "&Block4$ 

PRINT  USING  "3/,#"  !  Three  more  line  feeds 


8820! 

8830 

8840 

8850 

8860 

8870 

8880 

8890 

8900 

8910 

8920 

8930 

8940! 

8950 

8960 

8970 

8980 

8990 

9000 

9010 

9020 

9030 

9040 

9050 

9060 

9070 

9080 

9090 

9100 

9110 

9120 

9130 

9140 

9150 

9160 

9170 


Cl$="Time" 

C2$="Temperature" 

C3$="Frequency" 

C4a$="  Shear  modulus” 
C4b$=”Young's  modulus" 
C5$="In  Phase  Component" 
C6$="minutes" 

C7$="deg  C” 

C8$="Hz" 

C9$="Pa*10A" 

CIOS^DCV'  . 


SELECT  UPC$(Block35) 

CASE  "TORSIONAL” 

PRINT  USING  8980;Cl$,C2$,C3$,C4a$,C5$ 

IMAGE  4A.6X.1 1A,5X,9A,5X,13A,5X,18A 
PRINT 

PRINT  USING  9010;C6$,C7$,C8$,C9$vPower,C10$ 
IMAGE  5A,8X^A,11X^A,12X,6A2Z,14X^A,12X 
PRINT 

FOR  1=0  TO  5*Argl  STEP  Argl/2 

PRINT  USING  9050;Time(I).Thennistor(I)JResfreq(I), 
Gtor(I),Inphmag(I) 

IMAGE  X^D,8X,SDDDDD,4X,  DDDDDDDE^XJ).DDD,13X, 
SDDDDDDDE 
NEXT  I 

CASE  "FLEXURAL" 

PRINT  USING  9140;Cl$,C2$,C3$,C4bS,C5$ 

PRINT 

PRINT  USING  9010;C6$,C7$,C8S,C9$J>ower,C10$ 
PRINT 

FOR  J=0  TO  5*Argl  STEP  Argl/2 

PRINT  USING  9050;Time(J),Thermistor(J)Resfreq(J), 
Eflex(J),Inphmag(I) 

IMAGE  4A,6X,11A^X^A^X,15A,3X,18A 
NEXT  J 

CASE  "LONGITUDINAL" 

PRINT  USING  9140;Cl$,C2S,C3S,C4b$,C5$ 
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9180 

9190 

9200 

9210 

9220 

9230 

9240 

9250! 

9260 

9270 

9280 

9290 

9300 

9310 

9320! 

9330 

9340 

9350 

9351 

9360 

9361 
9370 

9380 

9390 

9400 

9410! 

9420 

9430 

9440 

9450 

9460 

9470 

9480 


PRINT 

PRINT  USING  9010;C6$,C7$,C8$,C9$JPowcr,C10$ 

PRINT 

FOR  K=0  TO  5*Argl  STEP  Argl/2 

PRINT  USING  9050*,rime(K),Thennistor(K),Resfreq(K), 
Elong(K),Inphmag(I) 

NEXTK 


END  SELECT 


PRINT  Formfeeds 
PRINT  USING  "3/,#" 

PRINT  "Filename  (Rod  ID): 
PRINT  "Run:  "&Block2S 
PRINT  "Mode:  ”&Block3$ 
PRINT  "Date:  “&Block4S 


!  Advance  to  top  of  next  page 
!  Three  more  line  feeds 
"&Blockl$;TAB(70);"Page  2  of  2’ 


PRINT  USING  "3/,#"  !  Three  more  line  feeds 

PRINT  "Physical  properties:” 

PRINT 

PRINT  USING  "3X,12A,14XJ)DDDDD";"Total  mass,  grams:", 
Totalmass 

PRINT  USING  "3X,12A,14XJ)DD.DDD";"Mass,  grams: ".Mass 
PRINT  USING  ”3X^0A,06X,DDD J)DD" ; "Length,  cm:",Length 
PRINT  USING  "3X^0 A,06XJ)DD.DDD” ; "Effective  Length,  cm:", 
Leff 

PRINT  USING  ”3X^2A,04XJ)DD.DDD"; "Diameter,  cm:"JDiameter 
PRINT  USING  "3X^lA,07XJ).DDDE”;"Volume,  cubic  meters:", 
Volume 

PRINT  USING  "3X,16A,09XJDDDD.D";"Density,  kg/mA3:"JDensity 


PRINT  USING  "3/,#"  !  Three  more  line  feeds 

PRINT  "Least-squares  fit  results  [freq  vs  temp]:" 

PRINT 

PRINT  USING  "3X,12A,13XJ)DDD.DDD";"Slope,  Hz/C:",SIope 
PRINT  USING  "3X,18A,07XJ)DDDJ)DD";"SIope  error,  Hz/C:", 
Slope.error 

PRINT  USING  "3X,14A,1 1  X,DDDDDDD";  "Intercept,  Hz:", 
Intercept 

PRINT  USING  "3X^0A,05X,DDDD.DDD";"Inteicept  error,  Hz:", 
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Intrcptcrr 

9490  PRINT  USING  "3X,12A,13XJ)DDD.DDDDD";"CorTeiaQon:", 

Correlation 

9500  PRINT  USING  "3X,20A,05X,DDDD.DDD";"Mean  temp,  C:",Tmean 

9510  PRINT  USING  "3X,19A,06X,DDDD.DDD";"Mean  frcq,  Hz:"  .Fmean 

9520  PRINT 

9530  PRINT 

9540  PRINT  "Other  statistics:" 

9550  PRINT 

9560  Tmin=Actual_temp_min 

9570  Tmax=Actual_temp_max 

9580  Tave=(Tmax+Tmin)/2 

9590  Fmin=Actuar_freq_mm 

9600  Fmax=Actual_freq_max 

9610  Fave=(Fmax+Fmin)/2 

9620  PRINT  USING  "3X,23A,02X,DDDD.DDD";"Min  temp,  C:",Tmin 

9630  PRINT  USING  "3X^3A,02XJ)DDD.DDD";"Average  temp,  C:",Tave 
9640  PRINT  USING  "3X^3A,02XJ)DDDJ)DD";"Maximum  temp,  C:",Tmax 
9650  PRINT 

9660  PRINT  USING  "3X,22A,03X,DDDD.DDD";"Min  frequency,  Hz:”,Fmin 

9670  PRINT  USING  "3X^2A,03X,DDDD.DDD";"Ave  frequency,  Hz:",Fave 

9680  PRINT  USING  "3X,22A,03X,DDDD.DDD";"Max  frequency,  Hz:",Fmax 

9690  PRINT 

9700  Ll$=  "Percent  change  in  shear  modulus  per  degree  C:  " 

9710  L2$="Percent  change  in  Young's  modulus  per  degree  C:" 

9720  L3$=Two  point  max-min  approach:" 

9730  L4$="Multi-point  least-squares  approach:" 

9740  SELECT  UPC$(BIock3$) 

9750  CASE  TORSIONAL" 

9760  Gmin=Scale _gtor_min 

9770  Gmax=Scale_gtor_max 

9780  PRINT  USING  "3X,30A,ZZ^X,D  J)DD";"Min  "&C4a$, 

Power,  Gmin 

9790  PRINT  USING  "3X,30A,ZZ,5X,D.DDD’,;,,Maximum  "&C4a$, 

Power,Gmax 

9800  PRINT  USING  "3/,#"  !  Three  more  line  feeds 

9810  PRINT  LIS 

9820  PRINT 

9830  PRINT  USING  "3X^7A,8X,DDD.DDD";L3$,Prcnt_per_c_2pt 

9840  PRINT  USING  "3X,35A,DDD.DDD";L4$,Prcnt_per.c_ls 
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9850  CASE  "FLEXURAL" 

9860  Emin=Scale_eflex_min 

9870  Emax=Scale_eflex_max 

9880  PRINT  USING  "3X^0A,ZZ,5XJ).DDD";  "Minimum  "&C4b$, 

Power, Emin 

9890  PRINT  USING  "3X,30A,ZZ,5XJ).DDD";"Maxiinum  "&C4b$, 

Power.Emax 

9900  PRINT  USING  "3/,#"  !  Three  more  line  feeds 

9910  PRINT  L2$ 

9920  PRINT 

9930  PRINT  USING  "3XJ7A,8X,DDDDDD";L3$,PrenLper_c_2pt 

9940  PRINT  USING  "3X,35A,DDDDDD";L4$,Prcntj>er_cJs 

9950  CASE  "LONGITUDINAL" 

9960  Emin=Scale_elong_inin 

9970  Emax=Scale_elong_max 

9980  PRINT  USING  "3XJ2A^Z^X,DDDD";  "Min  "&C4bS, 

Power  ,Emin 

9990  PRINT  USING  "3X^2A2Z,5XJ)JDDD";"Maximum  "&C4b$, 

Power.Emax 

10000  PRINT  USING  "W*  !  Three  more  line  feeds 

10010  PRINT  L2S 

10020  PRINT 

10030  PRINT  USING  "3X^7A,8XJ3DD.DDD"X3$J»rcnt_per_cJpt 

10040  PRINT  USING  "3X^5A,DDD.DDD";I^,Pnmt_per_cJs 

10050  END  SELECT 

10060  PRINT  Formfeeds  !  Advance  to  top  of  next  page 

10070  PRINTER  IS  CRT 

10080  RETURN 

10090!****************************************************************** 


10100  Freq_time:!  Produces  frequency  vs  time  graph  w/curve 


10120  GOSUB  Draw_freq_time 

10130  CALL  Generic_curve(Pen3,Time(*),Resffeq(*),5*Argl) 

10140  RETURN 


10160  Output_device:  !  Permits  user  to  route  graphs  to  screen  or  plotter 
101701****************************************************************** 


10180  OUTPUT  KBD;Clear$;  !  Gear  the  CRT 

10190  OUTPUT  KBD;Home$;  !  Home  display 
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!  Turn  off  the  graphics  display 


10200  GRAPHICS  OFF 

10210  SELECT  Plot_device$ 

10220  CASE  o"Plotter" 

10230  INPUT  "Graphs  appear  on  the  screen,  OK?  (Y/N) ", 

Response} 

10240  IF  UPCS(Response$)o"Y"  THEN 

10250  Plot_device$="Plotter" 

10260  DISP  "Graphs  will  be  sent  to  the  Plotter" 

10270  ELSE 

10280  DISP  "Graphs  will  remain  on  the  screen" 

10290  WAIT  1 

10300  DISP 

10310  GRAPHICS  ON 

10320  END  IF 

10330  CASE  "Plotter" 

10340  INPUT  "Graphs  are  sent  to  the  plotter,  OK?  (Y/N) ", 

Responses 

10350  IF  UPC$(Response$)o"Y"  THEN 

10360  Plot_deviceS=’’ Screen" 

10370  DISP  "Plots  will  be  sent  to  the  Screen” 

10380  ELSE 

10390  DISP  "Plots  will  stay  routed  to  the  plotter" 

10400  END  IF 

10410  END  SELECT 

10420  WATT  1 
10430  DISP 
10440  RETURN 

104501****************************************************************** 

10460  Xsca!e_option:  !  Permits  user  to  auto  or  manually  the  X  axis 
10470!  **^*^™****&****************+*********+****************************** 


10480 

OUTPUT  KBD.QearS; 

!  Gear  die  CRT 

10490 

OUTPUT  KBD;Home$; 

!  Home  display 

10500 

GRAPHICS  OFF 

!  Turn  off  the  graphics  display 

10510 

SELECT  XscaleS 

10520 

CASE  "Auto  scale  X” 

10530 

INPUT  "X  axis  is  automatically  scaled,  OK?  (Y/N) ", 

Responses 

10540 

IF  Response$o"Y"  THEN 

10550 

Xscale$="User  scale  X" 

10560 

INPUT  "Minimum  X  axis  value  ?",Xmin_manual 
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10570 

10580 

10590 

10600 

10610 

10620 

10630 

10640 

10650 

10660 

10670 

10680 

10690 

10700 

10710 

10720 

10730 

10740 

10750 

10760 

10770 

10780 

10790 

10800 

10810 

10820 

10830 

10840 

10850 

10860 

10870 

10880 

10890 

10900 


INPUT  "Maximum.  X  axis  value  ?",Xmaxjnanual 
INPUT  "Major  X  axis  increment  ?"tXmajor_manua] 
INPUT  "Minor  X  axis  increment  ?"tXminor_manual 
DISP  "The  X  axis  will  be  scaled  manually" 

ELSE 

DISP  "The  X  axis  will  remain  auto  scaled" 

WATT  1 
DISP 

GRAPHICS  ON 
END  IF 

CASE  "User  scale  X" 

INPUT  "X  axis  is  man  scaled,  OK?  (Y/N) ",  Responses 
IF  Response$o"Y"  THEN 
Xscale$="Auto  scale  X" 

DISP  "The  X  axis  will  be  scaled  automatically” 

ELSE 

INPUT  "Change  manual  limits  ?  (Y/N)  ".Responses 
IF  Response5="Y"  THEN 

INPUT  "Minimum  X  axis  value  ?",Xmin_manual 
INPUT  "Maximum  X  axis  value  ?",Xmaxjnanual 
INPUT  "Major  X  axis  increment  ?",Xmajorjnanual 
INPUT  "Minor  X  axis  increment  ?",Xminor_manual 
DISP  "X  axis  will  be  scaled  with  new  values” 

ELSE 

DISP  "The  X  axis  will  remain  manually  scaled” 
WATT  1 
DISP 

GRAPHICS  ON 
END  IF 
END  IF 
END  SELECT 
WATT  1 
DISP 
RETURN 


10920  Yscale_option:  !  Permits  user  to  auto  or  manually  scale  the  Y  axis 
10930!  *&************************************************************** 


10940  OUTPUT  KBD;Clear$;  !  Clear  the  CRT 

10950  OUTPUT  KBD,Home$;  !  Home  display 
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10960 

10970 

10980 

10990 

11000 

11010 

11020 

11030 

11040 

11050 

11060 

11070 

11080 

11090 

11100 

11110 

11120 

11130 

11140 

11150 

11160 

11170 

11180 

11190 

11200 

11210 

11220 

11230 

11240 

11250 

11260 

11270 

11280 

11290 

11300 

11310 

11320 

11330 

11340 

11350 


GRAPHICS  OFF  !  Turn  off  the  graphics  display 

SELECT  YscaleS 
CASE  "Auto  scale  Y" 

INPUT  "Y  axis  is  auto  scaled,  OK?  (Y/N) ",  Response! 

IF  Response$o"Y"  THEN 
Yscale$="User  scale  Y* 

INPUT  "Minimum  Y  axis  value  ?",Ymin_manual 
INPUT  "Maximum  Y  axis  value  ?",Ymax_manuai 
INPUT  "Major  Y  axis  increment  ?",Ymajor_manual 
INPUT  "Minor  Y  axis  increment  ?",Yminor_manual 
DISP  "The  Y  axis  will  be  scaled  manually" 

ELSE 

DISP  "The  Y  axis  will  remain  auto  scaled" 

WATT  1 
DISP 

GRAPHICS  ON 
END  IF 

CASE  "User  scaler 

INPUT  "Y  axis  is  man  scaled,  OK?  (Y/N)  ",  Responses 
IF  ResponseSoT  THEN 
Yscaie$="Auto  scale  Y* 

DISP  "The  Y  axis  will  be  scaled  automatically” 

ELSE 

INPUT  "Change  manual  limits  ?  (Y/N)  "  .Response! 

IF  Response$="r  THEN 

INPUT  "Minimum  Y  axis  value  ?”,Yminjmanual 
INPUT  "Maximum  Y  axis  value  ?”,Ymaxjnanual 
INPUT  "Major  Y  axis  increment  ?",Ymajor_manual 
INPUT  "Minor  Y  axis  increment  ?",Yminorjnanual 
DISP  "Y  axis  will  be  scaled  with  new  values" 

ELSE 

DISP  "The  Y  axis  will  remain  manually  scaled” 
WAIT  1 
DISP 

GRAPHICS  ON 
END  IF 
END  IF 
END  SELECT 
WAIT  1 
DISP 
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11360  RETURN 


11380  Temp.time:  !  Produces  Temperature  vs  Time  graph  w/curve 


11400  Main_title$=Label$(7) 

11410  Sub_dtle$=Label$(8) 

11420  X_axis_name$=LabelS(3) 

11430  Yjtxis_name$=Labei$(5) 

11440  Xmin=0  !  Minimum  dme  is  0  hours 

11450  Xmax=5  !  Maximum  time  is  5  hours 

11460  Y min=Actual_temp_min  !  Actual  minimum  temperature 

11470  Y max=Actual_temp_max  !  Actual  maximum  temperature 

11480  GOSUB  Generic_plot 

11490  CALL  Generic_curve(Pea3fTime(*),Thennistor(*)t5*Argl) 

11500  RETURN 


11520  Freq_temp:  !  Produces  Frequency  versus  temperature  graph  w/curve 


11540 

11550 

11560 

11570 

11580 

11590 

11600 

11610 

11620 

11630 

11640! 

11650! 

11660! 

11670! 

11680 


Main_title$=Label$(6) 

Sub_title$="" 

X_axis_name$=Label$(5) 

Y_axis_name$=Labei$(4) 

Xmin=Actual_temp_mm  !  Actual  minimum  temperature 

Xmax=Acmal_temp_max  !  Actual  maximum  temperature 

Y min=Actual_freq_min  !  Actual  minimum  frequency 

Ymax-Actual Jreqjnax  !  Actual  maximum  frequency 

GOSUB  Generic_plot 

CALL  Generic_curve(Pen3,Thermistor(*)JResfreq(*)^*Argl) 
PEN  Pen4 

MOVE  Actual_temp_min^Actual_temp_min*Slope+Intercept 
DRAW  Actual_temp_max,Actual_temp_max*Slope+Intercept 
PEN  0 
RETURN 


11700  Create_new_file:  !  Creates  a  file  for  a  new  data  run 


11720  Pass$="Initiar 

1 1 730  Grid_type$=Default_grid$ 
11740  XscaJe$="User  scale  X" 


11750  Yscale$="User  scale  Y" 

11760  PloLdeviceS*"  Screen" 

11770  BEEP  1000^1 

11780  ON  ERROR  GOTO  Fix2 

11790  INPUT  "Filename  to  store  the  data  ?",Filename$ 

11800  CREATE  BDAT  FilenameS^O 

11810  ASSIGN  @Path.l  TO  Filenames 

11820  OFF  ERROR 

11830  OUTPUT  @Path_l;Blockl$31ock2$31ock3$31ock4S 

11840  OUTPUT  @PaftJ;MassXengthJ)iameter,I)ensity 

11850  BEEP  1000..1 

11860  DISP  "Output  will  be  stored  under  Filename:  "filenames 

11870  WAIT  1 

11880  DISP 

11890  IF  Flag=l  THEN 

11900  FOR  1=0  TO  5*Argl 

11910  OUTPUT  @Path.l^Iod(I),Thennistora)3esfieq(I), 

lnphmag(I),Redfreq(I) 

11920  NEXT  I 

11930  END  IF 

11940  RETURN 


11960  Open_oid_file:  !  Opens  and  retrieves  data  from  an  existing  file 


11980  PassS="Follow  on" 

11990  Grid_typeS=Default _gridS 

12000  XscaleS="Auto  scale  X" 

12010  Yscale$="Auto  scale  Y" 

12020  Plot_device$="  Screen" 

12030  BEEP  1000,. 1 
12040  ON  ERROR  GOTO  Fixl 

12050  INPUT  "Filename  to  retrieve  the  data  ?"3ilename$ 

12060  ASSIGN  @Path_l  TO  Filenames 

12070  ENTER  @Path.l;Blockl$31ock2S31ock3$31ock4$ 

12080  ENTER  @Path_l  ^dassXengtiuDiameter density 

12090  BEEP  1000,1 

12100  DISP  "Retrieving  data  stored  under  Filename:  "filenames 

12110  FOR  1=0  TO  5*Argl 

12120  ENTER  @Path_l  ;Mod(I),Tbennistor(I)3esfteq(I), 

Inphmag(I),Redfreq(I) 
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12130  Time{D=I/Argl 

12140  NEXT  I 

12130  OFF  ERROR 
12160  ASSIGN  @Path_l  TO  * 

12170  RETURN 

12190  Fixl:  ! 


WyJMi 


12210 

12220 

12230 

12240 

12250 

12260 

12270 

12280 

12290 

12300 

12310 

12320 

12330 

12340 

12350 

12360 


VXJii 


12380  Fw2:  ! 


**  *  iTEIIZIIIIZIIIIII 


12400  SELECT  ERRN 

12410  CASE  53 

12420  DISP  "Limit  file  names  to  10  characters.  No  punctuation” 

12430  CASE  54 

12440  DISP  "Duplicate  filename!  Try  another  name." 

12450  END  SELECT 

12460  BEEP  300,.  1 

12470  WATT  1 

12480  DISP 

12490  WATT  1 

12500  GOTO  Create_new_file 

12510!****************************************************************** 


SELECT  ERRN 
CASE  53 

DISP  "Limit  file  mum  to  10  characters.  No  punctuation" 
CASE  56 

DISP  "Tins  file  doesn't  exist  on  the  data  disk” 

CASE  58 

DISP  "This  file  is  not  a  BDAT  file" 

CASE  ELSE 

DISP  "  This  file  can  not  be  processed  by  RT20B” 

END  SELECT 
PRINTER  IS  CRT 
BEEP  3004 
WATT  1 
DISP 
WATT  1 

GOTO  Open_old_file 
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12520  Qosejile:  IQoses  new  data  file  after  data  collection  is  over 


rix't'i 


12540  GCLEAR 

12550  ASSIGN  @Path_l  TO  * 

12560  BEEP  1000,1 

12570  DISP  "Data  is  stored  under  Filename:  "filename! 

12580  WATT  1 

12590  DISP 

12600  RETURN 


12620  Program.end:  !  Shuts  down  shop,  plays  a  little  melody 


tom 


12640  OUTPUT  KBD;Clear$; 

12650  GRAPHICS  OFF 

12660  CONTROL  1,12,*0 

12670  BEEP  157,1 

12680  BEEP  201,1 

12690  BEEP  178,1 

12700  BEEP  272,1 

12710  WATT  J 

12720  BEEP  272,1 

12730  BEEP  178,1 

12740  BEEP  157,1 

12750  BEEP  201,1 

12760  DISP 

12770  DISP  "Press  RUN  when  you  are  ready  for  another  try 

12780  END 


12800  SUB  Label(Csize,Asp.radoJLdirXorg^en^Cpos,Ypos,Text$) 


12820!  This  sub  defines  systems  variables  (Csize,  LDIR,  etc.), 
12830!  and  labels  the  text  (if  any)  accordingly. 

12840  DEG 

12850  CSIZE  Csize^sp.ratio 

12860  LDIR  Ldir 

12870  LORG  Lorg 

12880  PEN  Pen 

12890  MOVE  Xpos.Ypos 

12900  IF  TextSo”  THEN  LABEL  USING  "#4C";Text$ 


12910  PENUP 
12920  SUBEND 
129301" 

12940  SUB  Lbl_axes(Csize,Asp_rauo^enJ(min,Xmax,Xstep,Yaiin,Ymax,Ystep) 


12960  DEG 

12970  CSIZE  Csizc,Asp_ratio 

12980  PEN  Pen 

12990  CUP  OFF 

13000  LDIRO 

13010  LORG  6 

13020  Yrangc=Ymax-Ymm 

13030  Yoffset=.02*Yrange 

13040  FOR  L=Xmin  TO  Xmax  STEP  Xstep 

13050  MOVE  L,Ymia-Yoffset 

13060  IF  ABS(LK001  THEN  L=0 

13070  LABEL  USING 

13080  NEXT  L 

13090  LORG  8 

13100  Xoffset=.02*(Xmax-Xmin) 

13110  IF  Yrange<=l  THEN 

13120  Mmax=DROUND(Y range/Y step.l) 

13130  Yval=Ymin 

13140  FOR  M=0  TO  Mmax 

13150  IF  ABS(YvaI)<=.00I  THEN  YvaNO 

13160  MOVE  Xmin-Xoffset,Yval 

13170  LABEL  USING  "#,SDJ)D";Yval 

13180  Yval=Yval+Ystcp 

13190  NEXT  M 

13200  ELSE 

13210  FOR  M=Ymin  TO  Ymax  STEP  Ystep 

13220  IF  ABS(M)<=.001  THEN  M=0 

13230  MOVE  Xmin-XoffseU4 

13240  LABEL  USING 

13250  NEXT  M 

13260  END  IF 

13270  CUP  ON 

13280  PENUP 

13290  SUBEND 

13300!******************************"'*********************************** 
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13310  SUB  BlockJnfo(Blockl$,Block2$,Block3$,Block4$,Xgdumax,Pen,Pen2) 


13330  CALL  Ubt\QMa$m2XmKD&:  "ABlocklS) 

13340  CALL  Labcl(3,.6,0,5J>ciuXgduniax/3,2,'RuD:  "&Block2S) 

13350  CALL  Label(3,A0,5.Pen.Xgduinax^,2,"Mode:  "&Block3$) 

13360  CALL  Label(3,.6,0,8j>en,.97*Xgdumax,2,Block4$) 

13370  MOVE  0,4 

13380  PEN  Pea2 

13390  DRAW  133,4 

13400  PENUP 

13410  SUBEND 


13430  SUB  Lcast_squares(Imax^((*),Y(*),SlpJnl,Cor,Slp_erJnt_er, 
Xmcan.Ymcan) 


13450  DISP  "Computing  least-squares  fit _ " 

13460  Sumx=0 

13470  Sumy=0 

13480  Sumxx=0 

13490  Sumxy=0 

13500  FOR  1=0  TO  Imax 
13510  Sumx=Sumx+X(I) 

13520  Sumxx=Sumxx+X(I)A2 

13530  Sumy=Sumy+Y(I) 

13540  Sumxy=Sumxy+X(I)*Y(I) 

13550  NEXT  I 

13560  Delta=(Imax+l)*Sumxx-SumxA2 

13570  Int=(Sumxx*Sumy-Sumx*Sumxy)/Delta 

13580  Slp=((Imax+l  )*Sumxy-Sumx*Sumy)/DeIta 

13590  Sumerrerr=0 

13600  FOR  J=0  TO  Imax 

13610  Sumerrerr=Sumerrerr+(Y  (J>Int-Slp*X(J))A2 

13620  NEXT  J 

13630  Sigmayy=Sumerrerr/(Imax+l-2) 

13640  Sigmay=SQR(Sigmayy) 

13650  Int_er=Sigmay*SQR(Sumxx/Delta) 

13660  Slp_er=Sigmay*SQR((Imax+l)/Delta) 

13670  Xbarsum=0 

13680  Ybarsum=0 
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13690  FOR  K=0  TO  Imax 
13700  Xbarsum=Xbar$um+X(K) 

13710  Ybarsum=Ybarsum+Y(K) 

13720  NEXT  K 

13730  Xmean=Xbarsum/(Imax+ 1 ) 

13740  Ymean=Ybarsum/(Imax+l ) 

13750  Sigmaxxsum=0 

13760  Sigmayysum=0 

13770  Sigmaxysum=0 

13780  FOR  L=0  TO  Imax 

13790  Sigmaxxsum=Sigmaxxsum-KX(L)*Xmean)A2 

13800  Sigmayysum=Sigmayysum+(Y  (L)-Ymean)A2 

13810  Sigmaxysnm=Sigmaxysum+(X(L)-Xmcan)*(Y(L)-Ymcan) 

13820  NEXT  L 

13830  Sigmax=SQR(Sigmaxxsuin/(Imax+l)) 

13840  SigmaylsSQR(Sigmayy$um/(lmax+l)) 

13850  Sigmaxy=Sigmaxysum/(Imax+l) 

13860  Cor=Sigxnaxy/(Sigmax*Sigmayl) 

13870  DISP 
13880  SUBEND 

13900  SUB  ScaIe(L,R,B,T,Top,Left,Xcenter,Ycenter,Xgdumax,Ygduniax) 


13920  Top=T*Ygdumax 

13930  Bottom=B  *  Y gdumax 

13940  Left=L*Xgdumax 

13950  Righr=R*Xgdumax 

13960  Xcenter=(Ri  ght+Lef t)/2 

13970  Ycenter=(Top+Bottom)/2 
13980  VIEWPORT  Left JUght3ottom, Top 

13990  SUBEND 

14000!**i^^<^^**aW^'<******************************<t*******111******** 

14010  SUB  Yscale(Ymin,Ymax,Yminor,Ymajor) 

14020!**********  ******************************************************  *♦ 

14030  DIM  Diff(36),Mmor(36) 

14040  DATA  .lt.005^,.01,^5,.01,.3,.02,.4,.02,J,.02,.75,.05,l,.05 

14050  DATA  5,.2,10,.5,15,1,20,U5,1,30^,40,2,50, 2, 75^,100,5,125,5 

14060  DATA  150,10,200,10,250,10,300,20,400,20,500,20,750,50,1000,50 
14070  DATA  1250,50,1500,100,2000,100,2500,100,3000,200,4000,200,5000 
14080  DATA  200,7500,500,10000,500 
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14090  FOR  1=1  TO  36 
14100  READ  Diff(D.Mmor(D 
14110  NEXT  I 
14120  Ynmge=Ymax-Ymm 

14130  Index=l 

14140  WHILE  Yrange>Diff(Index) 

14150  Index=Index+l 

14160  END  WHILE 
14170  Yminor=Minor(Index) 

14180  Ymajor=Diff (Index  )/5 

14190  IF  Ymin<0  THEN 

14200  Newmin=Ymin-Ymajor+ABS(Ymin  MOD  Ymajor) 

14210  ELSE 

14220  Newmin=Ymin-(Y min  MOD  Ymajor) 

14230  END  IF 

14240  Newmax=Newmin+Diff (Index) 

14250  WHILE  Ymax>Newmax 

14260  Index=Index+l 

14270  Yminor=Minor(Index) 

14280  Y major=Diff (Index)/5 

14290  IF  Ymin<0  THEN 

14300  Newmm=Ymin-Ymajor+  ABS(Ymin  MOD  Ymajor) 

14310  ELSE 

14320  Newmin=Ymin-{Y min  MOD  Ymajor) 

14330  END  IF 

14340  Newmax=Newmin+Diff (Index) 

14350  END  WHILE 
14360  Ymax=Newmax 

14370  Ymin=Newmin 

14380  SUBEND 


14400  SUB  Xscale(Xmin,Xmax,Xminor,Xmajor) 


DIM  Diffx(40),Minorx(40),Majorx(40) 

DATA  5.,0^,l.,6.,.2,l.,7,^,1.0,8.0,.2,l,lO,.5Z0,12,J 
DATA  2., 14.,3 ,2.0, 16.,.5,2.,20., 1.0, 4.0,25.0, 1.03*030.0,1.0 
DATA  5.035.0,1.03.0, 40.,1.,5.30.,10,10.,60.3.0,10.,70., 2. 
DATA  10.0,80.0, 2.0, 10.0, 100.03.0, 20.0,120.,5.,20.,140.3.,20. 
DATA  200.0,10.0,20.0, 400.0,10.30., 600..20., 100., 800.,50.,100. 


14420 

14430 

14440 

14450 

14460 

14470 
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14480  DATA  1000.0,50.0,100.0, 1200.0,50.0,200.0, 1400.0,50.,200.,1 600. 

14490  DATA  50.0^00.,1800.,100j00.0000., 100.^.^500., 100.^00. 

14500  DATA  3000.,100.^00.0,3500.,100.,500.0,4000.,200., 400.0 

14510  FOR  Ial  TO  34 
14520  READ  Diffx(U,Mmoa(I),Majon(I) 

14530  NEXT  I 

14540  Xrange=Xmax-Xmin 

14550  Index*  1 

14560  WHILE  Xrange>Diffx(Index) 

14570  bdexslndex+l 

14580  END  WHILE 

14590  Xminar=Mmorx(Index) 

14600  Xmajor=Majorx(Index) 

14610  IF  XmiiKO  THEN 

14620  Ncwmin=Xmiii-Xmajor+ABS(Xmin  MOD  Xmajor) 

14630  ELSE 

14640  Ncwmin=Xmin-(Xmin  MOD  Xmajor) 

14650  END  IF 

14660  Newmax=Newmin+Diffx(Index) 

14670  WHILE  Xmax>Newmax 

14680  Index-Index+1 

14690  Xminor=Mmorx(Index) 

14700  Xmajor=Majorx(Index) 

14710  IF  Xmm<0  THEN 

14720  Newmin=Xmin-Xmajor+ABS(Xmin  MOD  Xmajor) 

14730  ELSE 

14740  Newmin=Xmin-<Xniin  MOD  Xmajor) 

14750  END  IF 

14760  Ncwmax=Ncwmin+Diffx(Indcx) 

14770  END  WHILE 
14780  XmaxsNewmax 

14790  Xmin=Ncwmin 

14800  SUBEND 

14820  SUB  Generic_curve(Pen,X(*),Y(*),Max_pomt) 


14840  PEN  Pea 

14850  FOR  Point=0  TO  Max_point 

14860  PLOT  X(Pomt),Y(Point) 
14870  NEXT  Point 
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14970  FOR  1=0  TO  Inux 

14980  Tcmp(I)=Thermistor(I)+273.15 

14990  X(I)=-129*(Temp(D-283.15)/(107+Temp(I)*283.15) 

15000  Y(I>=X(D*LOG(10) 

15010  Alphat(I)=EXP(Y(D) 

15020  Redfrcq(I)=Resfrcq(l)*Alphai(I) 

15030  NEXT  I 
15040  SUBEND 
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APPENDIX  X.  PR- 15 *2  DATA 


CONTENTS: 

-  PR1592  FIRST  TORSIONAL  MODE  DATA  (TABLE  E.l) 

-  PR1592  SECOND  TORSIONAL  MODE  DATA  (TABLE  E.2) 

-  PR-1592  DIRECT  Q  MEASUREMENT  DATA  AT  VARIOUS  TEMPERATURES 
FOR  CALIBRATION  OF  THE  IN-PHASE  VOLTAGE  (FIGURES  E.1-E.12) 
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TABLE  B.l  PR1592  FIRST  TORSIONAL  NODE  DATA 


SHEAR  (CORRE 


I  C 


RES  FREQ 


3.558E-3 

-14  571 

73355994 

3.440E-8 

-13.333 

721.35461 

3.340E-8 

-12.741 

710.72646 

3.249E-8 


3.175E-8 


3.1  15E-8 


3.047E-8 


2.97 1  E*8 


2.884E-8 


2.802E-8 


2.718E-8 


2.634E-8 


2.554E-8 


2.479E-8! 


2.399E-3 


2.30IE-8 


2  229E-8 


2.  '•  63E-8 


2  1 12E-8 


2.01  3E-8I 


1  944E-3 


1.875E-3 


1  803E-8 


1.734E-3 


1  663E-8 


1.614E-8 


1.553E-8 


1.499E-8 


1.445E-8 


1.395E-3 


1.350E-8 


1.305E-3I 


1.264E-8 


1.216E-8 


1.172E-8 


1.I30E-8 


1.039E-8 


1.049E-3I 


1.01  4E*3 


9.799E-7 


9.479E-7 


9.156E-7 


8.362E-7I 


S.583E-7! 


3.314E-7I 


3  057E-7I 


-12.262 


-1  1.877 


-1  1.553 


-1 1.039 


-10.518 


-9.887 


-9.335 


-8.761 


-3.145 


-7.596 


-7.0951 


700.97543 


692.981 1 


686.38914 


678.93261 


670.38667 


660.52178 


651.01918 


641.19773 


imkm 


INPHASE 


-0.204615 


felMI 


-0  203729 


■BSE3EE 

ill  il 


-5.901 


-5.179 


-4.545 


-3.963 


-3  243 


-2.588 


-1.99 


-1.2321 


-0.618 


612.32536 


60237463 


539.97692 


580.62669 


572.65712 


565.16869 


551.78365 


5*2.2671  1 


532.578481 


522.23082 


-0.1956S5 


-0.200562 


-0  209166 


-0.21  1506 


-0.207602 


-0.204488 


-0.206024 


-0  20645 


EK-SEM 


0.627 


1.423 


2.082 


2.734 


4.0321 


4.69 


5.3561 


6.011 


6.673 


7.339 


8.002 


8.697 


9.384| 


10.024 


10.702 


1  1.364 


12.028 


12.692 


13.3331 


14.012! 


512.10938 

-0.209376i 

502.29582 

-0.213031 

494.09934 

-0.216686! 

485.46042 

-0.218909 

476.17153 


467.56027 


459.405 


451.83432 


437.19854 


428.83433 


420.96193 


413.45877 


405.9057 


398.38314 


391.58948 


334.983 


78.6403 


372.14826 


366.12717 


360.31837 


354  6074: 


349.0953 


-0.220244 


-0.222705 


-0.229338 


-0.2261  1 


-0.227894 


-0.230282 


-0.234065 


-0.23729 


-0.242718 


3B2E32 


-0.252229 


-0.2584321 


-0.262851 


-0.2731391 


0.2786351 


-0.2S3906 


TABLE  B.l  CONTINUED 


A 

B 

c 

0 

7.8 1  3E-7 

14.669 

-0.290056 

■EB23E 

15.354 

336  7007 

-0.295068 

7.355E-7 

16.01 

333.53144 

7.138E-7 

16.678 

328.573 

-0.30658 

6.943E-7 

17.339 

324.05237 

-0.31346 

6  761E-7 

18 

319.78516 

Ml  IT  T 

6.585E-7 

18.672 

31559445 

MffiRPfF! 

6.425E-7 

19.328 

31 1.73893 

—a  it  n 

6.256E-7 

20.003 

307  50354 

-0.32995 

6.049E-7 

20.651 

302.47952 

-0J3045 

5.870E-7 

21.312 

297  97926 

■EFIM! 

S.721E-7 

21.99 

294.1769 

TIlTTll 

5.586E-7 

22.652 

290.66802 

-0.360781 

5.457E-7 

23  301 

287.30684 

Twnri 

5.334E-7 

23.959 

284.04529 

-0.37485 

5.220E-7 

24.623 

230.9944 

-0.38154 

5. 11 IE-7 

25.259 

278.03196 

0  38702 

5.007E-7 

25.925 

275.19593 

-0.3925 

4.8896-7 

26.595 

271  94384 

-0.38424 

4.75IE-7 

27.248 

268.07988 

4.622E-  7 

27.908 

264.40884 

-0.4035 

4.51SE-7 

28.58 

261.32993 

-0.41436 

4.420E-7 

29.225 

253.57502 

■■EXEQE 

4.332E-7 

29.873 

255.98542 

■1  Mil  1 1 

4.24IE-7 

30.532 

253.27751 

-0.43964 

4.I38E-7 

31.179 

250. 1 834 

-0.45577 

4.065E-7 

31.882 

247  95498 

-0.46195 

3.988E-7 

32.534 

245.60855 

-0.46908 

3.917E-7 

33.178 

243.41338 

-0  47417 

3.856E-7 

33.355 

241.51547 

■HXEEEB 

3.790E-7 

34.508 

239.41813 

■■CEEHE 

3.720E-7 

35.193 

237.2095 

■KXEEH3 

3.649E-7 

35.853 

234.93002 

-0.49373 

3.571E-7 

36  521 

232.42163 

-0.504741 

3.508E-7 

37.188 

230.33902 

-0  51896 

3.450E-7 

37.866 

228.43861 

-053036 

3.391E-7 

33534 

226.48174 

-0.54533 

3.336E-7 

39.194 

MKES331; 

3.28SE-7 

39.862 

223.00513 

-0.56788j 

3.236E-7 

40.535 

221.24825 

3.194E-7 

41.218 

219.79415 

-0.59178 

3.151E-7 

41.876 

213.32857 

-0.60326 

3.1 1 0E*7 

42.535 

216.8855 

-0.61481 

3.070E-7 

43.236 

215.47867 

-0.62977 

3  034E-7 

43.899!  214.21148 

-0.63979 

2.996E-7 

44  562 

212.9569 

■KXSEESj 

2.965E-7 

45.237 

21  '  77474 

-0.66305i 

TABLI  S.l  CONTINUED 


■El _ yss! _ 

A 

8 

C 

D  1 

94 

2.933E-7 

45.907 

210.64333 

95 

—  III  1  III 

46.589 

209.5051  1 

-0.679921 

96 

2.S69E-7 

47.264 

208.33231 

■■Enrrmn 

97 

2.840E-7 

47.958 

207.27333 

— ESEEHI 

96 

2.810E-7 

48.628 

206.151 1  1 

MHFffl 

99 

2.783E-7 

49.284 

205.16444 

— S3HEE 

100 

2.756E-7 

49.976 

204.16496 

■— SFFTR 

101 

2.729E*7 

50.637 

203.16787 

wm*zm\ 

102 

2.705E*7 

51.31 

202.27369 

103 

2.682E-7 

51.964 

—3 BiflB 

— x*mn 

104 

2.658E-7 

52.645 

200.50693 

—run  m 

105 

2.635E-7 

53.271 

199.64121 

-0.78971 

106 

2.613E-7 

53.937 

198.79755 

-0.80084 

107 

2.591E-7 

54.617 

197.98237 

— H0EH3 

108 

2.571E-7 

55.277 

197.1887 

-0.82831 

109 

2.549E-7 

55.947 

196.35229 

IMEXQiEE 

1  1C 

2.529E-7 

56.6 1 9 

■■EfTKTXFE? 

1 1 1 

■ms 

57.287 

194.85541 

niiia^w  i'll  ii  i 

112 

2.492E-7 

57.978 

■HlH-fcgi 

■KTTVi.Vi.1 

113 

2.474E-7 

58.626 

—HZ ! 

1 14 

2.457E-7 

59.296 

IIHHIHI 

■KEEEDD 

115 

2.439E-7 

59.956 

192.08169 

—r mm 

1 16 

2.424E-7 

60.622 

191.47189 

— EXEZ333 

117 

2.407E-7 

61.291 

190.80439 

-0.950091 

118 

2.392E-7 

61.973 

190.21615 

119 

— T1TTT1T1 

62.65 

189.59101 

— ra 

120 

2.361E-7 

63.308 

188.97801 

-0.99487 

121 

2.346E-7 

64.458 

188.38921 

-1.01  129 

122 

2.333E-7 

65.  J 15 

187.85093 

-1.02519 

123 

— TW'lTn 

65.758 

187  26996 

-1.04168 

124 

2.304E-7 

66.417 

186.69455 

-1.75918 

125 

2.293E-7 

67.073 

186.21917 

-1.07342 

126 

2.280E-7 

67.736 

185.63471 

-1.09103 

127 

2.268E-7 

68.377 

185.21015 

-l.lt  157 

128 

2.253E-7 

69.031 

Ml  1  % T 1  U  ' 

-1.12665 

129 

2.247E-7 

69.694 

184.36136 

-1.I4J4 

130 

2.236E-7 

70.345 

183.88961 

-1  15939 

131 

2.225E-7 

71.007 

183.44837 

-1.1 749 

132 

2.215E-7 

71.671 

183.04204 

-1.1895 

133 

2.206E-7 

72.344 

18254761 

-1.19659 

134 

2.195E-7 

73.028 

182.2121  1 

-1.2236 

135 

2.193E-7 

73.683 

182. 11  702 

-1.24404 

136 

2.130E-7 

78.393 

179  47319 

-0.64647 

137 

2.121E-7 

79.063 

179. 11  135 

-0.65615 

138 

2.1  15E-7 

79.704 

175.54323 

-0  66269 

139 

2  107E-7 

50.379 

175.51625 

-0.67139 

140 

2  100E-7 

31.048 

1  ”5.2 1 67 1 

-0.66125 

150 


TABLE  B.l  CONTINUED 


Label 

A 

B 

C 

_ D _ 

141 

2.093E*7 

81.707 

177.93266 

9MM2X3lXEI3 

142 

2.0866*7 

82.393 

177  63697 

■DXZiniE 

143 

2.0796*7 

83.048 

177.34139 

■gram 

144 

2.0746*7 

83.718 

177.09982 

145 

2.0776*7 

82.41 1 

177.22667 

-0.718571 

TABLE  1.2  PR1592  SECOND  TORSIONAL  MODE  DATA 


ms 

late! 

A 

 B 

C 

0  1 

Ljtet 

corrtctw 

SHEAR  MOOUIU 

RESONANT  FRE< 

IBM**.' 

1 

stear  modulus 

3.581E-8 

-13  765 

1479  4419 

_ 2. 

2nd  torsional 

3  466E-8 

-13.05 

1455.3742 

-0  2035691 

3 

mod*  3-16-93 

3.369E-8 

-12.491 

1434.9878 

4 

3.278E-8 

-1  1  962 

1415.4678 

5 

3.181  E*8 

-1 1.335 

1394.265 

5 

3.093E-3 

-10.843 

1374.8478 

■BEEH-S3 

7 

2.990E-8 

-10.132 

1351.8373 

-0.19346 

9 

2.903E-8 

-9.649 

1331.9333 

-0.192132 

9 

2.805E-8 

-8.967 

1309.3382 

10 

2.724E-8 

-8.453 

1290.231 16 

■Pirrarfl 

1  1 

2.631E-8 

-7.74 

1268.17618 

12 

2.552E-8 

-7.265 

1248.96705 

-0.1875931 

13 

2.465E-8 

-6.59 

1227.32184 

14 

2.376E-8 

-5.896 

1205.01801 

■HMH 

15 

2.294E-8 

-5.238 

1  184.12597 

^wnnr 

16 

2.215E-8 

-4.61  1 

1  163.43138 

■Bmagi 

17 

2.140E-8 

-4.012 

1 143.55856 

-0.186475 

IS 

2.067E-3 

-3.394 

1 123.99536 

-0.186919 

19 

1.999E-8 

-2.802 

1  105  33972 

-0.186537 

20 

1.936E-8 

-2.243 

1087.7452 

-0. 1 86059 

21 

1.372E-8 

-1.602 

1069.73184 

22 

1.812E-8 

-1.076 

1052.48239 

MEMFFFm 

23 

1.75IE-8 

-0.468 

mpmmi 

24 

1  688E-8 

0.179 

1015.83959 

-0.187108 

25 

1.626E-8 

0.817 

996.8548 

-0.187935 

26 

— W  1  ill  I 

1.494 

978.62057 

-0.189358 

27 

1.51  IE-8 

2.M6 

961.12765 

-0. 1 90594 

28 

1.457E-8 

2.745 

■KHCH3E 

29 

1.405E-8 

3  394 

926.73233 

30 

1.354E-3 

4.044 

909.76509 

-0.193422 

31 

1.302E-8 

4.692 

892.05993 

-0. 1 97496 

32 

1.260E-8 

5.342 

877.50896 

-0.2001 18 

33 

1.21 6E-8 

6.002 

361.99127 

-0.201633 

34 

1  172E-8 

6.642 

846.46443 

MEEZBEE 

35 

1.132E-8 

7.301 

831.82277 

-0.206897 

36 

1.093E-8 

7.95 

817.17401 

-0.209815 

37 

1.054E-8 

8.664 

802.76867 

■J’HH'IH 

38 

1.018E-8 

9.313 

788.80904 

kfum] 

39 

9.335E-7 

9.98! 

775.30535 

40 

9.S01E-7 

10.628 

762.03628 

41 

9.185E-7 

1 1.263 

749.23633 

42 

ms 

1 1.938 

736.72315 

-0.231776 

43 

3.536E-7 

12.605 

724.4171 

-0  236012 

44 

—  H  1 1 1 1  li 

13.275 

-0  240624 

45 

3.038E-7 

13.925 

"00  90944 

46 

7  7S5E-7 

14  569 

639  78353 

-0  2500751 
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TABU  B.2  CONTINUED 


A 

8 

7.542E-7 

15.261 

7.307E-7 

15.931 

7.084E-7 

16.589 

6.872E*7 

17.26 

6.669E-7 

17.918 

6.475E-7 

18.595 

6.294E-7 

19.229 

6.I24E-7 

19.905 

5.960E-7 

20.554 

5.805E-7 

21.222 

5.655E-7 

21.88 

5513E-7 

22.532 

5.378E-7 

23.216 

5.250E-7 

23.839 

5.I28E-7 

24.51 1 

BHB3B 

25.174 

4.90  IE-7 

25.825 

4.794E-7 

26.503 

4.691E-7 

27  166 

4.594E-7 

27  855 

4.503E-7 

28  489 

4  410E-7 

29.156 

4.324E-7 

29.816 

4.242E-7 

30.475 

4.J63E-7 

31.15 

4.088E*7 

31.8 

4.0I2E-7 

32.444 

3.942E-7 

33.1 13 

1 873E-7 

33.793 

J.804E-7 

34.463 

3.734E-7 

35. 1 37 

3  67 1 E*  7 

35.794 

3.61  IE-7 

36.448 

3.554E-7 

37. 1 06 

3.500E-7 

37.778 

3.447E-7 

38.439 

3.396E-7 

39.1  1 

3.348E-7 

39.76 

3.302E-7 

40.427 

3 256E-7 

41.063 

3  2ME-7 

41.722 

3.172E-7 

42.394 

3. 1 3  IE*  7 

43.029 

3.091E-7 

43.715 

3.053E-7 

44.387 

3.016E-7 

45.043 

2  979E-7 

45.693! 
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2  943E-7 


2 910E-7 


2.877E-7 


2.844E-7 


2.8I3E-7 


2.782E-7 


2.753E-7 


2.725E-7 


2.697E*7 


2.671E-7 


2.645E-7 


2.621E-7 


2.S98E-7 


2.575E-7 


2.554E-7 


2.534E-7 


2.5 1 4E*7 


2.494E-7 


2.476E-7 


2.458E-7 


2.440E-7 


m* Wl 


2.407E-7 


2.39JE-7 


2.375E-7I 


2.362E-7I 


2.348E-7 


2.334E-7 


2.32IE-7 


2.308E-7 


2.297E-7 


2.285E-7 


2.274E-7 


2.263E-7 


2.252E-7 


2.241E-7 


2  23 1E-7| 


2.221E-7 


2.21  1 E*  71 


2.202E*7| 


2.192E-7 


2  183E-7 


2.!75E*7 


2.I66E-7 


2.!57E*7 


4 


2. 14  ;£*7j 


46.359 


47  0241 


47  696 


48.34 


49.0  1  I 


49.655 


50.315 


50.98 


51.661 


52.32 


52.992 


53  636 


54.295 


54.969 


55.642 


56.286 


56.96 


57.638 


58.277 


58.949 


59.628 


60.336 


60.949 


61.6)8 


62.264 


62.898 


64.088 


64.758 


65.416 


66.031 


66. 66 


6 


68.6231 


69.302 


69.932 


70.579 


71.229 


71.892 


72.53 


73.15 


73.816 


74.431 


75.1931 


77.  147! 


424.13577 


421  71895 


3EEE3Z! 


416.94878 


414.62369 


412.3848 


410.2027 


406.03963 


404.0668 


402.10598 


400.26866 


398.4769 


396.74903 


395.08846 


393.54082 


391.99084 


390.44766 


388.99518 


387  58625 


386.17226 


384  78178 


383.52576 


382.2503 


381.00867 


379.9138 


378.79252 


377.68567 


376.61856 


375.60847 


374.66032 


373.69562 


372  79373 


371.8839 


370.96992 


370.09256 


369.24057 


368.41917 


36763375 


366.83531 


366.05546 


365.28087 


36458518 


363  37037 


63.07272! 


62.46718 


EEffigEI 


■1255261 
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FREQ  RESP 


Curve  Fit 

gRlElt _ And 

Zenos 

POL.ES  O 

ZEROS  O 

1  -248.  SSI 

l.S3400K*J  1.  188k 

2  -82.  7332  ±J  783.  093 

-1.37D19K±J  1.  44754k 

3  -155.849  ±J  1.G2947K 

1.71379K-J  3.  4440 IK 

4  —240.  038  ±J  2.  51O60K 

5  -328.  79  *J  3.  4S04K 

O  -3S8.  574  *J  4.  42724k 

7  -193.  778  ±J  S.  20974K 

Tims  delay—  O. O  S  Cain— 

2.  2E*2 1  Scale-  1.  O 

figure  E.l  PR-1592  direct  measurement  of  Q,  T=-16.75°c. 
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ba.  a 


Curve  F 

rit 

Po  1 aa  And 

Zar-oa 

POLES  O 

ZEROS  O 

1  —2.  22207K 

21  1.  049 

2  -121.  302  ±J  731.  465 

-450.  802  *J  051.  184 

3  —225.  384  -J  1.  3035 IK 

778.  248  *J  1.  77795K 

4.  -885.  635  *J  2.  54Q92K 

048.  452  *J  3.  18269K 

5  -005.  81  *J  3.  7530 4k 

890.785  -J  4.  3HOOK 

8  -863.  AAA  *J  5.  01007K 

244.473  ±J  5.  74301 K 

Tim*  dal  ay—  O.  O  S  Gain—- 

-22.25m  Seal  a-  1.0 

rigors  B.2  PR-1592  direct  neasurement  of  Q,  T=-10.16°c. 


156 


Curve  F 

Polaw  And 

POLES  O 

-it 

Zeroe 

ZEROS  O 

l  -349.  007 

ooe.  8oi 

2  -BS.  1385  ij  573.  008 

-502.  74  ±j  855.  055 

3  -173.07  *J  1.  20301 K 

942.  136  ±J  1.  524.77k 

A  -280.  335  1.  O07OOK 

-352.013  ±J  2.58523k 

5  -220.  170  ±J  2.  553k 

1.  05200k±J  3.  22073k 

6  -452.  841  *J  3.  18023K 

817.328  4.  631S2k 

7  -885.  508  A.  00707k 

8  -552.  735  ±J  A.  04008k 

Tlm«  daley-  O.  O  S  Gain"- 

-84..  E*S  Seel  a*  1.0 

Figure  B.3  PR-1592  direct  measurement  of  Q,  T=-5. 67°C. 
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X-l. 4205KH* 

FREQ  RESP 
Yb— 37.  908  dB 


Curvg  Fit 
P«?l8F _ Aod,,  .  Zar.s.9. 

POLES  O  ZEROS  O 


1 

-1.  22056K 

-5.  0090 IK 

2 

-130.  153 

-292.  017 

13 

-74.  3580 

* J 

420. 403 

534.  157 

£ 

-150.  292 

±J 

880.  051 

-440.  38 

±J 

774.  197 

■e 

-290.  404 

1.  43146K 

036.  092 

1.  20937K 

e 

-280.  021 

dj 

2.  1 OSOSK 

500.  907 

±J 

2.  24287k 

7 

-320.  211 

2.  50137K 

-317.  421 

±J 

2.  32042K 

B 

-558.  007 

±-J 

3.  1 8883K 

484.  453 

±J 

3.  348 15K 

9 

-27.  1782 

3.  787 18K 

-33.  9843 

±J 

3.  78828K 

io 

-587.  205 

±J 

4.  04075K 

213.  683 

±J 

4.  21743k 

11 

190.  088 

4.  29113k 

211.  525 

=tj 

4.  47332k 

i 

Tim*  d*l ay»  O.  O  S  Gain*-1.  25S  Scale*  1*0 


Figure  E.4  PR-1592  direct  measurement  of  Q,  T=4.50 °C. 
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FREQ  RESP 


Curva  Fit. 

polQg _ And  _ 

Zanoe 

POLES  O 

ZEROS  O 

l  -071.  752 

235.  AOA 

2  218.  SOS 

437.  171 

3  -00.  73 A0  ±j  322.  043 

-207.  052  *J  590.  228 

A  —  1 OO.  1S0  ij  073.  981 

003.713  *J  1.00809K 

5  -285.  250  1.  1DSA2K 

050.187  *J  1 . 08526k 

O  — AQ8.  192  *J  1.  7037 IK 

520.  AA7  2.  73A2AK 

7  — A57.  802  ±J  2.  35K 

300.  245  -J  3.  50935K 

8  -AA8.  82  ±J  3.  D3A1K 

141.  169  ±J  3.  80A5K 

9  -313.  9A9  ±J  3.  721Q5K 

-503.  038  ±J  A.  1  AAA  1 K 

lO  183.  135  3.  80795k 

Tima  delay—  O.  □  S  Coin—* 

-2.  003K  Scale—  1.0 

Figure  E.5  PR-1592  direct  measurement  of  Q,  T=14.23°C. 
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X-23A.  A  Hz 

FREQ  RESP 
Yb— 18.  783 
CURVE  FI 

o.  of 


FxdXY  130 


CurvQ  Fit 


POL.ES 


ZEROS 


1 

-IO.  0107K 

A98. 

AA6 

2 

-A3.  3636 

±J 

252.  968 

751. 

282 

3 

-63.  3033 

±J 

328.  1A1 

-156. 

706 

A98.  253 

A 

-196.  362 

787.  A61 

A3D. 

A72 

±J 

1 . 0965 1 H 

3 

-263.  369 

±J 

1.  23A13K 

A51. 

A52 

±J 

1.  B8847K 

iO 

-301.  6 A3 

±J 

1.  79133K 

7A8.  31 

=£j 

2.  66988K 

\7 

-331.  17 

2.  A1627K 

382. 

167 

2.  83A2K 

is 

f 

— 39S.  A77 

*J 

3.  Q29A6K 

73. 

031 

±J 

3.  59495K 

Tim*  doloy"  □.  O  S  Gain-  2.  359  Scalar  1.  O 


Figure  S.C  PR-1592  direct  measurement  of  Q,  T=25.40°C. 
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X-210.  S  H* 


FREQ  RESP 
Yb— IB.  24.  dB 
CURVE  FIT 


1 

-5.  301  1  7k 

203.  077 

2 

-31.  3313 

210.  916 

-146.  036 

411.  195 

3 

-40.  8185 

456.  83 

300.  348 

621.  941 

4 

-160.  027 

722.  045 

290.  781 

1. 22433K 

s 

-230.  222 

1.  12182K 

293.  463 

±J 

1.  80606K 

e 

-201.  274 

1.  59458K 

267.  432 

±J 

2.  34425k 

7 

-322.  301 

±J 

2.  1  2139K 

205.  474 

2.  92437k 

8 

-340.  069 

2.  6 194 IK 

117.  710 

*J 

3.  10428K 

Tim®  daisy"  0.0  S  Cain  — —2.  046m  Seal  a—  1.  O 

Figure  1.7  PR-1592  direct  measurement  of  Q,  T*37.00°C 
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X-1BO.  S  Hz 


FREQ  RESP 
Yb— 15.  026  dB 


Curve  Fit 
Elslgg-  ..  A.nsj _ <Lar«?.e 


POLES  O  ZEROS  O 


1 

-23.  6577 

197.  56 

24.7.  644 

2 

-43.  8038 

±J 

421.  323 

-177.  172 

±J 

452.  319 

3 

-67.  916-4. 

±J 

638.  484 

371.  596 

±J 

663.  196 

A 

-151.  A3 

±J 

926.  971 

381.  245 

±J 

1.  30449K 

s 

-175.  643 

±J 

l.  267BK 

354.  371 

1.  95162k 

B 

-177.  762 

AJ 

1.  6BOS2K 

1.  G3454k±J 

2.  24309H 

7 

-160.  SIS 

±J 

2.  OS44K 

212.  117 

2.  62723k 

e 

-168.  365 

±J 

2.  45328K 

-147.  187 

2.  64644k 

s 

-107. 03 

2.  677B5K 

T 1 mm  deloy« 


O.  O  S  Coln“-0O7.  SK  Scale*  1.0 


Figure  E.8  PR-1592  direct  measurement  of  Q,  T=46.42°C. 
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X-104.  6  HZ 

FREQ  RESP 
Yb--13.  7X5  dB 


Cupvq  Fit 
- And - S«rqw 

POLES  O  ZEROS  O 


1 

-19.  4151 

185.  984 

-21.  1002k 

2 

-33. 5907 

±J 

393.  474 

1 88.  882 

3 

-80.  4895 

592.  397 

-182.  705 

*J 

426.  722 

4 

-84.  4992 

842.  274 

288.  78 

548.  223 

S 

-97.  3929 

±J 

1.  0725 1 k 

-212.  117 

±J 

999.  242 

,s 

-180.  788 

±J 

1.  3507 1 K 

285.  569 

±J 

1. 088 18k 

7 

-210.  336 

1.  70143K 

313.  107 

AJ 

1.  59841k 

a 

-271. 888 

±J 

2.  OS783K 

327.  824 

±J 

2.  031 18k 

9 

-275.  528 

±J 

2.  33614K 

218.  734 

2.  44387k 

io 

-80.  9125 

±J 

2.  72989k 

-55.  8713 

2.  73418k 

11 

19.  2138 

±J 

2.  81981k 

20.  9659 

±J 

2.  82133k 

12 

123.  508 

±J 

3.  15241k 

Tima  delay*  O.  □  S  Gain*-19.  Ol^i  Scale*  1.0 


Figure  B.9  PR-1592  direct  measurement  of  Q,  T=55.85°C. 
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fcizfkvl  d. 

FREQ  RESP 
Yb— 12.  005  dB 


Cupvq  Fit 

P.OlOT - And _ Zmcsa. 

POL.ES  o  zeros  o 


1 

-3.  03879H 

-516.  395 

2 

-14.  5652 

176.  623 

168.  584 

3 

-28.  0398 

373.  634 

-145.  559 

±J 

472.  83 

4 

-38.  8522 

±J 

558.  188 

287.  662 

±J 

498.  02 

5 

-70.  4523 

±J 

784.  576 

326.  113 

±J 

942. 796 

6 

-63.  6018 

±J 

982.  986 

-115.  065 

±J 

055.  716 

7 

-160.  787 

±J 

1.  1 80 1  OR 

362.  833 

±J 

1.  301 24k 

a 

-164.  017 

1. 47345K 

266.  773 

±J 

1. 86879K 

o 

-46.  2129 

±J 

1.  87859k 

-60.  6101 

±J 

1.  68025R 

to 

-148.  596 

1.  8797k 

209.  658 

2.  Q3393K 

Tina  dal  ay-  O.  O  S  Gain— —4.  133  Scale—  1.0 


Figure  B.10  PR-1592  direct  measurement  of  Q,  T*67.oi°c 
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x-i7A.  e  hx 

FREQ  RESP 
Yb— lO.  807  d8 


Curve 

Fit 

.rmpmpi 

_ ZT-ea 

POL.ES  O 

ZEROS  O 

1 

1. 20S83K 

15 A.  085 

2 

-14.  11 AB 

1 75.  266 

-105.  933 

3AQ.  865 

3 

-17.  679A 

±J 

358.  A5A 

22A.  096 

A38.  092 

A 

— A2.  A3Q1 

±J 

5A9.  382 

-105.  1A9 

±J 

75A.  852 

s 

-37.  AS  AS 

753.  013 

25a.  ao8 

ij 

856.  805 

e 

-86.  5606 

±J 

962.  877 

—7  A.  2A87 

1.  15128k 

7 

-52.  3113 

ij 

1.  16796k 

261.  3A6 

±J 

1.  251A1K 

a 

-1. 19508K-J 

1.  2689K 

120.  19A 

±J 

1. 54388k 

s 

-220.  560 

±J 

1.  363 Ak 

226.  AA7 

1.  69923k 

10 

-297.  2A6 

1.  61223K 

201.  572 

±J 

2.  019A7K 

1 1 

112.  AS3 

±J 

1.  66329k 

12 

-215.  929 

1.  8707K 

Tima  daley1*  □.  O  S  Go  1  n“27.  AE-*"9  ScolaB  1.  O 


Figure  2>11  PR-1592  direct  measurement  of  Q,  T=74.60°C 
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x-ieo.  o  Hz 


FREQ  RESP 
Vb— -tO.  54.0  08 


Curve  Fit 
_ And — gero* 


POL.ES  O  ZEROS  O 


t 

-1.  1SSOOK 

1 82.  008 

2 

-11.  504 

170.  750 

-140.  507 

±J 

300.  207 

3 

-17.  SOOO 

350.  331 

275.  oee 

470.  433 

4 

—28.  0103 

±J 

533.  241 

-104.  330 

±J 

745.  371 

5 

-30.  4818 

±J 

732.  250 

302.  04 

±J 

874.  078 

e 

-50.  0217 

934.  350 

-120.  903 

±J 

1. 14217K 

7 

-53.  4814 

±J 

1.  1  4138k 

281.  275 

±J 

1.  201 52k 

8 

-1  18.  151 

±J 

1.  31283K 

200.  384 

1.  0008 K 

O 

-120.  233 

±J 

1.  5583 1  K 

-122.  202 

±J 

1.  77031 K 

to 

-70.  700 

±J 

1.  7023 IK 

185.  811 

±J 

2.  0098 1 K 

11 

-lie.  350 

±J 

1.  93914k 

Tims  daloy 

O.  O  S  Cain™ 

-2.  48K  Seals-  1.0 

riguro  B. 12  PR-1592  direct  measurement  of  Q,  T=84.60°C. 
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APPENDIX  F.  FNMA  DATA 


CONTENTS: 

-PUMA  FIRST  TORSIONAL  MODE  DATA  (TABLE  F.l) 
•PUNA  FIRST  FLEXURAL  MODE  DATA  (TABLE  F.2) 
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SHEAF  MODULI) 


2.13E-9I 


2  1 2E*  9| 


2.  1  2E*9| 


2. 1 2E*9 


2. 1 2E*9 


2.  1  1  E*9l 


2.J  1  E»9 


2.10E-9 


2. 1 0£*9l 


2.10E-9 


2.09E-9 


2.09E-9 


2.08E-9 


2.08E-9 


2.08E-9 


2.07E-9 


2.07E-9 


FIRST  TORSIONAL  MODE  DATA 


2.05E-9! 


2.04E-9 


2.04E-9 


204E-9I 


2.03E-9 


2.03E-9 


2.02E-9 


2.02E-9 


2.0 1  E-9 


2.0 1  E*9 


2.00E-9 


2.00E-9 


1.99E-9 


1.99E-9 


1.98E-9 


1.98E-9 


1 .97E-9I 


1.97E-9 


1 .96E-9 


1.96E-9 


1  95E*9 


1.95E-9 


1.94E-9 


1  94E-9 


1.9-5L-9 


1 .9  jE*9| 


-13.584 


-12.428 


-1 1.962 


-1 1.585 


-11. 153 


-10.682 


-10.161 


-9.579 


-9.017 


-8.395 


-7.799 


-7.305 


-5.958 


-5.363 


-4.738 


-4.273 


-3.593 


-2.972 


-1.679 


-1.01 


0.206 


0.7721 


1.527 


2.2231 


2.844 


3.487 


4.121 


4.771 


5.385 


6.054 


6.718 


7.367 


8.052 


8.703 


9.367 


10.028 


10.688 


1 1.349 


12.013 


12.665 


13.322 


13.984 


14  644| 


C 


RES  FREO 


2170.2295 


2167.273 


2165.391 


2163.971 1 


2162.4858 


2160.8468 


2159.0281 


2157.0361 


2155.0013 


2152.9106 


2150.7247 


2148.7597 


2146.365 


2144.058 


2141.8177 


2139.6827 


2137.7027 


2135.4514 


2133.1573 


2131.3394 


2 1 28.6048 


2126.0178 


2123.7212 


INPHASE 


-0.186175 


-0.18587 


-0.186084 


-0  185038 


-0.184178 


-0.18234 


■OTWTHIIl 

M-lflllH 


-0.179076 


TUFfl-VHl 


-0.176412 


-0.175547 


-0.17404 


-0.17181 


-0.17015 


mmiL 


21  19.1345 


21  16.5967 


21 13.8699 


21  1 1.2689 


2108.8021 


2106.3344 


2103.8559 


2101.4502 


2098.8933 


2096.3149 


2093.7681 


2091.2058 


2068.6161 


2085.9914 


2083.3422 


2080.6743 


2077.9868 


2075.329 


2072.6075 


2069.8903 


2067.138 


2064.3723 


-0.16598 


-0.16653 


-0.163175 


-0.161541 


-0.160346 


■dESEEDl 


-0.157669 


-0.156204 


-0.154686 


-0.153449 


-0.151969 


-0.150656 


-0.149201 


-0.147574 


-0.146263 


-0.144589 


-0.143087 


-0.141561 


-0.140178 


-0.137134 


-0.13566 


-0.134306 


TABLI  7.1  CONTINUED 


latel 

A 

8 

C 

_ P  1 

47 

1  92E*9 

15  3 

MIlK-iE 

46 

1  92E-9 

15.948 

2056.7444 

-0.13157 

49 

1.91E-9 

16.615 

2055.9833 

-0.130222 

50 

1.91E*9 

17.272 

2053.1806 

mpmmi 

51 

1.90E-9 

17.93 

2050.4226 

■EXFiZEE 

52 

1.90E-9 

18.583 

2047.6075 

■ramaa 

53 

1  89E*9 

19.249 

2044.8976 

-0.1247181 

54 

1.89E-9 

19.892 

2042.0883 

■CTEgEE 

55 

-  - 

1  88E*9 

20.557 

2039.3161 

-0.1219 

56 

1.88E-9 

21.21 

2036.5754 

-0.120365 

57 

1.87E*9 

2t.856 

■EFXVT^i 

58 

1.87E*9 

22.495 

2031.0002 

59 

1  86E*9 

23.161 

2028.0434 

-0.1 15529 

60 

1.86E*9 

23.808 

2025.1252 

-0.1  14124 

61 

1  85E-9 

24.451 

2022.308 

-0.1  12963 

62 

-  -1.84E-9 

25.104 

2019.4712 

-0.1  12031 

63 

1.84E-9 

25.762 

2016.4934 

-0.1  10583 

64 

1  83E*9 

26.419 

2013.4772 

-0.109642 

65 

1.83E-9 

27.08 

■K1TOCHE 

-0.108044 

66 

1  82E*9 

27.74 

2007.6932 

67 

I.82E-9 

28.392 

2004.729 

-0.104024| 

68 

1 .8 1  E*9 

29.045 

2001.7857 

■BMHW 

69 

1 .8 1  E-9 

29.697 

1998.8531 

-0.101044 

70 

1.80E*9 

30.368 

1995.884 

-0.09968 

71 

1.80E-9 

31.015 

1992.903 

-0  09851  1 

72 

1.79E-9 

31.669 

1989  9731 

-0.097088 

73 

1.79E*9 

32.324 

1987.0484 

-0.095585 

74 

1.78E-9 

33.007 

1983.9916 

■BXB3ES 

75 

1.78E-9 

33.659 

1980.9691 

-0.092622 

76 

1.77E-9 

34.328 

1977.9937 

-0.091032 

77 

1.76E-9 

34.994 

1975.0318 

-0.08946 

78 

1.76E-9 

35.665 

1972.0582 

79 

■■ESS 

36.32 

1969.1  1 1  1 

-0.086077 

80 

1.72E-9 

36.974 

1948.335 

-0.092382 

81 

1.71  E*  9 

37.644 

1944.5879 

-0.091  188 

82 

1.72E-9 

WMVlVM.m 

-0.088995 

83 

1.71E-9 

38.98 

1943.7924 

-0.087663 

84 

1.70E*9 

39.642 

1940.21  1 

85 

1.70E-9 

40.317 

1936.5939 

-0.0851981 

86 

1  69E*9 

40.964 

1933.0179 

■BUFFET) 

87 

■d^e 

41.63 

1929.4515 

11-IH  FFftil 

88 

1.68E-9 

42.282 

1925.8672 

-0.081  149 

89 

1.67E‘9 

1922.2206 

-0.079857 

90 

1  66E-9 

43.636 

1918.5008 

-0.078538 

91 

1.66E-9 

44.284 

1914.7348 

-0.077124 

92 

1.65E-9 

44.953 

191 1.0058 

-0.075723 

93 

1.65E-9 

45.61  1 

1907.1974 

-0.074261 

169 


TABLE  7.1  CONTINUED 


I.64E-9I 


I  63E-9 


I  62E-9 


1.62E-9! 


I.5IE-9 


I.60E-9I 


46.272 


46.941 


47.606 


48.269 


48.941 


49.607 


50.251 


50.895 


52.2751 


52.945 


53.616 


1903.14641 


1899.0635 


ezesei 


1890.032! 


1885.6676 


1881.5396 


1377.5627 


1873.4947 


1369.3214 


1365.1634 


1861.0747 


1856.98861 


1352.3529 


1848.7837 


1344.7561 


1840.93341 


1337.1454 


1333.4599 


1829.993 


1826.443 


1823.1  1  15 


1819.8698 


1816.8279 


1813.3706 


131  1.0658 


1808.8188 


1806.9164 


1805.5463 


8049857 


1305.4492 


1807.72551 


1313.451 


1826.1035 


1851.6973 


1894 1251 


1936.0094 


1964  325! 


1991.0099 


1998  4142 


2004.34181 


2009  391 


7248 


-0  0727941 


071298 


-0. 070089 


-0.069288 


-0.0686  18i 


-0.067665! 


-0.0665661 


-0.065518 


-0.064159 


-0.061554 


-0.0601051 


-0.0586645 


-0  057202 


-0.05573 


-0.054342 


-0  0528781 


-0.051514 


-0.050028 


-0  048531 


-0.046935 


-0.0453581 


-0.043664 


-0.041856 


-0.040127 


-0.038251 


-0.036304) 


-0.0341  1 2i 


-0.031776! 


-0.029136 


-0.0260628i 


-0.0222559! 


-0.017304 


-0.01  143991 


-0.0063043 


-0.0033283! 


-0.0020376! 


-0.00144lj 


-0.0010852! 


-0  00087791 


-0.000683  I  i 


-0.0005391! 


-0.0002596! 


-0.0001  77! 


-0  0000873! 


TABLE  7.1  CONTINUED 


Label 

A 

8 

C 

D 

14) 

1.87E-9 

77  782 

2030.5228 

0  0001071 

142 

1.87E-9 

76.431 

2033.3167 

0  000 1963 

143 

1.87E-9 

79.083 

2030.7322 

0.0002269 

144 

1.87E-9 

79.761 

2032.6043 

0.0002942 

145 

1.87E-9 

80.415 

2035.2389 

0.0003895 

146 

1.88E-9 

81.085 

2037.8902 

0.0004928 

147 

1  88E*9 

61.824 

2040.4153 

0.0005475 

148 

1.89E*9 

82.58 

2043.1523 

0.0006504 

149 

1  89E*9 

83.231 

2045.2337 

0.0007243 

150 

1  90E*9 

83.886 

2047.54091 

0.0008198 

151 

1.90E-9 

82.851 

2047  6868 

0  000768 
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T&BL1  7.2 


label 


label  iPttiiA 


1ST  FlEXllftAl 


DATA 


m*mrz.zwi 


6.091E* 


6.073E* 


6.059E* 


6.048E- 


6.038E* 


6.028E* 


6.018E* 


6.004E* 


5.990E* 


5.979E* 


5.964E- 


5.949E* 


5.933E* 


5.9I9E- 


5.907E* 


5.892E* 


5.879E* 


5.865E* 


5.848E* 


5.832E* 


5.81  7E* 


5.800E* 


5.784E* 


5.767E* 


5.753E* 


5.738E* 


5.723E* 


5.705E* 


5.689E- 


5.672E- 


5.657E- 


5.640E- 


5.623E* 


5.607E* 


5.590E- 


5.574E* 


5.557E- 


5.540E- 


5.523E* 


5.507E- 


5.49  IE* 


5.474E* 


5.439E-9 


5.423E*9 


5  406E-9 


(C) 


-14.252 


-13  115 


-12.561 


-12.163 


-11. 794 


-11.349 


-10.83 


-10.282 


-9.645 


-9.102 


-8.509 


-7.876 


-7.281 


-6.703 


-6.192 


-5.615 


-4.996 


-4.387 


-3.702 


-3 


-2.445 


-1  674 


-1.092 


-0.414 


0.204 


0.838 


1  508 


2.176 


2.848 


3.479 


4. 1  17 


4.767 


5.419 


6.08 


6.75 


7.395 


8.102 


8.766 


9.434 


10.095 


10.76 


11. 426 


12.091 


12.754 


13.407 


14.087 


RES  FREQ  (Hz) 


271.49919 


271.09085 


270.77732 


270.54238 


270.31047 


270.0915 


269.8718 


269.54729 


269.23162 


268.98782 


268.64069 


268.30832 


267.95162 


267.63059 


267  37337 


267.02394! 


266.71982 


266.40057 


266.01735 


265.6565 


265.319 


264.92506 


264.56703 


264.18612 


263.86577 


263.5001 


263.16049 


262.7626 


262.381 14 


261.98618 


261.64662 


260.860321 


260.48888 


260.09807 


259.72091 


259.31066 


258.931551 


258.53054 


258.14459 


257.77471 


INPHASE  VOLTA 


-0.47969 


-0.48558 


-0.48373 


-0.47379 


-0.46971 


-0.45769 


-0.45504 


-0  44861 


-0.4461  1 


-0  44244 


-0.43413 


-0.42969 


-0.42585 


-0.4216 


-0.41785 


-0.4133 


FI  I 


-0.3942 


-0.3899 


m 


-0.37568 


-0.37125 


-0.36638 


-0.36183 


-0.3567 


-0.34804 


-0.34297 


-0.33888 


-0.33474 


1 1 1 1 1 11 1 


256.972451 


256.56169 


256.17371 


255.774061 


-0.32491 


-0.32122 


-0.31646 


-0.31469 
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TABLE  F.2  CONTINUED 


Lat>*i 

A 

8 

C 

D 

47 

5  389E-9 

14.74 

255  37819 

-0.31019 

48 

5.372E-9 

15.4 

254.95939 

-0  30622 

49 

5.354E*9 

16.088 

254.5452 

-0.30175 

50 

5.337E-9 

16.742 

254.14721 

-0.2983 

5! 

5.320£*9 

17.399 

253.73427 

-0.29445 

52 

5.303E*9 

18.058 

253.31341 

-0  29004 

53 

5.285E*9 

18.736 

252.88746 

■EX3331I 

54 

5.267E*9 

19.385 

252.47315 

-0.26265 

55 

5.250E-9 

20.067 

-0.2793 

56 

5.234E*9 

20.732 

251.66251 

■BEEiSE 

57 

5.2 1 6E*9 

21.379 

251.25042 

■B3EEE3-] 

58 

5.200E-9 

22.045 

250.85408 

-0.268172 

59 

5.183E-9 

22.692 

250.43 1 5 

-0.264334 

60 

5.  !65E*9 

23.379 

250  00432 

-0.260238 

61 

5.148E-9 

24.032 

249.58872 

-0.256662 

62 

5.13  IE-9 

24.693 

249.1881  1 

63 

5. 1 1 5E-9 

25.36 

248.78806 

64 

5.097E-9 

26.017 

248.35937 

-0.2425461 

65 

5.080E-9 

26.702 

247.94171 

■EffHIB 

66 

5.062E-9 

27.354 

247.51388 

"0.232402 

67 

■OSEE3E 

28.001 

247.00949 

-0.226773 

68 

5.01  3E*9 

28.686 

246  3)319 

-0.224406 

69 

4.983E-9 

29.358 

245.56749 

-0.227018 

70 

4.962E-9 

30.006 

245.04262 

-0.227396 

71 

4.°42E-9 

30.68 

■ESSES 

-0.224925 

72 

4.922E-9 

31.351 

244.06237 

-0.223304 

73 

4.904E-9 

31.997 

243.60181 

-0.22)045 

74 

4.885E-9 

32.687 

243.13515 

-0.219335 

75 

4.S66E-9 

33.349 

242.65932 

-0.216672 

76 

4.848E-9 

34.034 

242.20496 

-0.21464 

77 

4.829E-9 

34  684 

241.7266 

-0.212782 

78 

4.8I0E-9 

35.352 

241.26388 

■HHKtt) 

79 

4.790E-9 

36.035 

240.74994 

-0.210961 

80 

4.776E-9 

36.702 

240.41202 

-0.197931 

81 

4.764E-9 

37.36! 

240.09571 

-0.199636 

82 

4.731E-9 

38.027 

239  26357 

-0.195996 

83 

4.71 3E-9 

36.697 

238.82677 

-0  196073 

84 

4.695E-9 

39.388 

238.35665 

-0.195201 

35 

4.676E-9 

40  044 

237.86973 

-0.192627 

86 

4.656E*9 

40.702 

237.36183 

87 

4.637E-9 

41.375 

236.879 

-0.187734 

88 

4.6 1  7E*9 

12.035 

236.3735 

-0.135289 

89 

4.598E-9 

42.708 

235.89405 

-0.183189 

90 

4.579E-9 

43.36 

235.40172 

-0.180735 

91 

4.560E-9 

44.045 

234  90971 

-0.17835 

92 

4.541  E-9 

44.698 

-0.1  ^5398 

93 

4.521E-9 

45.348 

233.91512 

-0  174064 
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TAB LI  F.2  CONTINUED 


4.503E-9I 


4.482E-9 


4.463E-9 


4.444E*9 


4.424E*9I 


4.40SE*9 


4.385E-9 


4.365E*9| 


4.346E*9| 


4.327E*9| 


4.354E*9l 


4.339E*9| 


4.323E*9 


4.308E*9 


4.292E-9 


4.278E*9 


4.265E*9 


4.253E*9 


4.24  IE*  M 


4.230E*9 


4.22 1E*9| 


4.2 1 2E*9| 


4.205E*9 


4.200E*9 


4.20  lE*9f 


4.21 0E*9 


4.224E*9 


4.253E*9| 


4.297E*9 


4.363E-9 


4.429E-9 


4.495E-9 


4.567E-9 


4.652E*9| 


4.735E-9 


4.760E*9 


4.8S0E-9 


4.966E-9 


4.998E-9 


4.687E-9 


4.668E*9 


4.71  IE*9| 


4.920E*9 


5.502E-9 


5.59JE-9 


5.610E-9 


6  ( 


46.018 


46.691 


47.367 


48.024 


48.685 


49  356 


50.036 


50.701 


51  378 


52.063 


52.744 


53  399 


54.053 


54.75 


55.405 


56.088 


56.764 


57.431 


58.1 13 


58.773 


59.443 


60.137 


60.789 


61.465 


62.142 


62.806 


63.737 


64.58 


65.248 


65.885 


66.548 


67.179 


67.796 


68.434 


69.077 


69.736 


70.421 


71.095 


71.303 


72.454 


73.101 


73.781 


74.45 


75  102 


75.601 


76  468 


77  166 


233  43427 


232.89825 


-0.171791 


-0.169821 


"i 


231.90219 


231.3790 


230.88204 


230.36055 


229.83951 


229.32451 


228.81621 


229.54222 


229.15715 


228.72323 


228 33618 


227.89321 


227.53613 


227.17343 


226.85622 


22653992 


226.25756 


225.99864 


225.77141 


225.59249 


225.44929 


225.44688 


225.48359 


225.70207 


226.09955 


226.85468 


228.0447 


229.7831  1 


231.51  122 


233.21795 


235  10142 


237.26624 


239.38204 


mmm\ 


243  01838 


245.14248 


245.935921 


238.15062 


237.67969 


278.7577 1 


243.99396 


258  04527 


260. 11  36 


260.5576 


-0.165549 


lam 

TUnTTFll 


-0.155281 


MBHEEZE 


-0.107914 


-0.103455 


-0.09865 


-0.091179 


-0.088552 


-0.083125 


-0.076889 


-0.073502 


-0.069912 


-0.066015 


-0.061 175 


-0.056175 


-0.050956 


-0.035828 


-0.03804 


-0.0235041 


-0.027936 


-0.052924 


-0.035818 


-0.036457 


-0.0250032 


-0  0240697 


0.0239374 


TABLE  7.2  CONTINUED 


■PS 

Label 

A 

8 

C 

D 

Ml 

5.62  IE-9 

77.742 

260  82056 

-0.0236572 

M2 

5.630E-9 

78.336 

261  02845 

-0.0233475 

143 

5.640E-9 

79.077 

261.24059 

■B.mnin 

144 

5.650E-9 

79.784 

26 1. 48 127 

WEE3333L E 

1 45 

5.659E-9 

80  501 

261.68881 

-0.0224509] 

146 

5.669E-9 

81.273 

261  93298 

147 

5.679E-9 

81.969 

262 14919 

-0.0220627 

>46 

5.688E-9 

82.66 

262.36215 

-0.0216973 

149 

5.697E-9 

83.321 

26255937 

-0.0215842 

150 

5.706E-9 

83.999 

262.77687 

151 

5.71 3E*9 

82.698 

262.93856 

-0  0213681] 
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APPENDIX  0.  POLYCARBONATE  DATA 
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-  POLYCARBONATE  FIRST  TORSIONAL  MODE  DATA  (TABLE  6.1) 

-  POLYCARBONATE  FIRST  FLEXURAL  MODE  DATA  (TABLE  G.2) 

-  POLYCARBONATE  FREQUENCY  RESPONSE  CURVES  AND  POLE  ZERO 
INFORMATION  FOR  BOTH  THE  TORSIONAL  AND  FLEXURAL  MODES  OF 
THE  HOLLOW  SAMPLE  ROD  (FIGURES  G.l  AND  G.2) 
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TABLE  0.1  POLYCARBONATE  FIRST  TORSIONAL  MODE  DATA 


LaMl  a  IB  ! 


LSMl  IPOLYCTI  SHEAR  (-OOULUjTEMP  iCi 


1ST  TORSIONAL  9 724E-8 


9.695E-8  -13.2321 


9.675E-8)  -12.671 


9.66  IE-3!  -12.221 


9.65 1 E- 31  -11.9091 


9.640E-8  -11.47 


9.626E-8I  -10.983 


9.609E-8I  -10.435) 


9.59 1 E-81  -9.7721 


9.575E-8I  -9.269 


9.553E-8  -3.673! 


9.54 1  E-3|  -3.024i 


9.525E-3I  -7.472! 


9.477E-3I  -5.7411 


9.464E-3I  -5.208 


9.445E-3I  -4.425 


9.430E-3!  -3.772! 


9.41 5c-3i  -3.141 


9.400E-.8I  -2.5411 


9.386E-8I  -1.872 


9.372E-3I  -1  2«w 


9.359E-3 


9.347E-8I  0.03 


IRES  FREQ  I  INPHASE 


9.322E-8I 


9.310E-3 


9.298E-3 


9.277E-8I 


9.266E-5! 


9.255E-  3| 


9.246E-3I 


9.236E-  3! 


3.978 


4.645 


5.3021 


5.964| 


6.635! 


1221.36923 


1220.75375 


1220.0566 


1219.15391 


1213.12276 


1216.94905 


1215.9571 1 


1214.84478 


1213.76645 


1212.78635 


121  1.70231! 


1210.72757 


1209.70394 


1206.36828 


1207.67921 


1206.674491 


1205.71632 


1 204.78427 


120386839 


1203.01206 


1202.17127 


1201.34901 


1200.561 1 


1 199.76499 


1 199  01288 


1 198.24817 


1 197.54846 


1 196.35019 


1 196.14473 


1  195.47763 


1 194.34059 


1 194.18978 


Inmi 

■BggZHi 


-0.121538 


-0.122435 


9.2 !  5E-3 

7.943  1192.92552 

9.206E-3 

3.645  1192.30088 

9.197E-8 

9.308  1191.6863 

9.188E-3 

9.976  1191.08726 

9  179E-3 

10.631  1190  4998 

-0.123072 


-0. 1 35792 


ESEB5I I 


-0.133279 


-0.153375 


-0.133156 


-0.138401 


-9.134883 


EIESZSfl 

[■MEii-q 

Milifri'i-il 

mnm 


-0.164835 


-0.154226 


-0.211296 


1 139  2206! 


-0.171644 


-0.172017 


DXCiESHEXEEEEEl 


TABLE  0.1  CONTINUED 


■E 

Libel 

A 

b 

c 

D  1 

47 

9.123E-3 

14.623 

1 186.89373 

■gram 

48 

9.1 14E-8 

15.302 

1 186.28628 

EBEEEEH 

49 

9.103E-8 

15.961 

1185.60096 

-0.181255 

50 

9.094E-8 

16.634 

1 1 84.99556 

-0.183068 

51 

9.085E-8 

17.294 

1 184.42276 

-0.189827 

52 

9.076E-8 

17.957 

1183.32581 

53 

9.066E-8 

18.615 

1183.20417 

E'XIIUC 

54 

9.057E-8 

1 1 82.59345 

m&mzki 

55 

9.047E-3 

19.953 

1 181.95599 

EfBI 

56 

9.037E-8 

20.617 

1 181.3102 

57 

9.02SE-3 

21.266 

1 180.6884 

58 

9.019E-8 

21.928 

1 180.09464 

-0.216299 

59 

9.010E-8 

22.588 

1179.48887 

-0.220503 

60 

9.000E-8 

23.254 

1  178.8623 

BM5-] 

61 

8.99 1 E-S 

23.917 

11 78.24662 

62 

8.981E-3 

24.572 

1 1 77.64259 

BM 

63 

8.972E-8 

25.227 

1 177.01261 

-0.2324841 

64 

3.96  IE-8 

25.915 

1 176.31068 

Ml  1  1  PI  1 

65 

8.951E-8 

26.586 

1 175.66305 

atm 

66 

8.942E-3 

27.252 

1175.03053 

67 

8.932E-8 

27.917 

1 174.39662 

-0.302391 

63 

8.922E-8 

28.59 

1173.75079 

69 

8.912E-8 

29.269 

1173.09212 

70 

8.903E-8 

29.93 

1172.49072 

■BJRTF1 

71 

8.894E-8 

30.586 

1171.88196 

72 

8.884E-8 

31.263 

1171.22765 

■b m 

73 

31.921 

1170.62065 

-0.282191 

74 

6.863E-3 

32.6 

1169.86199 

75 

8.355E-3 

33.253 

1 169.31279 

-0.29601 

76 

8.344E-3 

33.946 

1 168.61764 

-0  3«732 

77 

8.335E-8 

34.6 1 1 

1168.02641 

BETtm 

78 

8.825E-8 

35.289 

1167  3726 

■SM'ZH-i 

79 

8.319E-B 

35.96 

1166.93232 

80 

8.509E-8 

36.625 

1 166.271 15 

-0.318291 

81 

8.799E-3 

37  305 

1 165.65739 

■HHFIFti 

82 

3.790E-3 

37.967 

1165.01225 

KVHEi'l  v'.T.I 

83 

8.730E-3 

38.651 

nom 

-0.323091 

84 

8.771E-3 

39.316 

1 163.75581 

wnmm 

85 

3.761E-3 

1 163.1  155 

86 

3.752E-6 

40.668 

1 162.49039 

■HWU  1 

87 

3.744E-3 

41.32 

1161.95399 

RHECS&EE2 

88 

3.738E-8 

41.985 

1 161.57054 

89 

3  716E-3 

42.654 

1 1 60  1  1262 

-0.328061 

90 

5.725E-5 

4J.  3H 

■ar-MBTim 

EB-FR-tH 

91 

5.707E-3 

43  975 

1  159.465851 

-0.322261 

92 

5.682E-  5 

44  549 

wirri'lil; 

93 

5.sT5E-5 

-0  31346 
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TAELS  G.l  CONTINUED 


label 

A  ! 

8  1 

C 

_ 0 _ I 

94 

8.653E-8 

45.9841 

■gaHEE 

SUSjflTFE 

95 

8.63 1  E-8i 

46.6581 

1 154  43383 

■■fncrrm 

96 

8.626E-8j 

47.3181 

1 154.07846 

97 

8.687E-8) 

47.9991 

1 158.17008 

■HTTCTinrD 

98 

8.648E-8I 

48.664) 

■BEES 

99 

■HHrMfl-l 

49.223) 

1 152.2222 

-0.3189) 

100 

S.588E-6) 

49.9891 

■E53EZH 

101 

B.530£-8| 

50.6531 

1 151.0558 

MFFIE FT! 

102 

8.6I6E-8) 

51.321) 

1153.4352 

T1  fill 

103 

8.619E-8I 

51.9621 

1153.62245 

104 

8.589E-9i 

52.6411 

1  151.61744 

■BEEBES 

105 

8.60IE-5I 

53.332 

1 152.41556 

-0.30518) 

106 

8.590E-ei 

53.9961 

1 151.63543 

■EE3ZE 

107 

3.55  1 E -si 

54662! 

1 149.10699 

—m  1  a  il 

108 

8.575E-3i 

55.329) 

1 150.67028 

-0.35881) 

109 

3.484E-3 

55.973) 

■■EEffl 

110 

8.478E-8) 

56.6661 

1 1 44. 1  764 

-0.30526) 

1 1 1 

8.497E-8 

57.30el 

■IEBEME 

Trrrir 

112 

8.531E-8I 

57.979) 

1 147.73343 

■CBHE 

113 

3.542E-3I 

53.6431 

1 146.5016) 

■MU 

114 

8.520E-3) 

59.309 

1 146.97564 

1 15 

8.49  IE-31 

59.9631 

1 145.04673 

HBV(.r  r->i 

1 16 

8.425E-3I 

50.6361 

1  140.58997 

■EMC 

117 

3.436E-8I 

61.293) 

1 141.33875 

SBESUZiE 

118 

8.413E-8 

61.9451 

1 139.77071 

-0.384141 

119 

8.418E-8 

62.63 1| 

1140.10536 

■KTQFTi 

120 

8.465E-3 

63.28) 

1 143.28282 

HBfiU-VI-T.I 

121 

8.399E-8 

64.3791 

1 138.80922 

—ronni 

122 

8.398E-3 

65  038) 

1 138.74222 

-0.3702  «| 

■  4»« 

3.354E-8 

65.5931 

1 135.74022 

-0.37692 

124 

3.414E-8 

66.258) 

1 139.8722 

■BTCTTE 

125 

3.406E-8 

66.999) 

1139.21633 

126 

3  212E-3 

67.5491 

1 132.93717 

-0.368421 

127 

3.380E-3 

66  302) 

1 127  55476 

-0.37449) 

128 

8  342E-8 

63.974! 

1134 97596 

-0.33204 

129 

3.358E-3 

69.612) 

1 136.01375 

-0.371221 

130 

8.337E-8 

70.2391 

-0.373521 

131 

■■BXEEU 

1 1 36.78 1 47 

■■EHrTlCH 

132 

3.309E-3 

71.5671 

1  132.70025 

Tlillil 

133 

3.310E-3 

72.232) 

'  132.77157 

-0.36628 

134 

8.255E-8 

72.39li 

1 129.01553 

-0.37203 

135 

3.235E-3 

72  526) 

-0  27161 

136 

8.235E-3 

74  197! 

1127  54286 

-0.38178 

13“ 

5.223E-3 

74.544) 

1126.55469)  -0.27833) 

133 

S  2 ' 3E-3 

-5  ^gei 

1  126  156461  -0.37793) 

139 

5  Zzlim 5 

*C  '7{ 

l 5  30  94534 

-0.38) 

140 

■BHBE— iSEI 

1  '■  24  5-:929i  -0  354331 
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TABLE  0.1  CONTINUED 


8. » 796*6 

6.1606*6 

8.1976*6 

6.2416*8 

6.1466*6 

6.1576*6 

8.1336*8 

8.1556*6 


TABLl  0.2  POLYCARBONATE  FIRST  FLEXURAL  MODE  DATA 


Label 

A 

s  i 

_ C _ 

0 

Label 

POLYC 

YOUNG'S  (CORF 

TEMP  (C)  I 

INPHASE  VOLTa 

i 

FLEXURAL 

2.642E-9 

127  000391| 

-0.21429 

2 

2.637E-9 

-12  328 

wmssmi 

-0.222156 

3 

3-19-93 

2.632E-9 

-1  1  731 

126.773926 

-0.228169 

4 

2  629E*9 

-1 1.29 

126.696091 

-0.231512 

5 

2.627E-9 

-10.936 

126.639234 

-0.233968 

6 

2.624E-9 

-10.546 

126.564049 

-0.236452 

7 

2.622E*9 

-10.132 

126.516773 

-0.237755 

8 

2.618E-9 

-9.648 

126.436356 

-0.227368 

9 

2.615E-9 

-9.141 

126.353962 

MPWTlTfl 

10 

2.61  1  E*9 

-8.583 

126.267921 

BBBEBE 

1  1 

2.608E-9 

-8  046 

126.192124 

-0.2363281 

12 

2.605E-9 

-7.49 

126.107905 

^ 11 

13 

2.60  IE-9 

-6  992 

126.026784 

-0.246815 

14 

2.598E*9 

-6.451 

125.951 162 

-0.249454 

15 

2.595E*9 

-5.869 

125.861302 

-0.252972 

16 

2.59 1 E  -  9 

-5.167 

125.7651 

-0.257177 

17 

2.587E-9 

-4.632 

125.68581 1 

-0.259309 

18 

2.584E-9 

-4.009 

125.595644 

-0.26481 1 

19 

2.580E-9 

-3.327 

125.502299 

20 

2.577E-9 

-2.792 

125  426931 

-0.267575 

21 

2.573E-9 

-2  157 

125.331196 

-0.268595 

22 

2.569E-9 

-1.461 

125.240551 

-0.27233 

23 

2.565E*9 

-0.791 

125.154687 

24 

2.562E*9 

-0.23 

125.068641 

-0.279032 

25 

2.559E-9 

0.514 

124.984882 

-0.264202 

26 

2.555E-9 

1.188 

124.904864 

-0.283489 

27 

2.553E-9 

1.794 

124.842905 

-0.286021 

28 

2.549E-9 

2.486 

124.757551 

-0.288916 

29 

2.546E*9 

3.133 

124.679181 

-0.291323 

30 

2.543E-9 

3.72 

Mmmm 

-0.298152 

31 

2.539E-9 

4.421 

124.517856 

-0.30487 

32 

2.536E-9 

5.025 

124.439624 

-0.30729 

33 

124.368327 

-0.30638 

34 

2.530E-9 

6.348 

124.296501 

-0.31  12 

35 

2.528E-9 

6.993 

124.227743 

■K1ZEEE3 

36 

2.525E*9 

7.624 

124.159476 

-0.32253 

37 

2.522E-9 

5.314 

124.093972 

-0.32685 

38 

2.520E-9 

8.996 

124.028582 

-0.3315 

39 

2.517E-9 

9.616 

123.964034 

-0.33747 

40 

2.514E-9 

10.273 

123.89981 1 

-0.34123 

41 

2.51  1E*9 

10.914 

123.827324 

-0.34563 

42 

2.509E-9 

1  1.56 

123.769061 

-0.351  1 

43 

2.506E-9 

12.239 

123.699074 

-0.35691 

44 

2.504E-9 

12.868 

123.639647 

-0.36307 

45 

2.501E-9 

13.524 

123.576923 

-0.36493 

46 

2  499E*9 

14.2 

123.522434 

-0.37141 
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TABU  0.2  CONTINUED 


■El _ Lab«l _ 

A 

B 

c 

D 

47 

2.496E-9 

-0.38072 

48 

2.494E*9 

15.54 

123  40346 

-0.38335 

49 

2.492£*9 

16  181 

123.348461 

-0.38767 

50 

2.490E-9 

16.866 

123.294596 

-0.39282 

51 

2.487E*9 

17.519 

123.235722 

-0.39996 

52 

2  485E-9 

18.167 

123.176167 

-0.39963 

53 

2.483E-9 

18.85 

123.1 15736 

-0.40691 

54 

2.48 1  E*9 

19.492 

123.071968 

-0.41  1 

55 

2  479E-9 

20.176 

123.017078 

-0.41765 

56 

2.476£*9 

20.839 

122.953831 

-0.4266 

57 

2.475E-9 

21.509 

122.916744 

-0  4301  1 

58 

WMESESL 

22.139 

122  866088 

-0.4353 

59 

2.470E*9 

22.823 

122.6141 12 

-0.44086 

60 

2  468E*9 

23.478 

122.747125 

-0.44543 

61 

2.466E-9 

24.15 

122.703381 

-0  45089 

62 

2  464E*9 

24.796 

122.644444 

-0.45669 

63 

2.462E‘9 

25.472 

122.596565 

-0.46061 

64 

2.459E-9 

26.152 

122.539603 

-0.46553 

65 

2.457E*9 

26.807 

-0.47061 

66 

2  455E*9 

27.463 

122.422573 

-0.47455 

67 

28.128 

122.365124 

-0.48164 

68 

2  451E-9 

26.769 

122  325329 

-0.48608 

69 

2.447E-9 

-0.49149 

70 

2  445E-9 

30.127 

122.176972 

-0.49603 

71 

■EEEBE 

30  778 

122.127921 

-0.49936 

72 

2.440E-9 

31.451 

122.047595 

-0.50246 

73 

2.435E-9 

32.1 16 

121.932692 

-0.4972 

74 

2.43 1E*9 

32  792 

121.823426 

-0.49018 

75 

■KIES3 

33  462 

121.921701 

-0.48777 

76 

2  441E-9 

34.126 

122.079184 

-0.48064 

77 

2.440E-9 

34  796 

122.055235 

-0.48739 

78 

2.429E-9 

35.451 

121.777248 

-0.48912 

79 

2.41 8E*9 

36.126 

121.507604 

-0.481  1  1 

80 

2.430E-9 

36.784 

121.801276 

-0.4793 

81 

2.434E-9 

37.428 

121.914386 

-0.48025 

82 

2.410E-9 

38.106 

121.300488 

-0.47127 

83 

2.407E-9 

38.772 

121.240164 

-0.47712 

84 

2.430E-9 

39.417 

121.808303 

-0.47637 

85 

2.41 9E*9 

40.082 

121.538649 

-0.47266 

86 

2.40  !E*9 

40.76 

121.084137 

-0.47439 

87 

2.41  4E*9 

41.44 

-0.481  16 

88 

2  415E-9 

42.085 

121.433903 

-0  48265 

89 

2.404E*9 

42.781 

121  13976 

-0.48636 

90 

2  402E*9 

43.425 

121.1  11747 

-0.48191 

91 

2.396E-9 

44.063 

120.947438 

-0.47932 

92 

2.395E*9 

44.742 

120.927635 

-0.4975 

93 

2.390E*9 

45.405 

120.788192 

-0.483 
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TABLE  0.2  CONTINUED 


WF 

Label 

A 

5  | 

C 

D 

94 

2.389E-9 

46  077 

120  763558 

-0.50156 

95 

2  388E-9 

46  7441 

120.75862 

-0  49943 

96 

2.400E-9 

47.386 

121.04571 1 

-0  49404 

97 

2  395E-9 

48.072 

120.924066 

-0.47946 

98 

2  400E-9 

48.735 

121.058061 

-0.4789 

99 

2  394E*9 

49.392 

120.908332 

-0.48401 

100 

2.388E*9 

50.064 

120.756667 

•0  48031 

101 

2.378E*9 

50.716 

120.490645 

-0.48812 

102 

2.374E*9 

51.376 

120.397917 

-0.4964 

103 

2.373E*9 

52.051 

120.375247 

-0.48714 

104 

2.382E-9 

52.71 1 

120.609069 

-0.4934 

105 

2.377E*9 

53.365 

120.480052 

-0.49199 

106 

2.380E-9 

54.037 

120.539759 

-0.48619 

107 

2  382E-9 

54.697 

120.585249 

-0.47884 

108 

2.379E*9 

55.365! 

120.516678 

-0.46473 

109 

2.379E-9 

56.024 

120.522072 

-0.46798 

110 

2  376E-9 

56.688J 

120.451713 

-0.49324 

1  1  1 

2.364E-9 

57.371 

120.134821 

-0.49446 

112 

2.366E*9 

58.044 

120.182231 

-0.49044 

113 

2.366E-9 

58.685 

120  198239 

-0.48324 

1  14 

2.366E-9 

59.345 

120.197276 

-0.48414 

1  15 

2.357E-9 

60.019 

1  19.951637 

-0.48158 

116 

2.357E-9 

60.7 1 1 } 

1 19.9621 19 

-0.47913 

1  17 

2  350E-9 

61.354} 

119.773243 

-0.48353 

1  18 

2  350E-9 

62.012 

119.785193 

-0.47557 

1  19 

2.347E-9 

62.677 

119.70865 

-0.48729 

120 

2.357E*9 

63.344 

1 19.951269 

-0  488 

121 

2  349E*9 

64.4 13i 

1  19.764163 

-0.46064 

122 

2.355E-9 

65.067j 

1 19.910461 

-0.45743 

123 

2.336E-9 

65.742! 

1 19.429309 

-0  4912 

124 

2  326E-9 

66  4oei 

1  19.170168 

-0.47085 

125 

2.349E-9 

67.052 

1 19.762901 

-0.45371 

126 

2  327E-9 

67.7141 

1 19.204645 

-0.48143 

127 

2.335E-9 

66.369 

1 19.399575 

-0.47939 

128 

2.335E*9 

69.024| 

PKKg&i 

-0.49043 

129 

2.334E*9 

69.69} 

1  19.379054 

■KXEZnil 

130 

23 1 4E*9 

70.367} 

1 18.853069 

-0.48013 

131 

2.323E-9 

71.022! 

1  19.231 173 

-0.48014 

132 

2.335E-9 

7 1 .69 1 1 

1 19.388844 

-0.45177 

133 

2.312E-9 

72.343} 

1 16.804162 

-0.45707 

134 

2.305E*9 

73.0141 

1 16.643191 

-0.4669 

135 

2.304E-9 

73.626} 

-0  47016 

136 

2  307E*9 

74.286 

1 18.681741 

-0.47157 

137 

2  3 1 7E*9 

74.96} 

1 18.932521 

-0.47403 

136 

2.304E-9 

75.6361 

1 18.615157 

-0.47286 

139 

2.307E-9 

76.296! 

1  18.693026 

-0.47248 

140 

2.306E-9 

76.92! 

1  18.6485 

-0.48032 
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Figure  6.1  Polycarbonate  frequency  response  curves  and 
pole  zero  information  for  the  torsional 
mode  of  the  hollow  sample  at  T=21.7°C. 


184 


FREQ  RESP 


Curva  Fit 

KgAftP  Ansi.. 

.2, Tom. 

POLES  6 

ZEROS  B 

1  -1.  40001  *J  xaa.  ess 

-300.  07 A 

2  -2.  73381  3B2.  252 

34..  B517 

3  -S.  1  SSI -4  *J  SOS.  SIS 

353.  353 

4 

7.  55234K 

5 

-04.  3374  — J  536.  34 

S 

B4.  272B  *J  562.  002 

Tim*  delay*>  O.  O  S  Cain— 

-500.  3p  Scale-  1.  O 

Figure  8.2  Polycarbonate  frequency  response  curves  and 
pole  zero  information  for  the  flexural 
mode  of  the  hollow  sample  at  T=21.8#C. 
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APPENDIX  E.  RLC  CIRCUIT  DATA 


CONTENTS: 

-  FREQUENCY  RESPONSE  CURVES  AND  POLE  ZERO  INFORMATION  FOR 
THE  RLC  CIRCUIT  AT  VARIOUS  RESISTANCE  LEVELS  (FIGURES  H.l- 
H.7) 
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Figure  H.l  Frequency  response  curves  and  pole  zero 
information  for  the  RLC  circuit  (R*21H) 
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Flour#  H.2  Frequency  response  curves  and  pole  zero 
information  for  the  KLC  circuit  (R-710) 
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Figure  H.4  Frequency  response  curves  and  pole  zero 
information  for  the  RLC  circuit  (R®17in) 
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Figure  H.5  Frequency  response  curves  and  pole  zero 
information  for  the  RLC  circuit  (R«22in) . 
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Figure  H.6 
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rigor*  H.7  Frequency  response  curves  and  pole  zero 
information  for  the  RLC  circuit  (R*37in) . 
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